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THE INITIAL MAGNETIZATION OF NICKEL UNDER TENSION! 


By B. N. BROCKHOUSE? 


ABSTRACT 

In an attempt to resolve apparent discrepancies in the literature the initial 
susceptibility (x) and Rayleigh constant (a) of pure nickel wire were measured 
between room temperature and the Curie point under various tensions (Z). 
For annealed nickel the quantity xZ is independent of Z, as required by Becker’s 
theory, if the stresses are applied in order of increasing magnitude and if no 
plastic deformation occurs. x and a are found to be connected by Néel’s relation 
a =c x? Jo/Js. The same temperature dependence is shown by xZ as by the 
initial susceptibility of unstressed nickel given in the literature. At temperatures 
below 150°C. the temperature dependence of xZ is consistent with the recent 
theory of Déring. At higher temperatures xZ increases monotonically to the 
Curie point but at a smaller rate than Déring’s low temperature theory predicts. 
This temperature dependence is in apparent conflict with Becker's theory and 
magnetostriction measurements of Déring and Kirkham. It is possible, but 
unlikely, that the conflict arises from the presence of textures in the various 
specimens. An appendix gives the textures of the wire specimens as found by 
neutron diffraction. 


INTRODUCTION 


The interaction of magnetostriction with internal stresses is one of the 
principal factors influencing the magnetization curves of ordinary materials. 
In fact experiments indicate that for a great many materials this is the dominant 
process. Unfortunately no method of determining the internal stresses is known 
and it is not possible to correlate directly experiment and theory. Accordingly 
it is important to determine whether the effect of an external stress upon the 
magnetization is in agreement with theory. 


Becker’s Theory 

A theory of the magnetization curve of nickel under tension was first given 
by Becker (1) in 1930 and is fully discussed by Becker and Déring (2). His 
theory involves three assumptions: 

(1) The magnetostriction is isotropic. 

(2) The applied stress is so great that internal stresses and crystal anisotropy 
can be neglected. 

(3) Only the lowest (nonzero) order terms in the expansion of the free energy 
in the direction cosines of the magnetization are important. 
The magnetization of a nickel wire under tension is 
:' perste HH Z> Z, 
[1] + B= h) BRM 

1 Manuscript received July 16, 1952. 


Contribution from the Department of Physics, University of Toronto, Toronto, Ont. 
2 Now with the Atomic Energy of Canada Ltd., Chalk River, Ontario. 
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where \, is the saturation magnetostriction at T°K., 

J, is the saturation magnetization at T°K., 

Z is the applied (tensional) stress, 

Z;, is the internal stress, 

H is the applied field parallel to Z, 

K is the crystal anisotropy constant. 
According to the theory, therefore, the magnetization curve should be a straight 
line whose slope is inversely proportional to the applied stress and which is 
reversible. Since in practice the internal stresses and crystal anisotropy are not 
completely negligible boundary movements and irreversible rotations also 
occur and the magnetization curve has a small curvature. Because of this the 
initial susceptibility x = Lim. i is the quantity which must be compared with 
experiment. On the Becker theory 


<8 =e 

A careful experimental check of relation [2] has been made for nickel at room 
temperature by Kersten (8), who found the predicted inverse relation between 
x and Z and about the correct value for the constant of proportionality. Since 
the magnetostriction of nickel is not isotropic, A, being — 27 X 10-® in the 
111 directions and — 54 X 10~° in the 100 directions (2) at room temperature, 
Kersten used the mean value of 1/A, in comparing his results with the theory. 

In 1935 Scharff (15) measured x at temperatures from 20°C. to the Curie 
point (358°C.). She found that x was proportional to 1/Z at all temperatures 
but that the proportionality constant was strongly temperature dependent. 
Her values of xZ are shown in Fig. 1. 
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Fic. 1. Comparison of the quantities xZ and J2/3( — 24) as functions of the temperature: 
(1) xZ by Scharff, Reference (15). 
(2) J?/3( — As) by Kirkham (10). 
(3) J2/3( — As) by DGring (4). 
(4) Constant X Jo/Js normalized to (1) and (3) at 0°C. 
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No measurements of \,, the true magnetostriction, have been made at tem- 
peratures other than room temperature but the so-called longitudinal magneto- 
striction (A,’) has been measured from room temperature to the Curie point. 
This quantity is the fractional change in length of a specimen between the 
demagnetized state and saturation. It is proportional to \,, the average of \, 
over the crystal directions, if the domain distribution in the demagnetized state 
is the same for all temperatures and equal to X, if the distribution is isotropic. 
The quantity J,?/3( — \,’) computed from the results of Déring (4) and 
Kirkham (10) is also shown in Fig. 1. It is immediately obvious that equation 
[2] is not obeyed. Equality of the three curves is not to be expected because of 
possible preferred orientations in the specimens but they should have the same 
temperature dependence. DGring’s value of J,?/3( — A,’) at room temperature 
is about that expected for an isotropic specimen while Kirkham’s value suggests 
that his specimen had a preferred orientation in the 100 direction. This preferred 
orientation could also be the reason for the change of slope at about 180°C. 
The crystal anisotropy constant (K) becomes very small at about this tempera- 
ture, below it the distribution of domains in the demagnetized state is controlled 
by the crystal anisotropy and hence by the crystal orientation. Above this 
temperature the domain distribution is controlled by the internal stresses, and 
the two distributions may not be the same. 

No satisfactory explanation for the apparent failure of formula [2] has been 
offered. Becker and Doring (2) suggest that it may be due to \, having different 
temperature dependence for different crystal directions since then 1/X, and 
1/A, would not have the same temperature dependence and hence xZ and 
J ,?/3( — Xs) would also be different. But they consider this explanation unlikely. 
However, the difference between Déring’s and Kirkham’s results lends some 
support to this suggestion, if the preceding discussion is invalid. 


Déring’s Theory 

DG6ring (5) has recently put forward a new, fundamental, theory. He assumes 
that the individual spins are held in position by the exchange forces, the magnetic 
field, and stress anisotropy forces. In a low temperature approximation, and 
assuming isotropic magnetostriction, Déring finds 


[3] xZ = constant X Jo/J; 


where Jo is the saturation magnetization at 0°K. A plot of xZ as given by [3] 
is shown in Fig. 1. The agreement with Scharff’s results is improved but there is 
still a large discrepancy. 

Peppiatt and Brockhouse (14) attempted to check DGring’s theory at low 
temperatures where it is intended to apply. At low temperatures the crystal 


_anisotropy of nickel becomes very large and hence high stresses are necessary 


to make \,Z > K. In order to support the high stresses cold worked wire must 
be used, but these authors found that at room temperature cold worked nickel 
and annealed nickel gave different values for xZ and also showed a different 
temperature dependence. This disagreement was interpreted to mean that 
internal stresses affect xZ for cold worked nickel, and hence that no useful 
results could be obtained at very low temperatures. At temperatures from 
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20°C. to 100°C. xZ for annealed nickel had a significantly smaller temperature 
dependence than that found by Scharff and was in agreement with equation [3], 
but this agreement cannot be considered very significant in view of the small 
temperature range and the resulting uncertainty. 


EXPERIMENT AND CONCLUSIONS 


Apparatus and Method 

The work described in the previous section suggested the desirability of 
making new measurements of xZ up to the Curie point. Specimens of Johnson 
Matthey “specpure”’ nickel wire A, B, and C were prepared as described in 
Table I, in the form of rods 1 mm. in diameter and 250 mm. in length. The 
demagnetizing factor of these rods is sufficiently small to be neglected. 


TABLE I 


MECHANICAL AND HEAT TREATMENT OF SPECIMENS OF HARD-DRAWN NICKEL WIRE (JOHNSON 
MATTHEY “SPECPURE”’ NICKEL). TEMPERATURES ARE APPROXIMATE ONLY. THE TEXTURES OF 
THE SPECIMENS ARE GIVEN IN THE APPENDIX 





Specimen | Heat treatment 
A Annealed 1 hour at 830°C. in He 
B Annealed 2 hours at 1000°C. in He 
c Annealed 1 hour at 700°C., straightened and then 


annealed 1 hour at 830°C. in He 


The specimen was mounted in a vertical furnace and stresses applied by 
hanging weights. The main heater of the furnace was wound on a copper tube 
30 cm. long, and differential heaters at the top and bottom were separately 
controlled by variable transformers. Four thermocouples at different points 
on the specimen indicated the degree of thermal uniformity. When the heaters 
were properly adjusted the maximum measured temperature difference was 
about 2°C. All heaters were bifilar wound to reduce stray magnetic fields. 

The temperature at the center of the specimen was measured by a copper- 
constantan thermocouple which was calibrated at the beginning and end of the 
work. An internal reference temperature was provided by the Curie point 
itself which was found to be very sharp. It could be located to about 2°C. and 
was found to be 358°C., the accepted value, for all applied stresses. 

The magnetization curves were taken by the ballistic ‘‘method of reversals’’. 
The search coil consisted of 1500 turns of copper wire insulated by silicone-glass 
and impregnated with alundum cement. The coil was satisfactory at 370°C. 
though the wire itself is rated at 250°C. The effect of the applied field was 
eliminated by mounting a similar coil parallel to the search coil and connecting 
the two in series opposition. The turns on the compensating coil were adjusted 
to provide complete compensation at room temperature with no specimen 
through the search coil. Measurements above the Curie point indicated that 
the compensation was satisfactory at all temperatures. The deflection of the 
ballistic galvanometer was thus directly proportional to the magnetization J. 

The field coil was such that the applied field was very nearly uniform over the 
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length of the specimen. The vertical component of the earth’s-field was neutral- 
ized by use of an extra winding on the field coil. 

As the furnace contained air it was thought necessary to check the specimens 
for oxidation. One of the specimens was weighed before and after an extended 
series of runs. No change in weight greater than 1 part in 10‘ was detected. 


Determination of xZ 
The principal work was done on specimen A. A stress of 3.82 kgm. per mm.’ 


was applied and the temperature held at 360°C. for several hours. Normal 
magnetization curves were then taken at various temperatures below 360°C. 
with this applied stress. The process was repeated for successively increased 
stresses. Since plastic deformation occurs more readily at the higher tempera- 
tures this procedure ensured that, for any given stress, the specimen was in the 
same metallurgical state at all temperatures. 

Readings were taken for values of the magnetization J up to about 25 gauss. 

Measured values of J/H were plotted against H/ and fitted with a straight line. 
This involves the assumption that Rayleigh’s Law 
[4] J = xH+ aH 
is obeyed for small fields, as it is by nearly all materials. The initial susceptibility 
x is given directly by the intercept on the J/H axis. 

As the galvanometer scale could be read to 0.05 gauss J-values in the range 
5 to 25 gauss were useful for extrapolation to zero field. The uncertainty in x 
from this cause is thus of the order of $%. 

A typical set of results for stresses of 7.60 kgm. per mm.? is shown in Fig. 2. 
It will be seen that the fit by means of Rayleigh’s Law is satisfactory. 











H - QOersteds 


Fic. 2. Susceptibility of nickel wire under a stress of 7.60 kgm. per mm.? at various 
temperatures. 
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Fic. 3. Initial susceptibility (x) of nickel wire (A) under various stresses as a function of 
temperature. In this and similar figures the points at any particular stress were taken in order 
of decreasing temperature. 

Fic. 4. Initial susceptibility (x) of nickel wire (A) under various applied stresses. 


In Fig. 3 x is shown plotted as a function of the temperature for different 
stresses. The different symbols indicate that the measurements were taken on 
different runs. In the case of the small stresses runs were sometimes made after 
exposure to a larger stress. The values obtained were independent of the order 
in which the stresses were applied when the stress was less than 7.60 kgm. per 
mm.” Results for higher stresses are shown in Fig. 4. At these stresses the results 
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Fic. 5. Initial susceptibility plotted as a function of the inverse of the applied stress. 
The slope of the linear part of these curves is the quantity xZ of interest in this work. 
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were dependent on the order in which the stresses were applied. In the curves 
shown in Fig. 4 the maximum stress to which the specimen had been subjected 
after the heat treatment was the applied stress. The two symbols on the 8.85 
and 11.41 kgm. per mm.? curves represent successive runs at the same stress. 
The drooping of the curves in Fig. 3 for low stresses at low temperatures is due 
to the effect of crystal anisotropy, i.e., 4,2 ~ K in this region. 

In Fig. 5 x is plotted against the reciprocal stress 1/Z for a series of tempera- 
tures. Most of the points here plotted were obtained by interpolation from the 
curves of Fig. 3 and Fig. 4. The extent of the validity of the 1/Z law can be 
clearly seen. At low stresses deviations occur because of the internal stresses and, 
at the lower temperatures, crystal anisotropy. The proportionality between 
susceptibility and the reciprocal stress also breaks down at high stresses. Fig. 6 
shows the reduction in area of the specimen as measured after exposure to 
stresses. The breakdown of the 1/Z relation occurs at about the technical yield 
point of the specimen. 

The slopes xZ obtained from curves of the type of Fig. 5 are given in Table II 


0.03 


0.02 


Area Reduction 


0.01 


10 IS 20 
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Fic. 6. Area reduction of 1 mm. diameter nickel wire (A) as function of maximum applied 
stress. 


TABLE II 


THE INITIAL SUSCEPTIBILITY OF NICKEL UNDER TENSION. VALUES OF xZ WERE OBTAINED FROM 
oF X 
CURVES OF THE TYPE OF THOSE IN FIG. 5 








rc. xZ re. xZ ro xZ 

28 22.9 175 25.7 310 34.5 
50 23.3 200 26.7 320 35.7 
75 23.4 225 27.8 325 36.6 
88 23.9 250 29.2 330 37.5 
100 24.1 256 29.5 340 40.2 
125 24.6 275 31.0 345 42.5 
150 25.3 300 33.2 350 46 


354.7 53 
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_ Fic. 7. Open circles are values of xZ found in this work. The dashed curve is Déring’s 
formula [8], the constant having been adjusted for best fit at the low temperatures. The 
values of xZ found by Scharff are plotted for comparison. 


and are plotted as a function of the temperature in Fig. 7. The points are not 
experimental but points for which interpolations were made. Déring’s formula 
[3] gives the dashed line, the constant having been selected to fit the low 
temperature results. In the low temperature region the agreement is good but 
too much significance cannot be attached to this because the temperature 
effect is small in this region. At higher temperatures the experimental curve 
increases with temperature but more slowly than the theoretical curve. This 
increase appears to continue right up to the Curie point. Measurements very 
close to the Curie point (356°C.) were attempted and, though accurate readings 
could not be obtained, it was apparent that the susceptibility was still 
increasing. 

This is contrary to the result obtained by Scharff (Fig. 1) whose xZ vs. T 
curve shows a maximum a few degrees below the Curie temperature. The 
explanation is believed to lie in the relative purity of the specimens. Scharff 
used nickel containing 0.6% manganese. According to Marian (12) the Curie 
point of this alloy is about 350°C. instead of 358°C. for pure nickel. In general 
the Curie points of alloys are difficult to measure, being ‘‘smeared out’’, possibly 
because of localized nonuniformities. A spread in the local Curie temperatures 
around a mean of 350°C. could account for the observations of Scharff. 

In the lower temperature region there is also a large discrepancy between the 
curve of Fig. 7 and Scharff’s results, also shown in the figure. Normalization of 
the two curves does not make them agree. It is thought that part of this difference 
arises in the following way. The magnetization curves obtained by Scharff 
consisted of three or four points only, taken in the range of magnetization zero 
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to 25 gauss, and fitted by a law of the form J = xH. Thus any curvature in the 
magnetization curve was neglected. In every case examined by Peppiatt (14) 
and the author the H? term in [4] was appreciable, as exemplified in Fig. 2. 
In general the constant a@ in equation [4] increases with temperature at a faster 
rate than x itself. Assumption of a linear magnetization curve results in an 
erroneously high temperature variation of x. 

In order to determine, as far as possible, the extent to which the xZ vs. 
temperature curve of Fig. 7 is characteristic of nickel under tension as distinct 
from the particular specimen used, several experiments were performed. These 
subsidiary experiments are described in the two sections following. 


Effect of Cold Working on xZ 
Measurements were made at lower stresses after application of a higher 
seasoning stress. Fig. 8 shows one such set of measurements taken at various 


a- 10.4 kgm/mm? 





20 10 8 6 5 


100 200 300 0, 
Temperature °C. Z- (Reciprocal Scale) -kgm/mm! 


Fic. 8. Initial susceptibility of a nickel wire (A) at various stresses, after application of a 
seasoning stress of 19.5 kgm. per mm.? with resulting cold working. 

Fic. 9. Initial susceptibility x plotted against the reciprocal stress 1/Z, showing effects 
of cold working. The arrows indicate the direction of change in applied stress. 


stresses after application of a stress of 19.5 kgm. per mm.” This caused an area 
reduction of 3.4%. These values of x are not proportional to 1/Z as is shown 
in Fig. 9 for two temperatures. On these curves the arrows indicate the direction 
of change in applied stress. Though x is not proportional to 1/Z for cold-worked 
wire it appears possible that the curves are asymptotic to the curves obtained for 
annealed wire. 
Effects of Annealing Procedure 

Measurements were made on another specimen (B) which had been annealed 
at 1000°C. These results are shown in Fig. 10. The 1/Z relation is not obeyed, 
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Fic. 10. Initial susceptibility of nickel wire (B) annealed two hours at 1000°C. under 
various applied stresses. 


the initial susceptibility at the lower stresses being slightly too high. The 11.4 
kgm. per mm.? curve is identical with that for specimen A for the same stress, 
shown in Fig. 4. At this stress a small amount of cold working (0.32% area 
reduction) occurred. It is believed that at the lower stresses boundary movements 
made an appreciable contribution resulting in a higher susceptibility than expect- 
ed and departure from the 1/Z law. This is supported by the observation that 
magnetization curves for these stresses had more curvature than usual. The 
slight cold working at 11.4 kgm. per mm.* presumably suppressed the boundary 
movements. 

The situation is complicated by the fact that metallurgical changes ‘in the 
specimen can occur during the course of an experiment at temperatures below 
the Curie point. For certain purposes nickel is annealed at about 300°C. which 
results in small but definite changes in the mechanical properties. This is spoken 
of as ‘‘stress relief’’ or perhaps more properly “‘stress equalization”’. 

A specimen (C), which had been annealed at the same temperature as A 
(but which had had different treatment prior to annealing as indicated in 
Table I), was studied under an applied stress of 6.37 kgm. per mm.? at a 
temperature of 150°C. where the effects of crystal anisotropy are smaller than 
at room temperature. The curve of J/H vs. H is shown as curve I in Fig. 11. 
The specimen was then annealed for eight hours at 330°C. without removing the 
stress and the initial magnetization again measured at 150°C. (Curve II). 
The initial susceptibility is quite different from previously and is the same as 
that given in Fig. 3 for specimen A at the same stress. This experiment illustrates 
the necessity for making the measurements in a definite order of temperatures 
and stresses, and justifies the order of temperatures adopted. 
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Fic. 11. Susceptibility of nickel wire (C) at 150°C. under stress of 6.37 kgm. per mm.? 
I-Annealed without stress as shown in Table I. II-Then annealed eight hours at 330° under 
stress of 6.37 kgm. per mm.? 


Effects of Internal Stresses 

Thiessen (16) found that nickel wires under stresses up to several kgm. per 
mm.’ had susceptibilities much smaller than would be expected from equation 
[2] and that the susceptibility initially increased with applied stress. Thiessen 
found that repeated annealing eliminated the effect and he interpreted his 
results in terms of an anisotropic distribution of internal stresses. 

No such effect was shown by the specimens used in the present work, but 
Peppiatt and Brockhouse (unpublished) observed that the susceptibility of 
specimens annealed at about 1100°C. did increase initially with applied stress. 
The slight cold working produced by a stress of 11.4 kgm. per mm.? was sufficient 
to destroy the effect, which in these specimens was equally prominent at 
temperatures from 20°C. to 200°C. This temperature independence indicated 
that crystal anisotropy was not involved and that the effect was caused by 
internal stresses. The stresses deduced from equation [2] had values of about 
8 kgm. per mm.’ and stresses of this order may thus be present in apparently 
well-annealed nickel. 

In general internal stresses can affect the results in two ways (a) by causing 
deviation of the domain vectors from directions normal to the axis of the wire, 
hence permitting boundary movements to occur, and (d) by changing the mean 
stress over the volume of a domain, hence affecting the rotational susceptibility. 
A change in the mean stress will occur only if the internal stresses are coherent 
over volumes of the order of a domain volume and boundary movements can be 
appreciable only if the gradient of the internal stress is not too large. Both 
mechanisms thus require that the stress should be of long “‘period”’, i.e. a slowly 
varying function of position, as well as of large magnitude, in order to affect 
the results. It seems likely that annealing may increase the “period” of the 
internal stresses as well as decrease their magnitude, in the same way as it 
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promotes crystal growth. Hence it is not too surprising that specimen (A), 
annealed at a comparatively low temperature, seemed to show the most freedom 
from effects due to internal stresses, and that slight cold working eliminated 
most of the differences between specimens. In view of these considerations, the 
author does not believe that internal stresses appreciably affect the results 
given in Fig. 7 and Table II. 
The Rayleigh Constant 

The coefficient of the quadratic term in equation [4], known as the Rayleigh 
constant, was also computed from curves of the type of Fig. 2. One set of values 
is shown in Fig. 12. According to Néel (13) the quadratic term arises from 
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Fic. 12. Rayleigh constant of nickel wire A under stress of 8.85 kgm. per mm.? The 
dashed curve is the best fit to equation [6]. 


irreversible boundary movements. It can be shown, simply and generally, that 
if the number of domain boundaries does not change with temperature and if the 
irreversible movements take place between definite positions in the crystal, 
which also do not change with temperature, then 
~ Jo 2 
[5] a@ = constant X > xp 

Js 
where xg is the susceptibility due to boundary movements. This equation is 
given by Néel (13) but is not dependent on the details of his theory. The points 
in Fig. 12 are fitted with a curve of the form 


9 


. J 
[6] a = constant X + x. 


Within the rather poor accuracy the curves can be satisfactorily fitted though 
for some a better fit is obtained with an index slightly higher than 2. In general 
the deviation from the index 2 was greatest for those curves which failed to 
obey the 1/Z law. Since the initial susceptibility x is believed to be due almost 
entirely to a rotation process, this implies that the rotation and boundary 
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movement processes probably have the same temperature dependence, and 
hence probably arise from the same mechanism—the interaction of magneto- 
striction with stress. 
It is perhaps worth mentioning that over a considerable temperature range 
equation [6] turns out to be nearly 
a = constant X x‘, 


the relation between a and x found empirically by Radanovic and Weiss and 
de Freudenreich (17) for unstressed nickel. 


GENERAL DISCUSSION 
Initial Susceptibility of Unstressed Nickel 
A large number of investigations of the initial magnetization of nickel have 
been made (7, 10, 16, 17). The collected results for the initial susceptibility are 
shown in Fig. 13. The temperature scale may be divided into two regions 
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Fic. 13. Initial susceptibility of nickel by various authors. (a) Kirkham (10); (6) and 
(c) Kahan (7); (d) and (e) Thiessen (16); (f) Weiss and de Freudenreich (17); (g) and (4) author 
(unpublished). Generally speaking the initial susceptibility increases with annealing tempera- 
ture. 


. according to the sign of the crystal anisotropy constant K. In the upper region 
the crystal anisotropy is usually negligible, i.e. K « X, Z;. Where this condition 
is fulfilled the temperature dependence of x is very nearly the same as that 
found for nickel under tension in this investigation. Curves d, e, f in this region 
have the same temperature dependence as (xZ) to within about 5% (see Fig. 7). 
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Fic. 14. Plot of xZ found in this work against initial susceptibility at the same temperature 
given in curves of Fig. 13. In the high temperature region proportionality holds between all 
the curves, within the accuracy of the data, except in the immediate vicinity of the Curie point. 


This relationship is displayed in Fig. 14 in which (xZ) is plotted against values 
of x taken from the curves in Fig. 13 at the same temperature. 

The marked decrease in x at lower temperature is believed by the author to 
arise from a mixing of rotation and boundary movements caused by demagnetiz- 
ing fields produced at crystal and domain boundaries. Where crystal anisotropy 
is the dominant energy term (K > XA, Z;) the susceptibility for rotation 



























2 

Xr = constant X x 
This is usually only a small part of the total susceptibility, which comes 
principally from boundary movements. But the boundary movements are 
inhibited by demagnetizing fields produced at crystal and domain boundaries 
and at gross imperfections. These demagnetizing fields can be markedly reduced 
by the alternative rotation process and hence this has a profound effect on the 
over-all susceptibility. A detailed account of these ideas will be presented in 
another paper. 
If any such theory of the falling off of the susceptibility at low temperatures 
is accepted then the variation of the high temperature susceptibility can 
properly be investigated only with specimens which show no anomaly where K 
passes through zero (i.e. somewhere in the vicinity of 450°K.),* since only for these | 
specimens will the condition K KX, Z; be satisfied at higher temperatures. This 
excludes curves a, 6, c. The remainder of the curves have, as noted, the same 

temperature dependence as stressed nickel, shown in Fig. 7. 

*The measurements of Honda, Masumoto, and Shirakawa (6) show that above about 375°K. 
the (1, 1, 1) directions are no longer the easy directions of magnetization. It is believed (10, 3) that 
the anisotropy constant K changes sign somewhere above this temperature, but only the roughest 


estimate of K can be obtained, particularly in view of the fact that anisotropy constants of higher 
order and internal stresses may then be important. 
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Other Magnetic Quantities 
In addition to the initial susceptibility, expressions for other quantities have 
been derived from the strain theory, especially by Kersten (9). These include 
(a) Alteration of remanent magnetization of a wire by an applied stress 


(é22) ae 8 
dZ z=0 7 4Z; s 


(6) Alteration of magneto-resistance by an applied stress 


- (22) ik ee 
OZ \ p 15 Z; 
where D is found experimentally to be proportional to J,’. 
(c) Change in Young’s modulus E between the demagnetized state and 


saturation (AEF effect) 
1 2X; 
(5) " eee 


Where the crystal anisotropy K can be neglected the measurements of Thiessen 
(16) and Késter (11) are in agreement with theory for all three effects, if Z; 
is assumed to be independent of temperature. 

It will be noted that all three effects arise from competition between external 
and internal stresses, while the initial susceptibility and Rayleigh constant 
involve competition between internal or external stress and a magnetic field. 
That is, the theory describes correctly the temperature dependence of processes 
which involve only stresses but not processes which involve both stress and 
fields. However existing information is consistent with the view that the latter 
processes have always the same temperature dependence. This is evidence for a 
fundamental origin of the failure of the Becker formula [2] but unfortunately 
not conclusive evidence. A study of the true magnetostriction (and perhaps xZ) 
as a function of temperature and crystal orientation is needed. 

The author wishes to thank Prof. H. Grayson Smith of the University of 
Alberta for helpful correspondence. 

APPENDIX 
Texture Determination in Annealed Nickel Wires by Neutron Diffraction 

The texture of cold worked nickel wire resembles that of copper, having a 
strong [111] texture and a weaker [100] texture, i.e. the normals to the (111) 
and (100) planes respectively tend to lie in the direction of the axis of the wire. 
The texture of annealed nickel wire has not, to the author’s knowledge, been 
reported. By neutron diffraction studies with a crystal spectrometer, the tex- 
tures in the specimens A and B used in this work have been determined. The 
results are for wires which were annealed as stated in Table I and then subjected 
to stresses as described in the text. It is believed that the stresses applied after 
annealing did not affect the texture appreciably. A specimen which had been 
annealed at about 1100°C. under the same conditions as the specimen used by 
Peppiatt and Brockhouse (14) was also tested. This specimen had not been 
stressed beyond the elastic limit. 
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Neutron diffraction provides an easy method of determining textures in 
annealed materials. The large areas of the beams used in neutron spectrometers 
and the great depth of penetration of the neutrons ensure that the pattern 
obtained applies to a large volume of the specimen and hence averages over a 
large number of crystals, even if the crystals are fairly big. Because of the large 
depth of penetration surface conditions are unimportant, in contrast to texture 


determination with X rays. 
To achieve adequate counting rate the specimen wires were cut into short 
lengths and secured to thin aluminum foils as shown in Fig. Al. The wires 
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Fic. Al. (a) Form of specimen. Nickel wire segments are glued to a thin aluminum sheet. 

(b) Experimental arrangement. The counter rotates in plane of paper about O. 

The specimen is mounted in plane A-A, perpendicular to paper, and is rotated about the 
normal OB. 


formed a circular array with diameter smaller than the neutron beam. With 
this construction no corrections for geometry or absorption were required. 
The array was located in the plane A-A which normally bisected the angle 
between the incident and emergent beams but could be rotated out of this 
position. The spectrometer was set on the appropriate Debye Scherrer line 
(h, k, 1) and the array was rotated in the plane A-A about the normal OB. The 
intensity of the line was recorded as a function of the angle which the axes of 
the wires made with the spectrometer plane. This intensity is proportional, 
neglecting extinction, to the volume of material with planes (h, k, /) perpendicular 
to the bisector of the angle between the incident and the emergent beams 
(in Fig. Al line OA). Since the wire has cylindrical symmetry this intensity 
gives directly the distribution of the planes with respect to the axis of the wire. 
The distributions for the (111), (200), and (220) planes are shown in Fig. A2. 
The texture may be described as a [100] texture (25 + 10%) and a texture 
(65 + 10%) in which a 111 axis makes a small angle with the axis of the wire. 
This angle appears to increase with annealing temperature. The texture of cold 
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Fic. A2. Normalized intensities of 111, 200, and 220 lines of specimens annealed at 830°C. 
(A), 1000°C. (B), and 1100°C. as a function of the angle between the wire axis and the spectro- 
meter plane (plane of the paper in Fig. A1). The crosses represent intensities measured with 
the axes of the wires in the spectrometer plane as a function of the angle between the plane 
and the “‘half-angle”’ position of the specimen. 


worked wire (from another source however) was verified to be [111] with a 
weaker [100] texture. Further work is being done on textures of annealed mate- 
rials and will be reported in a separate publication. 

The work described in this appendix was carried out at the Chalk River 
Laboratory of Atomic Energy of Canada, Ltd. 
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COUNTING RATE MEASUREMENT BY INTERVAL SAMPLING! 









By M. F. Croucn Anp G. E. HAyNAM 


ABSTRACT 


Two methods are described for measuring high counting rates with a counting 
device having a long resolving time. Both methods are based on the deliberate 
introduction of a large but calculable counting loss to transform a high counting 
rate to a lower rate. Formulas giving the counting rate in terms of the reduced 
rate are presented, together with formulas for the standard deviation of the 
counting rate. The performance of circuits used to study one of the methods is 
described, and a comparison is made with conventional rate measurement 
methods. It is found that, while the relative simplicity of equipment is gained 
at the expense of speed of measurement, the latter sacrifice is not as severe as 
might be expected. 















INTRODUCTION 

In addition to the two conventional methods for measuring high counting 
rates with a Geiger—Miiller counter, viz. by means of scaling circuits and by the 
use of counting rate meters, a third method, which might be called ‘interval 
sampling”, can, in principle, be used. This method makes use of the fact that 
the distribution function (Equation [1], Fig. 1) for the time intervals between 
successive G-M counter pulses arriving at a mean rate R, 


[1] Po(t) = Rexp (— Rd), 



















is a known function possessing a ‘‘tail’’ at its upper end. By measuring the rate 
of arrival S of intervals longer than some predetermined value 7, corresponding 






R(t) 


PROBABILITY 


TIME INTERVAL %t% 


Fic. 1. Distribution function for time intervals between G-M counter pulses. 


to the shaded “tail” of Fig. 1, the counting rate R can be found by using the 
implicit relation Equation [2]. 
[2] S = Rexp ( — RT) 

1 Manuscript received September 26, 1952. 


Contribution from the Department of Physics, Case Institute of Technology, Cleveland, 
Ohio, U.S.A. 
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However large R may be, a value of T can always be found to make S 
sufficiently low to be measured with an electromechanical register. Thus, 
whereas the scaler passes on 1/64 or 1/100 of the G-M counter pulses to the 
message register, in a regular fashion, the interval sampling circuit passes on a 
fraction S/R = exp ( — RT) of the G-M counter pulses, in random fashion. 
It might therefore appropriately be designated a “‘pseudoscaler”’. 

Although Equation [2], regarded as a counting loss formula, has been verified 
(5, 6) up to a scaling factor R/S = 25, no previous attempt seems to have been 
made to utilize Equation [2] as the basis of counting rate measurements, which 
would require values of R/S ( = exp (RT)) up to 300 and larger. Although the 
apparent lack of any previous studies of this method is possibly attributable to 
several inherent disadvantages, it was nevertheless felt worth while to investi- 
gate certain aspects of the method. 

The investigation to be described had four objectives, which were (a), to 
ascertain the value or values of 7 needed to accommodate G-M counter rates 
from 0 to 50,000 counts per min., (0), to design and test circuits for performing 
the operation of interval sampling, such circuits to have 7-values readily calcul- 
able from the circuit constants, (c), to investigate how the time required for a 
counting rate determination by the interval sampling method varies with the 
counting rate, the accuracy desired, and the accuracy of time discrimination, 
and (d), to determine the advantages and disadvantages of this method over 
conventional methods. 

A related method for measuring rates with long resolving time equipment, 
which suggests itself when ways are considered to utilize the register at higher 
rates, will also be briefly discussed. This method is basically different, but it is 
equally well described as ‘interval sampling’. In this method, the counting 
device is rendered sensitive every 7 seconds for a brief interval of time 7. If the 
register records a single count whenever one or more input counts are received 
during the sensitive interval 7, then the relation between input counting rate 
R and output counting rate S is 
[3] S = (1/r)(1 — exp ( — RT)). 


This method will be discussed briefly in the following section, but no analysis 
of actual circuits will be presented. 
THEORY 
Equation [2], giving the relation between rates R and S for the first method, 
has been derived by several authors (1, 5). In addition, Feller (1) has worked 
out the rigorous expression for the standard deviation of the rate S. His result, 
to a sufficient approximation for the present purpose, is given in Equation [4]. 


[4] o, = (n,(1 — 2ST))} 


(n, is the total number of sampled counts). When Equation [2] is regarded as 
an implicit relation between R, S, and 7, the relation between the errors of 


these quantities is computed to be (7) 


[5] o/R = (1/1 — RT)V (Ci — 2ST)/n,) + (RT)*(07/T)?. 
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Although the error in R is seen to be a function of R, S, T, n;, and o7/T, each 
value of R in normal operation is associated with unique values of S and 7. 
Therefore by assuming a reasonable value for o7/7T, a one-to-one correspon- 
dence between n, and R exists for any desired accuracy op/R of rate measure- 
ment. After the circuits that were tested have been described, Equation [5] 
will be used to calculate typical values for the time required to determine R to 
a given accuracy. 

Equation [8], for the alternative method described in the introduction, is 
derived as follows: For G-M counter pulses occurring at a rate R and having 
a Poisson distribution, the probability of no count occurring in a given time 
interval T is easily found to be g = exp ( — RT). Therefore the probability of 
one or more G-M counts in a given time interval T is p = 1 — exp( — RT). 
Since g and p also represent respectively the probability of zero and one count 
recorded on the message register, the frequency distribution of the number of 
counts recorded in a fixed period of time is a binomial distribution. The standard 
formula (Fry (2)) for the mean of the distribution leads at once to [3]. Similarly 
the standard deviation formula becomes, in the present case, 


[6] o,=V 1, exp ( — RT). 


In [6] the coefficient of exp( — RT) should be the mean expected number of 
recorded counts. Ordinarily it is a reasonable approximation to use the actual 
number 7,. 

All of the expressions given here for rates and standard deviations are 
asymptotic expressions valid only for measurement times large compared to the 
critical times 7 and 7. 

EXPERIMENTAL METHOD 

As a working basis, maximum and minimum values for S of 10 counts per 
sec. and 2.5 counts per sec. were assumed —the former being the approximate 
upper limit of the Cenco dial-type message register and circuit used, for random 
pulses, and the latter being an arbitrarily imposed lower limit on the rate of 
data accumulation. The first circuit studied, designated ISC-1, has been 
described in a preliminary form (4). It contains a scale-of-two to permit analysis 
of every other interval, so that S and R are related as in Equation [7]. 


[7] S = 3 Rexp (— RT) 


With this arrangement, four values of 7 are sufficient to cover the range of 
counting rates desired. Fig. 2 is a plot of the family of curves (Equation [7]) 
of S vs. R for the four chosen values of 7. The range of counting rates over 
which each value of 7 is used is indicated by the heavy portion of each curve. 
Fig. 3 is a plot of the counting rate versus the scaling factor R/S plotted on 
a logarithmic scale. Plotted in this manner, Equation [7] becomes a straight 
line with slope 1/7, R/S intercept 2. Fig. 4 gives the same straight-line plot 
for a circuit described by Equation [2], five values of T being needed to cover 
the desired range of values of R. 
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Fic. 4. Input counting rate vs. scaling factor, ISC-2. 


The ISC-1 circuit (Fig. 5) was a straightforward approach to the problem, 
with little regard for economy of tubes. In brief, an exponential saw-tooth 
generator is turned on for the duration of each interval studied, so that the 
length of the time interval is indicated by the final voltage amplitude of the 
saw-tooth. These saw-tooth signals are then applied to a series diode amplitude 
discriminator, the cathode of which is maintained, in the absence of signals, at a 
reference potential corresponding to the critical time desired. An additional 
diode is used to clamp the cathode of the discriminator diode after it has risen 
3 volts above the reference potential, so that this potential cannot stray far 
from its proper value. The discriminator output is amplified and used to trigger 
a multivibrator message register driver. 

The saw-tooth generator used is a standard sweep generator using a ‘“‘hard”’ 
tube, whose plate potential rises exponentially with a time constant RC when 
the grid is driven below cutoff, but returns to its starting value very quickly 
when the grid signal is removed because the capacitor C is recharged by current 
flowing through the tube. It is to be noted that the nonlinear character of the 
saw-tooth doesn’t introduce any difficulty in this application. 

The second circuit tested (miniaturized), designated ISC-2 (Fig. 6), represents 
an attempt to perform the interval sampling operation with a near-minimum of 
vacuum tubes. Here the same saw-tooth generator is used, except that it is 
normally cut off in this case, and the grid signal consists of short positive pulses 
applied at the start of each interval. 
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Fic. 5. Block diagram, wave forms, and partial schematic circuit diagram, ISC-1. 
RESULTS 

The curves characterizing ISC-1 were verified along their entire length 
(circular data points, Fig. 3) using an auxiliary scaling circuit to measure the 
rate S. The complete ISC-1 circuit including message register was then checked 
(square data points) near the ends of the range of R-values covered with each 
value of 7. The agreement between experimental results and the theoretical 
curves was felt to be satisfactory, since, with the exception of the timing capaci- 
tors, the circuit was assembled from ordinary circuit components such 
carbon resistors with + 10% tolerance. 

The performance of ISC-2 was verified only for T = 0.04 sec., since exper- 
ience had shown it to be most difficult to obtain proper operation for the larger 


' T-values. Fig. 7 is an oscilloscope photograph of the wave form at the saw-tooth 


generator. The plate potential drops so rapidly at the start of each interval that 
the camera doesn’t record this part of the wave form. The discriminator output 
is used as a blanking signal, so that a gap in the trace should (and does) occur 
each time there is an interval long enough for the associated saw-tooth to rise 
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Fic. 6. Block diagram, wave forms, and schematic circuit diagram, ISC-2. 









to the discriminator reference potential. The clamping action which takes 
place as soon as the discriminator pulse is generated is also apparent from the 
photograph, the apparent slight drift after clamping being an oscilloscope 
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Fic. 7. Oscillogram of saw-tooth generator output voltage vs. time, ISC-2. 


The relationship between circuit parameters and the 7-value of the circuit 
is readily found for both circuits to be 


[8] exp (T/RC) = (Epp — €n0) (Esp — Ca), 


where E,, = saw-tooth generator plate supply potential, e,9 = saw-tooth 
generator plate potential (no signal), and eg = discriminator diode cathode 
potential. This expression is found to give 7-values agreeing closely with those 
inferred from the slope of the curves through the experimental points of Figs. 
3 and 4. As an example, for circuit ISC-1 operated with a nominal 7-value of 
0.006 sec., the circuit constants were measured as R = 0.68 MQ, C.= 0.015 uf., 
Ey = 110 v., e50 = 1.7 v., and eg = 48.0 v. Using these values, Equation [8] 
gives T = 0.057 sec. Thus the requirement that 7 be calculable from the 
circuit parameters is reasonably well satisfied. 
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Since the Schmitt circuit of ISC-1 responds readily to sinusoidal inputs, a 
simple alternative procedure for setting the discriminator level to correspond 
to a desired value of JT can be employed using an audio-frequency oscillator 
connected to the input. As the frequency is raised, a point is reached at which 
the message register suddenly stops counting. Calling this frequency f,, the 
T-value is given by T = 1/f,. Alternatively, the discriminator can be set by 
connecting the a-c. filament voltage to the input and using the appropriate 
7-value. 

DISCUSSION 

Having shown that circuits can be made to perform the operation of interval 
sampling (first method), it remains to consider the effect on speed of measure- 
ment of counting only the long intervals. As an example, Table I lists (for 
circuit ISC-2) the number of sampled counts m, and the average time ¢ to 
register these counts in order to give a standard deviation of 5% in the inferred 
value for R (i.e., ¢g/R = 0.05). It was assumed that o7/T = 0.02, requiring, 

TABLE I 
NUMBER OF SAMPLED COUNTS NEEDED TO DETERMINE R To 5% accuracy (ISC-2) 


No. of Avg. time | Conventional 

Input rate | r sampled | of measurement measurement 

R (counts/sec.) | (sec.) counts Ms (sec. ) time (sec.) 
833 .007 21.6 8.9 0.5 
700 .007 32.6 6.3 0.6 
540 007 61.3 §.2 0.7 
540 010 25.8 10.4 0.7 
460 | 010 39.0 | 8.2 0.9 
360 .O10 68.6 7.0 1.1 
360 014 30.4 13.1 1.1 
300 O14 47.2 10.5 1.3 
220 O14 102 10.1 1.8 
220 021 37.5 is 1.8 
160 .021 81.5 14.7 2.5 
110 021 252 | 23.0 3.6 
110 .035 58.0 24.8 3.6 
80 035 133 27.3 5.0 
60 035 385 52.4 6.7 


according to Equation [8], that the threshold potential be held to about 13%. 
Fig. 8 is a graph of time of measurement vs. rate. 

The most striking feature of the ISC method is apparent from Table I, viz. 
that, despite the randomness of arrival of the sampled pulses, the counting rate 
can be determined to within 5% by a total count of as little as 22 counts, 
instead of the usual 400. This surprising result is due to the fact that, except 
for a range of very low rates, a given percentage variation of R-values corre- 
sponds to a much larger range of S-values, percentagewise. That is, there is a 


‘sort of ‘‘band spread” effect, to borrow a term from radio communication. 


As a result of this effect, the measurement time is < 15 sec. for 220 < R < 833 
per sec., and it is < 50 sec. for 65 < R < 833 per sec. 

However the region of the maximum of the S vs. R curve (Fig. 2) is readily 
seen to be characterized by inherently low accuracy, and the limit of accuracy 
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Fic. 8. Measurement time for 5% error in R, ISC-2. 






of measurement depreciates to + 30% in the worst case. It is, of course, possible 
to remedy this difficulty by providing additional T-values, avoiding the use of 
this region of any of the curves. 

In evaluating the possible utility of the interval sampling method, one observes 







the following advantages: 

1. Simplicity (cf. Fig. 6). The improvement over conventional scaler circuits 
is quite considerable. Counting rate meters, although simple in principle, are 
usually somewhat more complex. 

2. Nonsusceptibility to multiple G-M counts, which provide the nonconstant 
contribution to the slope of the G-M “plateau”. 

3. Adaptability to the preset time method of rate measurement. For a wide 
range of rates (cf. Fig. 8), a fixed duration of measurement yields roughly a 
constant precision of measurement. 

4. The ‘“‘band spread” effect noted above leads to unexpectedly short times 
of measurement, although the speed of measurement is still far less than with 
conventional methods. 

5. All rates from 0-50,000 per min. and greater can be counted. This consti- 
tutes an advantage over the counting rate meter, which is not adaptable to very 
low rates, such as one encounters in cosmic ray measurements, etc. 

6. Poisson distribution histograms, for student laboratory experiments, are 
readily obtained without additional equipment. 

To summarize the disadvantages, they are: 

1. Speed of measurement is sacrificed, although the sacrifice is not so drastic 
as might be expected. (Cf. Table I.) 

2. Only input pulses obeying a Poisson distribution can be counted. 

3. Accuracy is limited by the precision with which the time discrimination 




























can be accomplished. 
4. The curve of S vs. R is double valued (cf. Fig. 2), although one can easily 
devise simple operational procedures for eliminating this ambiguity. 
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5. The circuit gives only rates, and cannot give number of counts. 
6. A calibration curve is required. 
7. The proper 7-value must be selected to give 2.56 < S < 10 counts per sec. 

It is acknowledged that the advantages of this method of rate measurement 
are probably not sufficiently compelling to make it a serious competitor to the 
scaler or counting rate meter at present, except possibly for low rate measure- 
ments. Hence the results are presented primarily for their own intrinsic interest, 
with the added thought that the results may be of practical use in some future 
rate measurement problem which may be especially amenable to solution by 
the interval sampling method. Faster message registers now becoming available 
may permit ISC-type circuits with improved characteristics. 

In this connection, the lack of a practical motivation means that the circuits 
described are far from being engineered circuits. Should an actual application 
arise, more elegant and reliable circuits, or possibly mechanical devices for 
accomplishing the same end, could undoubtedly be developed. 

Although the second method has not been studied by analysis of actual 
circuits, it is apparent that it does have the advantage of a single-valued relation 
(Equation [3]) between R and S. Moreover, using a reasonable value of 7, say 
1/30 to 1/60 sec., the entire useful range of input counting rates can be covered 
with only two different 7-values. The actual circuit to accomplish this type of 
counting could be quite simple—for example an a-c. operated phase-controlled 
shield grid thyratron circuit would exhibit the desired periodic intermittent 
sensitivity. It is also apparent that this method can be used for rate measure- 
ment of incoherent periodic phenomena. 

Recently Greenberg and Happ (3) have described another instrument for 
studying the statistical properties of counter pulses. With their method the mean 
deviation is continuously measured, in conjunction with a conventional counting 
rate measurement. 


The authors wish to express thanks to Prof. A. L. Hughes, Dr. O. S. Duffen- 
dack, Prof. Jonathan Townsend, and Dr. Harold Lyons for helpful criticisms. 
The helpful suggestions of the referee, particularly that pointing out an error 
in Equation [4] as originally written, are gratefully acknowledged. 
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PLAFONNEMENT DES AEROSTATS A GRANDE ALTITUDE! 
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SOM MAIRE 


On décrit un procédé simple combinant l’usage de ballons fermés de caoutchouc 
et d’un ballon ouvert de taille moyenne en Polythéne, qui a permis de faire 
plafonner pendant 9 a 10 heures, des émulsions sensibles aux rayons cosmiques a 
des altitudes comprises entre 24 et 32 km. 











1, INTRODUCTION 





L’étude des rayons cosmiques avec des ballons requiert de préférence une as- 
cension rapide, un plafonnement prolongé et une descente rapide. Ces conditions 
ne sont pas faciles a réaliser par l’emploi des seuls ballons de caoutchouc fermés a 






moins de mécanismes lance-ballasts. Ils gardent en effet leur force ascensionnelle 
et continuent de monter jusqu’a ce qu’ils crévent. Le plafonnement est impossible 
avec un seul ballon qui suffirait pourtant 4 porter une charge intéressante pour 
l’expérimentateur. Théoriquement un systéme de trois ballons ou davantage, 








gonflés de fagon a régler trés exactement les hauteurs d’éclatement et les vitesses 






ascensionnelles, assure le résultat cherché, mais ce réglage est difficile et 
aléatoire vu les propriétés variables des ballons. Des systémes de nombreux 
ballons avec une forte charge sont plus faciles 4 régler, et donnent parfois un bon 
plafonnement, mais les renseignements publiés font défaut 4 ce sujet. En somme 








avec des ballons de caoutchouc fermés, une faible charge peut bien étre portée a 






une grande altitude, mais elle n’y plafonne pas facilement. 
D’autre part, nous avons décrit un nouveau procédé de fabrication des ballons 
. de Polythéne ouverts (3). Les ballons de ce type sont capables de plafonner. Ce- 
pendant avec les matériaux actuellement disponibles pour leur construction, le 








plafonnement au-dessus de 25 km. exige des ballons vastes et lourds, dont la 






manipulation est peu commode. Ces ballons n’apportent pas une solution au 






probléme du plafonnement de petites charges 4 de plus grandes altitudes. 




















2. PROCEDE NOUVEAU 
Au cours de |’été 1952, nous avons essayé une solution qui s’est avérée excel- 
lente, et qui apparemment n’a pas été décrite ailleurs. Elle combine des ballons 
de caoutchouc fermés, de force ascensionnelle constante, et un ballon de Poly- 
théne ouvert, de taille moyenne, qui perd 4 mesure qu’il monte, de sa force 
ascensionnelle. Le procédé est équivalent 4 une ascension de ballons de caout- 


chouc avec un lance-ballast négatif, ou avec une charge croissant avec I’altitude. | 
Ou encore, a une ascension d’un ballon de Polythéne avec suppression constante 


d’une fraction choisie de son poids total et de sa charge. 

Il faut naturellement prévoir au départ la force ascensionnelle nette F, des 
ballons de caoutchouc pour qu’elle ne suffise pas 4 supporter le poids P du reste 
de la charge y compris le ballon de Polythéne supposé vide, sinon I’ascension 
ralentit une fois que celui-ci a pris forme, mais ne s’arréte pas avant qu’un ballon 

1 Manuscrit regu le 24 septembre, 1952. 
Contribution de I’ Institut de Physique, Université de Montréal, Montreal, Que. 
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de caoutchouc ait crevé. Au moment du plafonnement, il faut que le ballon de 
Polythéne ait 4 contribuer une force ascensionnelle F, positive: 


F, = P — F,. 


F, est proportionnelle a la pression au moment du plafonnement. I] faut en outre 
que l’altitude du plafonnement soit inférieure a l’altitude d’éclatement des 
ballons de caoutchouc dans |’état de gonflement qu’on leur a imposé. Pour un 
plafonnement a I’altitude maximum, il faut les gonfler peu, en mettre plusieurs 
et avoir une faible valeur de F,. Le lancement du systéme peut se faire 4 une 
grande vitesse ascensionnelle, le ballon de Polythéne pouvant étre fortement 
gonflé puisque F, ne dépend pas de son gonflement initial. 

Les ballons de Caoutchouc J8-800 Darex chauffés dans l’eau 4 60°C. pendant 
cing minutes ont été gonflés a I'hydrogéne jusqu’a ce qu’ ils tirent a leur base avec 


loo 
pieds 


305 
m. 


Fic. 1. Du haut en bas, cing ballons de caoutchouc peu gonflés, un ballon de Polythéne, 
queue de gaze, nacelle et barométre enregistreur, aspect au départ. Approximativement a 
l’échelle. 

‘une force nette de 700 a 800 gr. en général, |’altitude d’éclatement étant alors 
environ 37 km. Ces ballons étaient attachés par des rubans Padoue, (voir Fig. 1), 
grace a une patte a cing doigts en méme ruban, collée avec du ruban type Durex, 
au sommet d’un ballon de Polythéne de 10 fuseaux, préparé au laboratoire, et 
gonflé ensuite, par les procédés déja décrits (3). Le ballon de Polythéne était 
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gonflé de facon a assurer une force nette au départ de 143 kg. Une grande queue 
en gaze surmontait la nacelle pour éviter les oscillations. Le départ a eu lieu de 
l’Université de Montréal. 







3. RESULTATS 





Sur six ascensions, trois ont été retrouvées, et permettent les conclusions 
suivantes, en partie décrites dans la liste que voici: 







No 3 





No 2 





No 1 











































Date 7 VII 52 17 VII 52 7 VIII 52 
Lieu d’atterrissage Glen Buell, Ont. _ Pierreville, Qué. Courcelles, Qué. 
Distance de Montréal 205 km. 91 km. 208 km. 
0.S.0. N.E. E. par E.N.E. 
Nombre de ballons de caoutchouc 2 ! 5 
Poids sous les ballons de caoutchouc 3900 gr. 4150 4150 
F.., force nette des ballons de 
caoutchouc 1700 2950 3600 
f,, contribution du ballon de Poly- 
théne au moment du plafonnement 2200 1200 550 
Altitude et pression prévues 28.1 km. 31.8 36.4 
0.0163 atm. 0.0089 0.0041 
Altitude et pression atteintes 22.1 km. 26.6 31.1 
0.04 atm. 0.02 0.01 
Départ 10 h. 30 H.A.E. 10 h. 05 10 h. 16 
Arrivée 23-24 h. 23 h. 10 23 h. 40 
environ 
Durée du vol 13 h. 12 h. 05 13 h. 24 
environ 
Couche; du soleil au sol et fin du 20 h. 47 20 h. 41 20 h. 15 
crépuscule 4 Montréal 23 h. 36 23 h. 14 22 h. 24 


La pression a été déterminée par un barométre 4 minimum du type imaginé 
par Piccard. Des tubes de verre fermés dans le haut laissant échapper I’air qu’ils 
contiennent au cours de la montée par le bas qui plonge dans de I’huile pour 
pompe a vide. A la descente I’huile remplit les tubes sauf une petite bulle ren- 
fermant I’air résiduel, dont la longueur est proportionnelle a la pression minimum. 
D’aprés notre expérience de cet appareil trés simple, dans d’autres ascensions, la 
pression minimum trouvée peut étre trop faible d’environ 0.005 atm. Nous avons 
négligé cette erreur dans |’évaluation de la pression atteinte, l’altitude réellement 
atteinte est probablement légérement inférieure a celle indiquée. Les observations 
faites sur les émulsions portées dans ces ascensions, en particulier sur les noyaux 
lourds, confirment les altitudes indiquées. 

Le systéme atteint son plafond en 1h.} a 2h.3, d’aprés l’observation au théodo- 
lite du No 1. Son équilibre est ensuite trés stable et évidemment le coucher du 
soleil seul le fait descendre en le refroidissant. Le coucher du soleil a lieu a ces 
altitudes prés d’une heure plus tard qu’au sol, soit environ 4 9 h.30 au lieu de 
8 h.30, et la descente prend environ deux heures. On dispose donc ainsi d’un 
moyen trés sir d’obtenir un plafonnement d’environ neuf heures ou davantage 
si le lancement est fait plus t6t, a une altitude réglable, sans avoir 4 manipuler 
un trés grand ballon de Polythéne. | 

L’altitude atteinte est notablement inférieure 4 la prévision théorique. Les 


pertes de gaz par des trous ou Aa travers les surfaces ne peuvent expliquer ce fait, 
puisque si elles suffisaient 4 abaisser le plafond, elles abrégeraient aussi le vol. 
Aprés considération des facteurs mis en jeu, l’explication parait étre la suivante: 
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le ballon de Polythéne acquiert sensiblement son volume maximum quand il 
prend forme, mais il ne garde pas sa forme pleine aux trés basses pressions ov il 
est entrainé, par suite de la réduction des tensions transversales de l’enveloppe. 
La tension verticale qui lui est imposée par les ballons accrochés 4 son sommet, 
et par la charge a sa base, contribue a lui donner une forme étirée. Un reméde 
partiel serait de soutenir le ballon de Polythéne par le haut, et la charge par une 
corde de longueur appropriée le traversant. 
4. VENTS DANS LA STRATOSPHERE 

Une remarque concerne le lieu d’atterrissage. Le vent entre 0 et 10 km. d’alti- 
tude avait une composante vers l’est, comme les vents de surface dominants, et 
une allure de l’ordre de 40 km. par heure. La proximité du lieu d’atterrissage et sa 
position suggérent fortement que les vents au cours du plafonnement présentaient 
une inversion, i.e., étaient dirigés vers l’ouest, d’accord avec une étude bibliogra- 
phique aimablement préparée par M. D. P. McIntyre. Selon lui, les vents strato- 
sphériques au-dessus de Montréal vont surtout vers l’est en hiver et surtout vers 
l’ouest en été. De plus l’ascension No 1 a permis d’observer cette inversion, elle 
est partie vers le N.N.E. puis est revenue exactement au zénith aprés 101 min., 
puis est disparue a l’ouest aprés 150 min. 

Ces expériences avaient pour but d’irradier des émulsions photographiques 
orientées au cours du vol, leurs résultats seront décrits en détail plus tard. Voir 
aussi (1, 2, 4, 5,6). 
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SUMMARY 


A simple procedure is described, combining the use of closed rubber balloons 
with that of an open Polythene balloon of medium size, which has allowed us to 
keep cosmic ray sensitive emulsions for 9 to 10 hr. at a ceiling of 24 to 32 km. 
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INTERIOR! 


By J. A. JACOBs 


ABSTRACT 


Bullen and Ramsey have shown that the reciprocal of the compressibility is a 
linear function of the pressure both in the Earth’s core and in the mantle below 
1000 km. In view of this result, it seems reasonable to suppose that a similar 
relationship exists between the reciprocal of the volume coefficient of thermal 
expansion and pressure. Support for this hypothesis is obtained by two indepen- 
dent methods (a) using Uffen’s results based on the theory of solids, (6) using 
Murnaghan’s theory of finite strain. Assuming the validity of the hypothesis, 
an estimate is made of the adiabatic gradient throughout the Earth. Taking the 
temperature at 1000 km. to be 3600° K., values of the temperature at greater 
depths are estimated. In particular it is found that the temperature at the 
boundary of the core is 4350°K., and at the center of the Earth a little over 
4800°K. These results have considerable interest and bearing on Bullard’s theory 
of the transfer of heat from the core. 


1. INTRODUCTION 

From a detailed analysis of seismic data, Bullen (4, 5) has derived both the 
density p and the adiabatic compressibility 8 of the Earth as functions of depth. 
In later papers (6, 7) he found empirically that the adiabatic incompressibility 
k, which is the reciprocal of 8, is a linear function of the pressure / viz. 


1.1] k = ky + ap. 


This relation was found to be valid both in the core and in the surrounding 
mantle below a depth of 1000 km. Ramsey (10, 11, 12) has re-examined this 
relation from the point of view of the theory of solids, and his conclusions point 
to the validity of equation [1.1]. He found, however, that the constant ky had 
different values in the core and mantle, although the derivative dk/dp = a was 
essentially the same. In amore recent and comprehensive paper on the constitution 
of the Earth’s interior, Birch (1) disagrees with the work of Bullen and Ramsey, 
and regards the relationship [1.1] as at best approximately valid over limited 
ranges of pressure for the particular materials of the Earth's interior. It is the 
purpose of this paper to try to obtain some information on the behavior of 
the thermodynamic coefficients with increasing depth and hence to make an 
estimate of the temperature distribution in the Earth. In view of equation 
[1.1], it seems reasonable to suppose that a similar relationship exists between 
the volume coefficient of thermal expansion a and pressure p, viz. that the 
reciprocal of @ is a linear function of p 
| 

[1.2] ae const. + bp = C + bp. 

Birch himself (1) has pointed out the importance of the thermal expansion in a 
number of geophysical problems, and little attention has been given to the 
question of the change of @ with pressure. Experimental data are scanty and 
subject to relatively large uncertainties. 


1 Manuscript received October 14, 1952. 
Contribution from the Geophysics Laboratory, University of Toronto, Toronto, Ontario. 
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It is because of this lack of data that the possibility of a relationship [1.2] is 
investigated, in spite of Birch’s objections to [1.1]. It will be shown that excellent 
support may be obtained for the temperature—pressure hypothesis [1.2], and 
although it is not suggested that it is a fundamental physical law, it is neverthe- 
less a good approximation over large ranges of pressure and should give an 
indication of the variation of a in the Earth’s interior—at least below 1000 km. 


2. ANALYSIS OF THE TEMPERATURE-PRESSURE HYPOTHESIS 
In a recent investigation on the thermal properties of the Earth (based on 
the theory of solids), Uffen (13) has obtained values of a within the mantle. 
These are reproduced in Table I. In Fig. 1, 1/a@ is plotted against p, and it can 
TABLE I 
VALUES OF @ IN THE MANTLE (AFTER UFFEN) 





Depth, | a 1/a, 

km. | i> af, 

200 | 36X10 | 2.78 x 10! 
600 -| 28 | 3.57 

1000 ~=—s|_—s(20 | 5.00 

1400. | «(17 | 5.88 

1800 | 14 | 7.14 

2200 | 12 | §.34 

2600 | 10 '10.00 


2900 Long (11.11 


‘ 


‘2200 Km. 


1800 Km. 


"1000 Km. 





0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 
p (10? dynes/cm?) 


Fic. 1. The reciprocal of the volume coefficient of thermal expansion as a function of the 
pressure in the mantle. 
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be seen that these results lend support to the hypothesis [1.2]. An alternative 
method was sought to corroborate this result if possible—this was made possible 
using Murnaghan’s (9) theory of finite strain. Birch (1) has used the same theory 
to obtain many of his results. Values of a are calculated from the equation 








YPC» 
2.1 a= 
[2.1] 2 
where y is Griineisen’s ratio and c, the specific heat at constant volume. We 
assume that we can neglect any changes in c,—a fact which is substantiated by 
the work of Uffen and Birch (ibid. 13, 1). 








Then % = 7? %°, where the suffix 0 refers to zero pressure. From Mur- 
ao Yo Po 






naghan’s theory, we can obtain expressions for p, 7, and p in terms of the strain 
€ or compression f, where we write f = — ¢€, a more convenient variable for 
compression. We have (see, for example, Birch (1), equation (15)) 


pb = 3kof (1 + 2f)°” (1 — 2f), 







2.2) Nap), = 2 + 49f — 282 + 32f + 81f)]/3[1 + Tf — 242 + OD], 
kr/p = ko/po. (1 + 2f) [1 + 7f — 28f (2 + OP)I, 


p/po = (1+ 2)”. 







These equations which are for isothermal compression express the pressure p 
and the ratio k7/p in terms of one parameter &. & depends on third order terms 
in the expansion of the Helmholtz free energy, i.e. in this case the strain energy, 
as a Taylor series in the strains. kp and & are functions of the temperature alone. 
The latest experimental work of Bridgman (2) shows that there is likely to be 
little error in taking £ to have its most probable value, viz. zero. To obtain 
values of y, we use equation 18 (Birch (1)), viz. 

1, ldlogk 3 do 


fe 
[2.3] Y= —L4idloek_ 2,42] ( — 26)(1 + «) | 

























which is obtained from Debye’s theory of the thermal expansion of solids. 

Little is known of the variation of Poisson’s ratio ¢ with p, but by plotting o 

against depth (found from seismic data), the graph is approximately linear 

between about 1000 and 2600 km. In this range o increases by about 0.024 and 

p by 0.86. We thus take do/dp to be approximately constant at a value of 

0.0279. The initial value of y is taken as 1.5. Finally d log k/d log p = (dk7/dp)r, 

so that we can construct Table II. | 
Values of f are calculated from the depth on the assumption that the compres- 

sion is adiabatic. 
The values of ao/a obtained in Table II merely reflect the changes in a due 

to compression and do not take into account the effect of temperature changes, 

since we have used the isothermal incompressibility k7 instead of the adiabatic 

incompressibility &. The effects of temperature will, however, be small com- | 

pared to those of pressure. 





SN Se ST SD 


| 
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TABLE II 
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CHANGE OF GRUNEISEN’S RATIO AND THERMAL EXPANSION WITH COMPRESSION (USING 


MURNAGHAN’S THEORY OF FINITE STRAIN) 

















| 
Depth, f p/Ro y/yo 

km. 

| 
1000 0.0199 | 0.186 0.851 
1200 0.0590 | 0.234 0.823 
1400 0.0688 | 0.285 0.797 
1600 0.0775 | 0.333 0.773 
1800 0.0853 | 0.380 0.752 
2000 0.0930 | 0.428 0.731 
2200 0.1007 | 0.478 | 0.714 
2400 0.1088 | 0.534 0.691 
2600 0.1167 | 0.592 0.676 
2800 0.1271 | 0.672 0.659 





Fic. 2. 











pko/ pok ao/a 
.678 1.735 
.632 1.920 
593 2.110 
561 2.305 
534 2.490 
512 2.675 
.488 2.865 
467 3.100 
447 3.310 
.422 3.595 





The reciprocal of the volume coefficient of thermal expansion as a function of the 


pressure in the mantle. 


Fig. 2 shows ao/a plotted against p/k» and again excellent agreement is found 
with the linear hypothesis expressed in equation [1.2]. Using Uffen’s values for 


a, we have from Fig. 1, 


[2.4] 


Using Murnaghan’s theory of finite strain, we have from Fig. 2 


jc 
Lo 





= 2.4 X 10‘ °K.., 





= 6.2 X 10-° °K. cm.? per dyne. 


% 143922, 
Qa 0 
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3.92 , 

aoko “" 
Values of ao and ko are not easy to estimate precisely. Using Birch’s (1) figures, 
we can take ap = 36 X 10-® (°K.)—!, ko = 2 X 10'? dynes per cm.? With these 
values, equation [2.5] can be written 


= C + bp where 
Qa 


(C = 2.78 x 10'°K,, 


2.6 Sis =f) Ore 
2,6] [2 = 5.45 X 107° °K. cm® per dyne. 
The agreement between [2.4] and [2.6] is considered good enough in view of 
the approximations in the theory and in the uncertainties in some of the elastic 
and thermal coefficients. 
3. THE ADIABATIC TEMPERATURE GRADIENT 

It is now possible to make an estimate of the adiabatic temperature gradient 

throughout the Earth. The adiabatic gradient may be written 


[3.1] 


where c, is the specific heat at constant pressure. From [1.2] and [8.1], we may 
write 
[3 2] dT =- ee 

= T pc (C + bp) 
Both Uffen (13) and Verhoogen (15) have obtained values of the ratio a/c, in 
the mantle using seismic data. Uffen found that the values of c, were approxi- 
mately constant (in agreement with Bridgman’s observations at relatively low 
pressures), and that the variations in a/c, were almost entirely due to variations 
in a. Assuming that we can take c, as constant, we have from [3.2] 


3 .T=t{[—e_— 

[3.3] log, T see aay 

It is estimated that variations in c, will not affect the validity of [3.3] by more 
than 7% if the temperature rise is adiabatic. Since p is known as a function of 
depth (and hence of pressure), equation [3.3] can be integrated numerically. 
The values of C and } were taken from equation [2.4]. In order to obtain a 
numerical value for T at any depth, its value at some depth must be known. 
Taking T at 1000 km. to be 3600°K. (Uffen (13)), values of T can be estimated 
at greater depths. In particular the temperature at the boundary of the core ~ 
and mantle is found to be 4350°K. 

If we assume in addition that the relationship [1.2] is also valid in the core 
as Bullen postulated [1.1] to be, then we can integrate equation [3.2] throughout 
the core as well. It is found that the temperature at the center of the Earth is a 
little over 4800°K. The increase throughout the core is thus only 500°. The 
results of the integration are presented graphically in Fig. 3. We thus find that 
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5,500 
5,000 


4,500 637) Km. 
4800 Km. 
Temperature T| 
oe 


4,000 2900 Km. 
! 


=< Mantle: Core — 


3,500 


1000 Km. 


3,000 
0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 


Pressure p ( 10"? dynes/cm?) 


Fic. 3. Temperature distribution below 1000 km. 


the temperature at the center of the Earth is of the order of 5000°K. rather 
than the accepted value of about 10,000°K. Valle (14) has also obtained values 
of the ratio a/c, in the mantle using seismic data, and has calculated the adia- 
batic variation of the ratio of the temperature at any depth in the mantle to the 
temperature at 33 km. He has extended his theory to the core to obtain values 
of the ratio of the temperature at any depth in the core to its value at the 
boundary of the core and mantle. He obtains over 30% increase in the temper- 
ature at points halfway across the core which seems very excessive. 


4. CONCLUSIONS 


Bullard (3) and Elsasser (8) have considered the question of the transfer of 
heat from the core and examined the possibility of convection currents. Bullard 
found that the heat flow required to maintain convection currents was about 
four times greater than that available from a core with a radioactive content 
comparable with that of iron meteorites. He used a temperature of 10,000° in 
his calculation of the adiabatic gradient, and a value of 5000° would thus reduce 
his discrepancy by a half. Moreover, Bullard used the value of 1.6 for y in his 
estimation of a. Both Uffen’s work and the above analysis based on the theory 
of finite strain show that y decreases with depth so that Bullard’s value of 1.6 
is too high. The value of y is more likely to be about 0.8, and this would remove 
any discrepancy in Bullard’s estimate of the heat flow. 

In conclusion it must be emphasized that no physical proof of the validity of 
the relation [1.2] throughout the Earth’s interior has been given. Two indepen- 
dent analyses however agree in suggesting that it may be a very fair approxima- 
tion to the truth—at any rate in the mantle below 1000 km. In view of the 
scarcity of data on this question, the results of this paper should be of value in 
giving at least some indication of the behavior of a and y with pressure. 
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A DETERMINATION OF THE HALF LIVES OF SOME MAGNETIC 
DIPOLE y-RAY TRANSITIONS! 


By R. L. GRAHAM AND R. E. BELL 


ABSTRACT 


A number of magnetic dipole y-ray transitions have been studied using a 
coincidence circuit of short resolving time (27% = 2 X 10~* sec.), a two lens 
single 8-ray spectrometer, a pair of lens spectrometers placed end to end with 
coincidence counting of the focused radiations from a single source, anda scintilla- 
tion spectrometer. Lifetimes have been measured using the delayed coincidence 
method and where feasible conversion coefficients and K/L ratios obtained. 
Comparison is made with theoretical estimates of the lifetime-energy relation 
for M1 y-ray transitions. 

The new results are as follows: : 











Nucleus -Ray | Half life, sec. a | K/L 
energy, kev. | K 

Tels | 159 } (1.9 +0.3) X 10-10 — oe 
Tels | 35.4 (1.58 + 0.15) X 1079 | — poem 
Cs183 |} 81 } (6.0 +0.4) xX 10-9 | 1.77 7.5 
Css 248 | (2.8 +0.8) xX 10-” — 7.0 
Pm? | 91.5 2 Ps + 0.08) X 10-° 1.8 | 7.3 
T1208 40 (A line) | 7 X10 — | _ 
(Thc” | 

Bit? | 238 (F line) <2 xX 107 — | _ 
(ThC) | | 


I. INTRODUCTION 


Upper limits have been reported for the lifetimes of a number of magnetic 
dipole, M1, y-ray transitions (31) but until now actual lifetimes have been 
determined for only three M1 y transitions. These are the 478 kev. transition 
in Li? with Ty. = (5.2 + 1.7) X 10~" sec. (12), the 14 kev. transition in Fe*’ 
with 71,2 = 1.1 X 10-7 sec. (9), and the 80 kev. transition in Xe!*! with 
Ti2 = (4.8 + 2.0) X 10" sec. (17) 

This paper reports the results of experiments on seven other J/1 transitions, 
five of which were found to have measurable half lives. Where the multipolarity 
was uncertain assignment has been made by one or more of the following 
methods: (a) comparison of the experimental K conversion coefficients with 
the theoretical values deduced from the calculated values of Rose et al. (30) and 
Spinrad and Keller (35), (6) comparison of the experimental K/Z ratio of 
conversion coefficients with empirical curves of Goldhaber and Sunyar (16), 
(c) comparison of the experimental L;: L,;: L;;;, conversion coefficient ratios 
with the theoretical values calculated by Gellman ef al. (15) which have been 
verified in part experimentally by Mihelich and Church for y-rays of known 
multipolarity (26, 27). The delayed coincidence method of measuring short 
half lives is reviewed. It is seen that scintillation counters employing fast 
organic phosphors and recently developed delayed coincidence techniques make 
it possible, under favorable circumstances, to measure lifetimes as short as 
10-'° sec. with some accuracy and to set upper limits of ~ 10-"™ sec. thus 
bringing many M1 transitions within range of experimental detection. 


1 Manuscript received October 10, 1952. 
Contribution from Physics Division, Atomic Energy of Canada Limited, Chalk River, 
Ontario. Issued as A.E.C.L. No. 11 
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The experimental half lives of 1/1 transitions in odd A nuclei are ~ 100 times 
longer than those predicted theoretically for ‘‘allowed” transitions and appear 
to be slightly longer for odd proton nuclei than odd neutron nuclei. The discre- 
pancy between theoretical and experimental lifetimes of ‘forbidden’ M1 
transitions (A/ = 2 in shell model orbital assignment) in heavy odd neutron 
nuclei has been used by Sachs and Ross (31, 33) to provide evidence for the 
nonadditivity of nucleon moments and the existence of exchange moments. 


II. METHOD 


The measurement of short half-lives by the delayed coincidence method has 
been discussed fully in earlier publications (3, 17). When the lifetime is long 
compared to the resolving time of the coincidence circuit, the half-life 74/2 
may be deduced most accurately from the slope of the tail of the coincidence 
resolution curve F(x) plotted on a logarithmic scale, where x is the artificially 
inserted delay time. When the lifetime is short compared to the resolving time, 
the shift of the centroid of the delayed resolution curve F(x), with respect to 
that of a “prompt”’ curve, P(x), is a direct measure of the mean life, 7, and may 
be determined from the experimental data by the methods of Bay (1) and 
Newton (28). Radiations are considered to be “prompt”’ if their lifetimes are 
too short to be measured by the delayed coincidence technique. 

For the shortest lifetimes it is often convenient to use the more sensitive 
“self comparison’? method in which the delayed resolution curve F(x) is com- 
pared with its own inverse F(— x). This method was first used (17) to place an 
upper limit on the lifetime of the 411 kev. excited state in Hg!®* and is briefly 
as follows: A radioactive source of strength 10 to 100 uc. is placed at the center 
of the pair of B-ray spectrometers mentioned below. The delayed curve F(x) is 
then obtained with the conversion line of the y-ray transition focused in one, 
the north, spectrometer and a point on the 6-ray spectrum just below the line 
focused in the other, the south, spectrometer. The inverse curve F(— x) is 
obtained by increasing the currents in both spectrometers slightly so that the 
north spectrometer focuses on the 8 spectrum just above the line while the 
south focuses on the line. Among the advantages of this method are (1) it 
doubles the effect (the shift of the centroids of the two resolution curves is then 
27 = 2.89 Ti,2), (2) it cancels out systematic errors, and (3) it eliminates the 
need for obtaining the prompt resolution curve. 


III. APPARATUS 

In most of the investigations described below the radiations were selected 
magnetically in a pair of identical lens B-ray spectrometers placed end-to-end, 
which have a common vacuum chamber (2, 18). This apparatus is quite flexible; 
with the source at the center and the two coils excited separately, radiations of 
different energies can be focused simultaneously on the stilbene-1 P21 scintillation 
counters at either end of the vacuum chamber. With the source at one end, the 
coils connected in series, and a change in the baffle system it behaves as a single, 
two-lens 6-ray spectrometer, focusing radiations on one counter only. As 
described elsewhere (2) the current in one or both coils can be made to assume 
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up to 35 different fixed values successively by means of a stepping switch 
control which may be operated automatically. 

The scintillation counters were connected to a coincidence circuit (3, 4) 
whose resolving time, 279, can be varied down to 2 X 10~° sec. or less without 
reducing the coincidence efficiency below 0.9. Artificial time delays were inserted 
between both counters and the coincidence mixer by means of the motor driven 
helical delay unit described in an earlier publication (3). 

Almost all of the data in these experiments were recorded on Esterline Angus 
charts (2). A timing switch was connected to the delay unit to change the 
inserted delay every 10 min. After counting had proceeded for 10 min. at each 
of the chosen delay settings the delay unit reset itself automatically and repeated 
the sequence in the same order. A switch actuated by the delay unit drive 
mechanism and connected to the current stepping control allowed the current 
setting to be changed while resetting the delay unit. This was used in the self 
comparison determinations to rock the current back and forth after each sweep 
of a resolution curve. 

A scintillation y-ray spectrometer consisting of a 1 in. dia. by 1 in. Nal (TI) 
crystal, a 5819 photomultiplier, and a single channel pulse height analyzer was 
used to measure the K conversion coefficient in two cases. 






















IV. ERRORS 


The accuracy with which short half lives can be measured is limited by the 
statistics and also by certain systematic effects not all of which are completely 
amenable to calculation. The propagation velocity in the air spaced coaxial 
transmission line in the helical delay unit is very nearly that of light. A slight 
reduction of velocity due to the presence of the Lucite pegs supporting the 
central conductor is calculated to be 1%. Inaccuracies in calibration due to the 
mechanical dimensions of the delay unit are estimated to be less than 1% and 
due to the uncertainty in positioning the movable carriage less than 10"? sec. 
The average time interval, ¢, between the excitation of the phosphor and the 
emission of the first photoelectron from the photomultiplier cathode is depen- 
dent upon the energy E (kev.) of the incident electron expended in the phosphor 
(29). With trans stilbene as a phosphor and fairly energetic electrons we expect 
(3) for our counters f ~ (3/E) X 10-8 sec. f is measured experimentally as the 
reciprocal of the slope of the sides of a prompt resolution curve plotted on a 
logarithmic scale and is found experimentally to be of the order estimated above 
indicating that the circuit responds to the first photoelectron. Corrections due 
to variations in @ are avoided wherever possible by comparing delayed and 
prompt radiations of the same energy. In the case of the self comparison method 
the changes in { in the two counters roughly cancel since the energies of the 
radiations focused are increased by approximately the same small amount. 

‘For the highest accuracy, however, it is necessary to make second order correc- 
tions if the ?’s for the two counters are different and errors of ~ 10-" sec. or 

greater may arise. Transit time of the electrons in the spectrometer must also 

be considered when the current is changed in making a comparison. Since the 
momentum (and hence velocity) of the electrons being focused and their path 
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length are known this correction can be made with an accuracy usually better 
than 10-" sec. In order to minimize errors due to other effects such as slow drift 
in photomultiplier characteristics and thermal changes, comparisons were made 
as close in time as possible. In self comparison experiments the runs were 
sandwiched so that these corrections would tend to cancel out. With these 
precautions using the self comparison method it would seem possible to set 
limits as small as a few times 10—" sec. on lifetimes and to measure lifetimes with 
standard deviations of ~ 10~-" sec. if the statistical accuracy is sufficiently 
good. The errors assigned to lifetimes determined by a centroid shift analysis 
are the combined standard deviations due to (a) statistics, (0) corrections in 7? 
whose error is estimated from the standard deviations of the points on the side 
of the resolution curve from which é is measured, (c) slow drifts in time as 
measured by the shift of the centroid of a particular resolution curve from day 
to day, and (d) uncertainty in the spectrometer transit time correction which 
is appreciable only for low electron energies. The upper limits on lifetimes 
quoted below are taken to be twice the total combined errors described above. 
V. RESULTS 

(a) Te!*8, 159 Kev. y-Ray Transition 

The 159 kev. y-ray transition in Te!** has been shown to follow the 90 day, 
88.5 kev. y-ray transition in Te!** by Katz, Hill, and Goldhaber (20) and has 
been classified as M1 radiation (shell model orbital d3;2 — 51/2) on the basis of 
its conversion coefficient and K/L ratio (19, 20). Deutsch and Wright (10) 
found the half life to be < 4 X 10~-® sec. and McGowan found it to be < 10-° 
sec. (24). 

A source of neutron irradiated Te of strength ~ 20uc. was placed at the 
center of the pair of B-ray spectrometers. With the K 88.5 line focused in the 
north spectrometer and the K 159 line focused in the south spectrometer the 
coincidence resolution curve F(x) in Fig. 1 was obtained. Two prompt compari- 
sons were made using sources of thorium active deposit. The thorium F line 
transition is shown below to have a half life 71,2. << 2 X 10-" sec. In one com- 
parison, no changes were made in spectrometer currents so that the north 
spectrometer focused 8-rays from ThB and the south spectrometer focused the 
F line which was degraded by an Al absorber placed over the source so that it 
corresponded in energy to the K 159 line. The results are shown as the prompt 
resolution curve P(x) in Fig. 1. In the other case the F line was not degraded 
and the south spectrometer current was raised to bring the F line into focus. 
The centroid of the resolution curve obtained in the second case (not shown) 
was found to agree within the limits of statistical accuracy with that of the 
first comparison after correcting for the changes in transit time in the south 
spectrometer and in / in the south counter. The shift between the centroid of 
F(x) and a weighted average of the two centroid determinations for P(x) gives 
for the mean life of the 159 kev. y-ray transition in Te’, 7 = (2.7 + 0.3) X 
10-!° sec. The total error in the mean life including the uncertainty in the mean 
life of the F line transition (~ 3 X 10-" sec.) and drift of the centroid of the 
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Fic. 1. Delayed resolution curve F(x) for the 159 kev. y-ray transition of Te! with a 
prompt curve P(x) for comparison. The half life of this transition is measured by the shift 
of the centroid of F(x) to the right. The result is Ti. = (1.9 + 0.3) & 107° sec. The standard 
deviations of the points are shown by vertical bars. 


resolution curves with time (~~ 2 X 10-' sec. per day) is ~ 4.7 X 107" sec. 
Hence we conclude that the half life of the 159 kev. transition in Te!?* is 


T12 = (1.9 + 0.3) K 107! sec. 


(6) Te!®> 35.4 Kev. y-Ray Transition 

The 35.4 kev. y-ray transition has been shown to follow the 58 day 109.7 
kev. y-ray transition in Te!*> by Bowe and Scharff-Goldhaber (8) and has been 
classified as predominantly 7/1 radiation (shell model orbital d3;2 — s1;2) on the 
basis of its conversion coefficient and K/L ratio by Bowe and Axel (7) and also 
by Friedlander et al. (14). The half life for this transition has been shown to be 
<5 X 10-* sec. by Deutsch and Wright (10) and < 2 X 10~° sec. by McGowan 
(24). 

The delayed resolution curve F(x) shown in Fig. 2 was obtained using the 
source mentioned in (a) above with the north spectrometer focused on the 
L 35.4 line and the south spectrometer focused on the Z 109.7 line. The prompt 
resolution curve P(x) shown in Fig. 2 was obtained when the Te source was 
replaced by a source of ThB and the current in the south spectrometer was 
increased in order to focus the F line. To improve the signal to noise ratio in 
the north counter in this experiment the photomultiplier was operated at a 
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Fic. 2. Delayed resolution curve F(x) for the 35.4 kev. y-ray transition of Te™ with a 
prompt curve P(x) for comparison. The half life of this transition as measured from the shift 
between the centroids of F(x) and P(x) is Tis2 = (1.58 + 0.15) X 10~® sec. The standard 
deviations of the points are shown by vertical bars. 


lower voltage (1200 v.) than is customary (3) and the extra gain required for 
proper performance of the coincidence circuit was obtained by two wide band 
(200 Me. per sec.) amplifiers, each of gain 10, in cascade. The half life for the 
35.4 kev. transition as deduced from the slope of the tail of F(x) is ~ 1.6 X 10° 
sec. A centroid analysis of F(x) and P(x) yields, after correcting for the change 
in transit time and 7 in the south spectrometer, 71,2 = (1.58 + .11) X 10~-® sec. 
To allow for the uncertainty in the F line lifetime and a possible error in the ¢ 
correction of ~ 10~!° sec. we increase the limits of error and obtain for this 
transition 
T12 = (1.58 + 0.15) X 10-® sec. 


Because of the thick (~ 2 mgm. per cm.*) source used in this experiment the 
L 35.4 line was not too well resolved in the 8-ray spectrum. In order to check 
that the coincidences were in fact due to this transition the north spectrometer 
was swept over the region of the Z 35.4 line and the line was clearly outlined in 
coincidence with the Z 109.7 line focused in the south spectrometer. 
(c) Cs'83, 81 Kev. y-Ray Transition 

5.3 day Xe!*’ decays mainly by 8 emission, 745 kev., to an 81 kev. 
excited state of Cs!** (6). The 81 kev. y-ray transition to the ground state of 
Cs153 has been classified by Bergstr6m and Thulin (6) as a mixture of M1 and 
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E2 radiation (d5;2 — gz;2) on the basis of the conversion coefficient, 2.9, and the 
K/L ratio of 5.9. 

As a further check on the multipolarity assignment of this transition the 
conversion coefficient and K/L ratio were determined using a source of fission 
product iodine of initial strength ~ 100uc. which was allowed to stand for some 
days until 5.3 day Xe!** was the predominant activity. By comparing the areas 
of the ~ 32 kev. K X-ray peak and the 81 kev. y-ray peak in the scintillation 
spectrum and making corrections for absorption, escape peaks, and fluorescent 
yield in the same manner as discussed elsewhere for the 84 kev. y-ray transition 
in Yb!7° (18) we obtain for the K conversion coefficient 

ax = 1.77 + 0.15. 


The 8-ray spectrum of the source, examined at this time, showed that there 
was no appreciable contribution to the K X-ray peak in the scintillation spec- 
trum due to internal conversion of other y-ray transitions in the source. The 
experimental value of ax is in good agreement with the theoretical figure for 
M1 radiation of 1.75 obtained by interpolating between the conversion coeffi- 
cients calculated by Rose et al. (30) and the threshold value of Spinrad and 
Keller (35). The K/(Z + M+ N) ratio as determined by the areas of the 
two peaks in the 8-ray spectrum (counts/momentum vs. momentum) is 6.0. 
Assuming that 0.8 of the ZL + M+ WN peak is due to Z conversion we obtain 


10° 
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Fic. 3. Delayed resolution curve F(x) for the 81 kev. y-ray transition of Cs with a 
prompt curve P(x) for comparison. From the slope of the tail of F(x) we get for the half life 
Tij2 = (6.0 + 0.4) XK 10-® sec. Where the standard deviations of the points are larger than 
the radii of the circles they are shown by vertical bars. 
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K/L = 7.5 in good agreement with the value of 7.5 predicted by Goldhaber and 
Sunyar’s empirical curve (16) and conclude that the 81 kev. y-ray transition 
in Cs!88 is predominantly M1 in character. 

In order to measure the lifetime of this transition the source was placed 
directly over the crystal of the north counter which was covered with a 7 mgm. 
per cm.? Al absorber to prevent conversion electrons from the 81 kev. y-ray 
transition from being counted. The delayed coincidence resolution curve F(x) 
Fig. 3 was obtained when the K81 line was focused on the south counter by the 
two-lens single spectrometer and only 6 rays were being counted in the north 
counter. The prompt curve P(x) shows, on a magnified ( X 60) scale, the results 
obtained when @ rays just above the K 81 line were focused on the south counter. 
The prompt coincidences observed here are due in part to other prompt radia- 
tions of low intensity in the source and in part to electrons back-scattered from 
the north crystal. From the slope of the tail of F(x) we obtain for the 81 kev. 
y-ray transition in Cs!*8 


Ti,2 = (6.0 + 0.4) X 107° sec. 
This figure was confirmed by another delayed coincidence experiment in which 


the source was placed directly between the crystals of the two counters which 
had been removed from the spectrometer. 


(d) Cs'*5, 248 Kev. y-Ray Transition 
9.2 hour Xe'® decays by 8 emission, 905 kev., followed by a 248 kev. 
y-ray transition in Cs'** (5). No estimates of the multipolarity of this transition 


have been reported. 

The M1 character of this transition was deduced by measuring the areas of 
the K and L + M + N peaks in the same fission product iodine source mentioned 
in (c) above after waiting a few hours to allow Xe!** to grow in from the decay of 
6.7 hour I'**. The ratio of peak areas gives K/(L + M+ N) = 5.6 + 1.0 and 
hence K/L = 7.0 + 1.2 if we assume that 0.8 of the L + M+ N peak is due 
to L conversion. This is to be compared with the values of 7.8 and 4.0 expected 
for M1 and £2 radiations respectively. No measurement of the conversion 
coefficient by the scintillation spectrometer method was attempted because of 
the presence of other y activities in the source at this time. The experimental 
K/L ratio indicates that this transition is predominantly M1 in character but 
is not accurate enough to exclude some admixture of £2 radiation. 

The lifetime of this transition was determined by the self comparison method 
using the K conversion line. The source was placed 1.5 cm. to the north of the 
center of the pair of spectrometers whose coils were connected in series. Data 
were collected over a period of about 60 hr. The delayed resolution curve F(x) 
in Fig. 4 shows the results obtained when the K 248 conversion line was focused 
in the south spectrometer and 6-rays just below the line were focused in the 
north spectrometer. The curve F( — x) was obtained with a slightly higher 
spectrometer current so that the south spectrometer focused 8-rays just above 
the line and the north spectrometer focused the K 248 line. From the shift 
between the centroids of F(x) and F( — x) we obtain for the mean life 7 = 
(4.0 + 1.0) X 10-!° sec. The estimated error due to differences in transit time 
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Fic. 4. Self comparison resolution curves for the 248 kev. y-ray transition of Cs!%, 
The curve F( — x) is the inverse of the delayed resolution curve F(x) and its centroid should 
lie to the left of the centroid of F( — x) by 2.89 T1,2. The result is T1/2 = (2.8 + 0.8) X 1071 
sec. The standard deviations of the points are shown by vertical bars. 


and 7? in this experiment is at most ~ 10-"™ sec. so that only the statistical 
error is important. Hence we conclude that the half life of the 248 kev. y-ray 
transition in Cs! is 


Ti2 = (2.8 + 0.8) X 10-!° sec. 


(e) Pm'*", 91 Kev. y-Ray Transition 

11 day Nd'4? decays in part ( ~ 60%) by 8 emission, 825 kev., to a 91 
kev. excited state in Pm!47 (11). The 91 kev. y-ray transition to the ground state 
has been classified as M1 or M1 + E2 radiation by Kondaiah (21), who found 
K/L = 6.5 + 1.5 and ax ~ 0.9, and also by Mihelich and Church (27) who 
found K/L = 6.4 + 2.0 and that the Z conversion was predominantly J). 
These findings are to be compared with the predicted K/L ratios of 7.7 for M1 
and 0.5 for £2 radiation (16), the theoretical K conversion coefficient of 1.9 for 
M1 radiation (30, 35), and the theoretical Z shell conversion coefficients of 
Gellman et al. (15), which show that M1 radiation is converted predominantly 
in the L, shell while E2 radiation is converted mainly in the L,;,; and L; 7; 
shells. Emmerich and Kurbatov (11) found that the lifetime of this transition 
is <5 X 10-7 sec. and McGowan (23) found no 6-y decays in the region 
3 X 10-' sec. to 3 X 107? sec. 

The 8-ray spectrum was examined in the 8-ray spectrometer using a source 
of Nd!47 prepared from neutron irradiated neodymium oxide about two weeks 
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after completion of the irradiation. In this way interference from shorter lived 
activities was avoided. From the areas of the K, Z, and M conversion lines, 
which were clearly resolved, we obtain K/L = 7.3 and L/M = 4.0. From the 
areas of the ~ 40 kev. X-ray peak and the 91 kev. y-ray peak in the scintillation 
spectrum we obtain after the usual corrections ag = 1.8 + 0.15. These values 
are in good agreement with the theoretical values for M1 radiation mentioned 
above and verify that this transition is predominantly M1. 

A preliminary measurement of the half life of the 91 kev. y-ray transition 
with a ~ 10uc. source at the center of the pair of 8-ray spectrometers gave from 
the slope of the tail of the delayed resolution curve (not shown) T1,;2 = (2.5 + 
0.4) X 10~-° sec. Having identified this lifetime with the 91 kev. transition the 
counters and source were removed from the spectrometer and the source was 
sandwiched between the stilbene crystals of the two counters. A 7 mgm. per 
cm.? Al absorber was placed over the north counter to prevent conversion 
electrons of the 91 kev. y transition from being counted. The delayed resolution 
curve of coincidences between conversion electrons counted in the south counter 
and 8 rays counted in the north counter is shown as F(x) in Fig. 5. The curve 
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Fic. 5. Delayed resolution curve F(x) for the 91 kev. y-ray transition of Pm'” with a 

prompt curve P(x) for comparison. From the slope of the tail of F(x) we obtain for the half 
life of this transition T1/2 = (2.44 + 0.08) X 10~* sec. Where the standard deviations are 
larger than the radii of the circles they are shown by vertical bars. 
P(x), obtained when the Nd!*" source was replaced by one of thorium active 
deposit, illustrates the performance of the system with prompt radiations. 
From the slope of the tail of F(x) we obtain for the 91 kev. y-ray transition in 
Pm!*7 


Ti. = (2.44 + 0.08) X 10-* sec. 


(f) TI? (ThC”) 40 Kev. y-Ray Transition 
ThC decays in part by a emission, 6.04 Mev., to a 40 kev. excited state in 
ThC”’. The 40 kev. y-ray transition to the ground state is highly converted in 
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the L shell. Ellis (82) found that the intensities of the Z;, L,,;, and L;,;; con- 
version lines are respectively ‘‘very strong”, ‘‘strong’’, and ‘‘medium”’ which is 
the sequence of intensities for M1 transitions and unlike that for £2 transitions, 
as shown by the conversion coefficients calculated by Gellman et a/. (15) and the 
experiments of Mihelich and Church (26, 27). We therefore tentatively identify 
thisasan M1 transition. The A line intensities measured by Flammersfeld (13) and 
the a-ray intensities of ThC measured by Lewis and Bowden (22) suggest that 
this transition is almost entirely internally converted. This is consistent with 
the high L; conversion coefficients (> 10) calculated by Gellman et al. (15) 
for M1 radiation or with those deduced for L; + L,;; + L,7, shells using the 
known M1 conversion coefficients for the K shell (30, 35), the K/L ratio found 
empirically by Goldhaber and Sunyar (16), and extrapolating down to 40 
kev. assuming an E~-* dependency. 

A source of thorium active deposit of strength ~ 10yuc. was placed over the 
stilbene crystal of the south counter and the A line conversion electrons were 
focused on the north counter by a two coil, single spectrometer. Despite the low 
energy of the A line electrons and poor signal to noise ratio due to the high 
voltage (2000 v.) applied to the photomultiplier it was possible to observe the 
line with a counting rate of 1000 to 4000 c.p.m. on a background of noise and 
continuous 6-ray spectrum of 2000 to 4000 c.p.m. The resolution curve F(x) in 
Fig. 6 shows coincidences between the a particles registered in the south counter 
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and the A line electrons focused in the north counter plotted on a logarithmic 
scale as a function of inserted delay. The prompt curve P(x) in Fig. 6 is a 
composite of coincidences observed between F line conversion electrons in the 
south counter and £ rays just below (solid circles) and just above (open squares) 
the A line after correcting for the change in spectrometer transit time and in 
i in the north counter. The coincidence counting rate off the line, P(x), was 
about 1/5 that on the line, F(x); the data of F(x) include this prompt coinci- 
dence contribution. Analysis of the data according to the method of Bay (1) 
gives for the mean life of the A line transition r = ( — 1.8 + 1.8) X 10~™ sec. 
which also includes the lifetime of the F line transition (see below). The largest 
uncertainty (~ 5 X 10-" sec.) is in the # correction applicable to the south 
counter because of the slightly different pulse heights in the two cases. We may 
conclude then that the mean life of the 40 kev. (A line) y-ray transition in 
T12°8 is less than 10~-!° sec. and hence 
Ti;2 <7 X 10-" sec. 

(g) Bi?!? (ThC) 238 Kev. y-Ray Transition 

ThB decays mainly by B~ emission, 0.33 Mev., followed by a 238 kev. y-ray 
transition, whose K conversion line is the well known F line, to the ground state 
of ThC (Bi?!*) (36). Martin and Richardson (25) have measured the internal 
conversion coefficient of this transition and find that it agrees reasonably well 
with the theoretical coefficient for M1 radiation calculated by Rose (30). 
Moreover estimates of the K/L ratio of this transition compiled by Way (36) 
cluster near the value of about 7.5 found by Goldhaber and Sunyar (16) for 
other M1 transitions. However the possibility of some admixture of £2 radiation 
remains until better agreement is obtained between measured and theoretical 
conversion coefficients. 

The upper limit in the half life of this transition was previously reported as 
Ti;2 < 107! sec. (3, 17). Four subsequent measurements have been made using 
the self comparison method and stronger sources ( ~ 200uc.) of thorium active 
deposit in order to improve the statistical accuracy. The two spectrometer coils 
were connected in series in all cases. The momentum of the F line electrons in 
one spectrometer was degraded by about one line width (5%) by an Al absorber 
of 3-5 mgm. per cm.? placed over the source. In one case the source position was 
also slightly asymmetrical. The data for each run were accumulated over a 
period of about 24 hr.,.and the shift between centroids of the two resolution 
curves was analyzed by the method of Bay (1) carrying the standard deviations 
of the points through the numerical computations. After making corrections 
for the small difference in ¢ in the two counters and for the difference in spectro- 
meter transit time in one case, we obtain a weighted average for the mean life 
of this transition of 7 = ( — 13+ 0.4) X 10-" sec. or 71,2 = ( — 0.9 + 0.3) 
X 10-!! sec. This result indicates the presence of systematic errors of ~ 1.3 X 
10-! sec. which are of the order suggested in the discussion of errors given in 
Section IV. We may conclude, therefore, that the mean life of the 238 kev. 
y-ray transition in ThC is less than 3 X 10-"! sec. and hence 
Ti)2 < 2X 10-" sec. 
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VI. DISCUSSION 


The information obtained above together with other available information 
about M1 y-ray transitions is compiled in Table I. Since the measure of the 
K/L ratio seems to be a sensitive test for the presence of an admixture of E2 
radiation, the experimental values were compared with Goldhaber and Sunyar’s 
empirical curves for M1 and £2 transitions (16). 7, is here defined as the mean 
life in seconds for photon emission and is deduced from the experimental half 
life after correcting for internal conversion using the adopted total conversion 
coefficient listed in the table. Where possible shell model orbital assignments for 
the initial and final states are given. 

Some additional information on M1 transitions has been included from the 
case of Hg!*® based on the upper limit for the half life of the 209 kev. excited 
state previously reported (3) and the disintegration scheme for the decay of 
Au!*® proposed by Sherk and Hill (84) which has been borne out by coincidence 
spectrometer studies in this laboratory (3). A 50 kev. and a 209 kev. y-ray 
transition both proceed from the same excited state; the upper limit in the mean 
life for each transition is deduced using the ratio of y-ray intensities found by 
Sherk and Hill. 

The estimate of the K conversion coefficient ( ~~ 14) used by Goldhaber and 
Sunyar (16) for the 14 kev. y-ray transition in Fe*’ is probably high. The value, 
a,x = 8.5, used here was deduced from the threshold value for Z = 26 (35) 
assuming an E~-* dependency. 

In Fig. 7 log 7, is plotted as a function of y-ray energy on a logarithmic scale. 
For comparison, Weisskopf’s relation (37) for magnetic multipole transitions, 


5 





¢ ODD Z EVENN 
© EVEN Z ODDN 

“6 6 OD0Z ODDN 
+ UPPER BOUND 


“7 


“8 


“9 a hs cs'35 
Te!23 Nigdest N 


“10 ODD PROTON 


LOGig Ty 


SS 
opD NEUTRON 
WEISSKOPF'S RELATION 


“ul 
FOR MI y-RAY TRANSITIONS 


-13 






~14 icles 8a 


10 20 40 60 100 200 400 600 1000 
y-RAY ENERGY (KEV.) 


Fic. 7. Lifetime vs. energy for M1 y-ray transitions. Where the standard deviations are 
larger than the dimensions of the points they are shown by vertical bars. 








GRAHAM AND BELL: MAGNETIC DIPOLEy-RAY TRANSITIONS 391 


which for M1 transitions becomes logie tr, = — 13.48 — 3 logy E (in Mev.), 
is shown as a solid straight line. Sachs and Ross (33) give a relation for M1 
transition probability (1/r,) which is proportional to a matrix element |u|? 
where » has the dimensions of a magnetic moment. For ordinary ‘“allowed” 
M1 transitions (shell model orbital Al = 0) u will be the magnetic moment of 
the odd neutron or proton and this relation becomes for the neutron case 
log tr, = — 13.39 — 3 log E and forthe proton case log r, = — 13.72 — 3 log E, 
which are in reasonable agreement with Weisskopf’s relation. 

Most of the experimental lifetimes are ~ 100 times longer than any of the 
theoretical estimates and may be called empirically forbidden. The transitions 
in odd neutron nuclei appear to obey an empirical relation log 7, = — 11.7 
— 3 log E while the three odd proton transitions indicate an empirical relation 
log ty = — 11.0 — 3 log EL, i.e. lifetimes about five times longer than in the odd 
neutron nuclei. These empirical relations are indicated as dashed lines in Fig. 7 
and must, of course, be used with caution in predicting M1 lifetimes in other 
nuclei. A similar plot of log [7,/(2 J, + 1)] against y-ray energy (not shown) 
did not improve the scatter of the experimental points from similarly drawn, 
E-* slope empirical lines. 

According to the shell model orbital assignments most of the empirically 
forbidden transitions are also theoretically forbidden, i.e. they involve a change 
of 2 in the orbital angular momentum / of the odd nucleon (A/ = 2). The 
transition in Li’ (proton orbital ~1;2 — p32) is both empirically and theoretically 
(Al = 0) allowed. Of interest are the two competing transitions in Hg!**. 
According to the shell model assignments (34) the 209 kev. transition is theoreti- 
cally allowed (neutron orbital p3;2 — 1/2) while the 50 kev. transition is theoreti- 
cally forbidden (neutron orbital p3;2 — d5,2). The relative transition probabilities 
are, however, very nearly the same after correcting for energy, and the two 
transitions are probably both empirically forbidden (see Fig. 7). The orbital 
assignments for the two (odd neutron) states in Fe®” (1,2 — 3,2) are the only 
ones which are consistent with the known decay scheme (36) and the shell model. 
This transition is thus another example of a theoretically allowed odd neutron 
transition which is empirically forbidden. It is interesting to note that the 
lifetimes of the two transitions in odd—odd nuclei are both distinctly shorter 
than any of the forbidden transitions. 
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THE SOUND FIELD OF A PISTON SOURCE! 
By ErRNEsT W. GuPTILL 


ABSTRACT 
An exact solution is presented for the sound field between two infinite walls 
when the source is a piston in one of the walls. A method is also outlined for the 
case of a source having any amplitude distribution which is a function only of 
radial distance from the center of the source. Several graphs are shown and they 
are intended to illustrate how the pressure on a receiving crystal would vary when 
the wall separation is changed. 


The sound field from an acoustical source set in a rigid and infinite wall has 
been investigated theoretically by King (3), Rayleigh (4), Stenzel (5), Williams 
(7, 8), and others. It seems quite improbable that an exact solution to this 
problem in terms of presently known functions is possible. Approximate solutions 
have’been obtained (5, 7, 8) and the velocity potential plotted for special cases. 
Not a great deal of emphasis, however, has been placed upon the difficulties facing 
the experimenter who attempts to measure the velocity of propagation with some 
precision. Williams and others (7) have given the deviation from the plane wave 
front by averaging the velocity potential over a circular area which would coin- 
cide with the receiving crystal. There is some question as to the usefulness of this 
procedure since it assumes that all the incident energy is absorbed by the receiver. 

The theory set forth here attempts to simulate more nearly experimental con- 
ditions by adding a second infinite and rigid wall (parallel to the first) in lieu of 
the receiver. The theoretical advantage of doing this is considerable since it makes 
possible an exact solution of the field between the two walls. In actual experiment 
the size of the ‘‘infinite’’ walls may be quite finite, depending on the directivity 


SOURCE 


Fic. 1. A schematic diagram of the circular source set in one of two parallel and rigid walls. 


1 Manuscript received October 6, 1952. 
Contribution from the Department of Physics, Dalhousie University, Halifax, N.S. 
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of the transmitter and the separation of the walls. Receiving may be effected by 
drilling a small hole in the second wall and placing a crystal one half wave length 
behind the first face of this wall. Preliminary tests indicate that this technique 
is very successful. 

Consider first the so-cailed piston source of radius a set in a rigid flange with a 
second flange a distance d in front of the source (Fig. 1). The velocity potential 
which is a solution of the wave equation is then subject to the following boundary 
conditions 


p> a, 
(1] =e pS; 
= 0O<p<o@ 


where p and z are the parameters of the cylindrical coordinate system. 

It is obvious that there will be waves travelling in the negative as well as the 
positive z-direction, and the integral solution of the wave equation given by King 
(3) must be altered to include both these waves, thus: 


[2] ¢ po a fla i 4. Boye | 


Jo(do) ddr 
1p 


where pp = V k?— )*, g < d, and k is a propagation constant with a positive 
imaginary part. The integration is conducted along the positive real axis but may 
be deformed into the complex plane of A, provided the contour remains on that 
branch of the function » which maintains the imaginary part of u positive. A (A) 
and B(d) are functions of the integration variable \ and must be adjusted so that 
the derivatives of ¢ satisfy the boundary conditions. If A (A) and B(A) are evalu- 
ated using the Fourier—Bessel transform (6) and Equation [1], one obtains for @ 
(the harmonic time variation will be omitted henceforth) 


me “cos u(z — d) 
[3] omen J 05 BE © 3,(0a) Juldp) ad 


which is a solution of the wave equation satisfying the given boundary conditions. 

The integrand in [3] is a single-valued function of \, and it will be shown that the 

integral may be evaluated by integrating along a properly chosen contour in the 

complex plane if the Bessel functions are expressed as a sum of Hankel functions. 

Two solutions of the integral are possible, one for each of the two substitutions: 
2Jo(Ap) = Ho(dp) + Ho(dp) 


and 


2J,(ka) = Hy(Aa) + Hi(da). 


The former will yield ¢ for p > a and the latter ¢ for p < a. In the first instance 
we have 


b= $f Pu) 10a) Hh) dd-+S fu) FiQa) Hi) dd 


where f(u) has been written for cos u(z — d)/u sin pd. 
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The second integral of Equation [4] may be written in terms of a Hankel 
function of the first kind for 


J, Js(na) Ha(dp) dr 2 ftw) Ji(da) HB (Ape) dd 
2 f fw) Ji (dae) Ho(Ape™) d(de™) 


; f f(u) Ji(da) Hi(Ap) dd 


and Equation [4] becomes 


[5] 6 =% [ Fu) S10) Hl) dx + $f Su) JiQa) FAQ) ad. 


The path of integration for the second integral of Equation [5] must be obtained 
from the path of the first integral by its reflection on the origin. The contour must 
then extend the entire length of the real axis. Since there is no contribution from 
the half circuit about the origin, the contour may be deformed into any path in 
the top half of the A-plane which does not cut the line of singularities occurring 
at ud = mz. Along this line » is wholly real and on crossing, the imaginary part 
changes sign. If the path is folded to the dotted path in Fig. 2, ¢:1/2m7 is the sum 
of the enclosed residues which are all at simple poles: 


[6] d= — ir >, cos nxz/d, 
n=0 


where A, = av 1 — (nx/kd)?. 
(The dash above the summation sign indicates henceforth that the ” = 0 
term is to be multiplied by 1/2.) 


A -PLANE 


BRANCH CUT 
For H> 


Fic. 2. 


The series converges rapidly for small values of kd but the convergence de- 
teriorates as kd becomes larger. For purposes of computation the series may be 
closed at that value of » which makes the argument of the Hankel function 
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imaginary, that is, when mm exceeds kd. Interpreted physically, this means that 
when the separation is less than one half wave length, only one of the family of 
cylindrical waves is propagated radially outwards. The others (m = 1, 2, etc.) 
are evanescent waves which attenuate rapidly with increasing p. When the wall 
separation increases to two half wave lengths, the cylindrical waves effectively 
propagated increase by one. 

If for 2J,(Aa) we now put Hi (Aa) + Hi (ha), a second integral for ¢ is obtained, 
namely, 
_2 


7] 2 =" J flu) Jodo) HiQa) dr + $ 


—o 


fw Jo(dp) Hi (da) dd. 


€ 


The Hankel function of order one has a singularity of the same order at the origin 
as well as a branch point. From this it follows that the contour must be indented 
upwards if it is to be closed in the top half of the A-plane (Fig. 3). Otherwise the 
contour is the same as for ¢;. On evaluation of the residues ¢2 becomes 


= ql CS F#E=—D _ LS cos nae /d Lenn Hides) | 
[3] d2 = of a imo cos nrz/d Xd ; 


A -PLANE 


BRANCH CUT 
' 
For Hy 


Fic. 3. 


The first term of Equation [8] is the velocity potential of a plane wave normally 
incident on a rigid wall. The series again represents a family of outward bound 
cylindrical waves with the attenuation of the evanescent ones increasing as ka 
increases. It is interesting to note the form of the ¢2 (which is proportional to 
the pressure) in the plane of the source, that is, at z = 0. If the wall separation 
is one quarter wave length, the plane wave component contributes nothing to 
this pressure, and the principal contribution from the series is the first term 
which gives a variation in the radial direction of Jo(kp). 

It is apparent from the form of Equations [6] and [8] that the former must be 
used for p > a and the latter for p < a, and it will now be shown that the two 
values of @ are identical when p = a. 

The Wronskian determinant of the two solutions Jo(z) and H}(z) of Bessel’s 
equation has the value 12/72, that is, 

Tidy) HQ) = Jo(dya) EE (dye) + 22 
TrA,a 
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which relation gives for ¢; at p = a, 


ores a ({ JoQxa) Hi Qua) as 
1 ima >, cos nxs/d| dd TAnAA,A 


n=0 


19} o| SB! _cosmas/d 5 cos nas/d Loa Hi) | 
ka t= (kd)” — (nm) n=0 And 
The first sum of [9] is the series expansion of [cos k(z — d)]/ka sin kd about the 
points kd = nz, and therefore the identity of ¢; and @2 at p = a is established. 
Equation [8] has been plotted for a few special cases. In Fig. 4, z has been put 
equal to d, ka = 10, p = 0, and \o/a| plotted for wall separations up to one wave 
length. It will be noticed that the velocity potential becomes large as kd ap- 
proaches nz. At this separation the height of the curve depends upon the imagi- 
_nary part of k, that is a, and in these curves it has been assumed that @ is much 
less than the real part of k. 


50 100 


WALL SEPARATION IN WAVELENGTHS 


Fic. 4. The piston source. A graph showing the modulus of velocity potential vs. wall 
separation for ka = 10, z = d, and p = 0 (Equation [8)}). 


The amplitude and frequency of the undulations occurring just beyond the 
half wave length separations are very sensitive functions of ka. As the value of 
ka increases, their amplitude and frequency increase, as may be seen in Fig. 5 
where the value of ka has been increased to 25 while the other parameters are the 
same as in Fig. 4. Apparently these undulations merge into a single smooth curve 
in the limiting case of ka > ~. 

In the case of a cylindrical boundary a solution with inward bound cylindrical 
waves is necessary and may be obtained by writing Equation [3] in terms of 
Hi(Ap) and Hi(da) instead of H}(Ap) and Hi(\a). The contours of Figs. 2 and 3 
must then be closed in the bottom half of the \-plane since the branch cut of the 
Hankel function of the second kind lies along the positive imaginary axis. 

The value of ¢2 averaged over a receiving surface of radius } at z = d is easily 
calculated: 
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1 b 
= +f $2 2mp dp 
rb 0 
1 


i RTT as UF Jub) Hi Qua) | 
2) 5— 5 . 
ate, (bd) — (nx) 


df sin kd 


WALL SEPARATION IN WAVELENGTHS 


Fic. 5. The piston source. A graph showing the modulus of velocity potential vs. wall 
separation for ka = 25, s = d, and p = 0 (Equation [8)). 
Fig. 6 shows a plot of |($2)/a| for ka = 25 and kb = 5. The points for this graph 
were easily computed. In fact, when ka > 10 and kb > ka/5 only one term in 
the series is important and even this may be omitted when kd is not near mr. 

The field between the walls may be calculated for a source having any ampli- 
tude distribution g(p) which is a function of p only. g(p) may be expanded in the 
range 0 < p < a thus: 


‘50 
WALL SEPARATION IN WAVELENGTHS 


Fic. 6. The piston source. A graph showing the modulus of velocity potential averaged over 
a surface which is coaxial with the source and which hasa radius 6; ka = 25, kb = 5,andz = 
(Equation [10]). 
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[11] g(e) = A,Jo(yip) + AzJo(y2p) +... 


where ya is the /™ root of the zero order Bessel function and A; a coefficient 
determined by (vide (1)) 


A,= J eo) Jo(vw) p dp. 


What this series yields for p > a is not relevant. We may now obtain the velocity 
potential for any term in the series by using the boundary conditions 


241) et 
(2 a aoa: 


= A, Jo(yp) p <a, 
(262) =0 O0<p<o. 
For the solution of ¢; one may proceed in the same manner as for the piston 
source. In place of Equation [3] we have (dropping the subscript / on y, ¢, and A) 


@¢=-AyaJ, (a) f f(u) 0lA@) Jo(Xo) dy 


and in place of [6] and [8] 
wee: ae iq — Jadu) Ho(dnp) 
$= — Biya) oy 408 88/4 Gay — Gon)” — Ge)? 
for p > a and 


[12] 


B sin Bd ea , (kd)? — (nw)? — (yd)? 
forp <a (s@= Vk — 7%). 


= A| Jove) cos B=) SR gq Jono) Hahn) | 


to 


° . 1:00 +50 


50 
WALL SEPARATION IN WAVELENGTHS 


Fic, 7. A graph showing the modulus of velocity potential from a source with an amplitude 
distribution equal to Jo(yp). ka = 10.2, ya = 5.52, p = 0, and z = d (Equation [12)). 
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Figs. 7 and 8 show a plot of the modulus and phase of ¢/a for ka = 10.2, 
p = 0, ya = 5.52, and z = d. Once again the height of the curves at the half wave 
position depends upon the imaginary part of k. One important difference between 
the graph in Fig. 7 and the corresponding case for the piston source (Fig. 4) is 
worth noting. In the former the half wave length points are always marked by a 
sharp rise in the velocity potential, whereas for the piston source the same places 
are marked by an outstanding maximum. 


0 Li) 25 


2-0 
WALL SEPARATION IN WAVELENGTHS 


Fic. 8. A graph showing the phase of the velocity potential. The modulus of this velocity 
potential is shown in Fig. 7. 


In Fig. 8 the slope of the curve is surprisingly linear for separations which are 
not close to an integral number of half wave lengths. The average slope is, how- 
ever, considerably less than that for a plane wave, the phase of which is repre- 
sented by the dotted line. This average slope gives a phase velocity very nearly 
equal to w/ V k*— 7°, and the distance between two consecutive points of equal 
phase is then [1 — (ya/ka)?]|-* times the free space wave length. If ya = 2.40 
< ka, the factor becomes 1 + 2.88/(ka)*. This result compares well with the 
approximate solution in (2), although it is doubtful whether the comparison 
means much since the boundary conditions in (2) were slightly different. 


The author is indebted to the University of Toronto’s Computation Centre 
which provided the computations used in Figs. 7 and 8. He also wishes to express 
appreciation for the benefit derived from frequent discussions with his colleague, 
Dr. A. D. MacDonald. 
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THE MAXIMUM AGE OF THE ELEMENTS AND THE 
AGE OF THE EARTH’S CRUST! 


By C. B. Cotiins, R. D. RUSSELL, AND R. M. FARQUHAR 


ABSTRACT 


Estimates are given of a maximum time of formation of the elements of 5.5 
billion years and of the time of formation of the earth's crust at 3.5 billion years. 
These estimates were based on the isotopic constitution of lead ores dated 
from isotopic analyses of radiogenic leads from uranium minerals. Methods of 
calculation are described. 

Experimental measurements of the isotopic constitution of lead from mass 
spectrograms of lead tetramethyl are reported with examples. 


A number of papers dealing with the age of the earth’s crust and setting a 
limit on the time of formation of the elements have been published in the past 
ten years. Most of these estimates have been based on variations in the isotopic 
constitution of leads which have been separated as lead ores from the earth’s 
crust at different times. These variations are believed to result from additions of 


100% 


Thorium 
Lead 


Uranium 
Lead 


450x10%yrs. 0.71x10°yrs, 13.9x10°yrs. (half life) 


Ordinary “= (2) 
Lead € om 


Uranium Thorium 
1.2 01.6% 
204 205 206 207 208 


Fic. 1. Common and radiogenic lead. (After Bullard.) 


radiogenic lead to ordinary lead. Ordinary lead is assumed to be a mixture of 
radiogenic lead from uranium and thorium combined with lead formed by other 
processes at the time of formation of the elements (see Fig. 1). 

Estimates of the age of the earth’s crust by Holmes (9), 3.35 billion years, 
and by Bullard and Stanley (4), 3.29 billion years, have been reported. Alpher 
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and Herman (2) assigned a value of 5.3 billion years ago as the earliest time of 
formation of the elements. These authors based their calculations on a set of 
isotopic analyses of common leads by Nier (12, 15). Nier also reported the ages 
of uranium minerals from various localities (13). These ages were used to date 
the time of mineralization of the common leads. 

All the investigators who have used Nier’s analyses have pointed out that 
more data were needed to avoid the possibility of gross errors in their calculations. 
We have analyzed a number of common leads from lead ores, some of which 
have been dated from Nier’s ages (13), and from isotopic analyses of radiogenic 
leads in our laboratory. New estimates of the age of the earth’s crust and the 
maximum limit on the time of formation of the elements have been calculated 
from our data combined with those reported previously by Nier (12, 15). 


EXPERIMENTAL 


The measurement of the isotopic constitution of lead (which is a solid at room 
temperature) has always presented a problem in mass spectrometer analysis. 
Since the usual procedure is to produce ions in a mass spectrometer by electron 
bombardment of gaseous atoms or molecules, ionization of solids has been 
carried out by raising the temperature of the solid until a vapor is given off 
which can be ionized in an electron beam like any other gas. Analysis of solids 
is generally much more difficult than that of gases. To avoid the problems 
associated with solid sample analysis recent investigations in the use of gaseous 
lead tetramethyl for the isotopic analysis of lead have been made by Collins 
et al. (5, 6). Dibeler and Mohler (7) reported similar studies of organometallic 
compounds including lead tetramethyl. Both groups arrived at the conclusion 
that the measurement of the isotopic abundances of lead from mass spectro- 
grams of lead tetramethy] is feasible. 

Lead tetramethyl dissociates on bombardment by 70 v. electrons to produce 
Pbt+, Pb(CH3)2+, and Pb(CH3)3+ ions. A low abundance of molecular Pb(CH3) 4* 
ions was also observed. Dibeler and Mohler (7) reported the abundances of all 
ions observed from lead tetramethyl from C+ to Pb(CHs)4*. Calculation of the 
abundances of the lead isotopes is complicated by the presence of hydrides, 
by the loss of from one to three hydrogen atoms, and by the C"* effect. For 
isotopic abundance measurements, spectrograms of the Pbt+ and Pb(CHs3)3+ 
ion groups are the most useful. 

The magnitude of the abundances of PbH* ions in the Pbt spectrogram, and 
of the C' effect in the Pb(CH;);* spectrogram given in this report are not the 
same as those reported by Dibeler and Mohler (7). Presumably the difference 
arises in the observed abundances of PbH+ ions because of different ion source 
conditions; their ion source was maintained at a temperature of 285°C., while 
ours was operated at room temperature. The validity of the analysis is not 
believed to be dependent on source conditions, providing that these conditions 
are maintained constant throughout the analysis. 

The interpretation of the Pb+ spectrogram is best illustrated by the following 


example. 
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PbH?* EFFECT IN THE Pbt SPECTROGRAM 





Mass number Ions 


204 Pb? 

205 Pb?4H 

206 Pb? 

207 | Ph? + Ph26H] 
208 | Pb8 + Pb207H 
209 | Ph28H] 


RELATIVE ISOTOPIC ABUNDANCES OF LEAD FROM THE Pbt SPECTROGRAM 





Relative 
abundance, % 


Relative ion 


Mass number intensity Lead isotope 





205 | 0.135 + 0.003 

206 20.31 +0.02 206 22.07 +0.03 
207 | 23.23 +0.02 207 23.33 +0.03 
208 50.65 + 0.03 208 53.02 + 0.05 


| 
| | 
—_ | 
| 

204 1.448 + 0.003 204 | 1.574 + 0.005 
| | 
209 | 4.23 +0.01 | | 
| | 


It will be noted that the PbH* ions present accounted for 8.3% of the total 
signal from the Pb*+ ion group. Furthermore, the fraction of PbH* ions present 
could be expected to vary depending on the temperature of the mass spectro- 
meter ion source (16) and the energy of the bombarding electrons, although 
under constant ionization conditions no change in this fraction was observed. 
No problem exists in the measurement of the isotopic abundances of ordinary 
lead where Pb?°’ is the most abundant isotope, and a good measure of Pb?°SH+ 
ions at mass number 209 is possible. But uranium lead has an abundant 206 
isotope, a low 207, and little or no lead 208. Here a satisfactory measure of the 
PbH?+ fraction is not possible and therefore radiogenic lead abundances calculat- 
ed from this spectrogram could be in error. These difficulties are largely avoided 
in the Pb(CH3)3+ spectrogram. 

In the interpretation of the Pb(CH3;)3* spectrogram, corrections have been 
made for the loss of one hydrogen atom, and for the C!* effect. Dibeler and 
Mohler (7) reported signals due to Pb(C3H10)* ions amounting to 0.33% of the 
total signal from the Pb(CH;3);* ions. This figure was obtained (by Dibeler and 
Mohler) by subtracting the contributions of ions containing C!* atoms and 
assigning the residue to Pb(C3H10)* ions. The relative abundances of the carbon 
isotopes vary according to the source of the carbon (11), so any correction 
based on the abundance of C' is uncertain. To account for our observed signal 
at mass 254 in the Pb(CH;);* spectrogram on the basis of C!* effect, an abun- 
dance of 1.125% for this isotope would be required. If the abundance of C'™ 
for carbon from limestone of 1.11% (14) is used, then the Pb(C3Hi0)+ ions 
would amount to about 0.08% of the Pb(C3;H,)+ ions. This figure applies to 
analysis when the mass spectrometer ion source is operated at room temperature. 
However, the calculation of the relative abundances of the lead isotopes is not 
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Mass number Ions 


| 

248 Pb?“(C2H;) o(C#H 2) 

249 Pb24(C?H3)3 

250 Pb2"(C#2H;).(C#H;) + Pb?(C?H;) (C#H») 
251 | Pb26(C12H)3 + Pb2°7(C!H;) o(C#H 2) 

252 | Pb?°7(C!2H;),; + Pb?(C?H;)2(C8H;) + Pb?°8(C?H;) o(C#H 2) 
253 Pb?°8(C#2H;)3 + Pb2°7(C!2H;) o(C3Hs) 

254 Pb?°8(C2H;) (CH) 


RELATIVE ISOTOPIC ABUNDANCES OF LEAD FROM THE Pb(CH3)3* SPECTROGRAM 


Relative 
Relative intensity Lead isotope abundance, “% 


| 
j 
| 
| 


Mass number 


248 0.01 +0.005 
249 1.505 + 0.005 204 1.569 + 0.005 
250 0.23 +0.01 
251 | 21.26 +0.02 206 | 21.97 +0.03 
252 | 23.56 +0.02 207 23.41 +0.03 
253 51.67 +0.03 208 53.04 +0.05 
254 1.77 +0.005 


The errors listed for Pb+ and Pb(CHs)st spectrograms described above apply to repro- 
ducibility of the measurements; the absolute error is believed to be less than 1%. 


changed by the presence of a small fraction of Pb(C3Hi»)* ions whose contribu- 
tion may be combined with that of C'* and measured from the ion current at 
mass number 254. 

Either the Pb+ or Pb(CHs;)3* spectrogram is suitable for measuring the 
abundances of the isotopes of common lead. The Pbt spectrogram is probably 
the less reliable because of Hg?** which adds to the Pb®** peak, and because the 
PbH+/Pb+ fraction is dependent on the energy of the ionizing electron beam 
and the temperature of the ion source. In the analysis of radiogenic lead made 
up chiefly of Pb?°* with a low Pb” and little or no Pb*** or Pb, no direct 
measure of the PbH* contribution is possible. In the Pb(CH3)3+ spectrogram 
the largest correction is for the contribution of ions containing a C!* atom, whose 
magnitude is not critically dependent on the energy of the ionizing electrons or 
the temperature of the mass spectrometer ion source. The contribution from 
ions which have lost one atom of hydrogen is of the order of 1%, and any small 
variation in this quantity would not change the calculated abundances of the 
lead isotopes by an appreciable amount. For these reasons the isotopic constitu- 
tion of uranium lead determined from the Pb(CH3;);+ spectrogram of lead 
tetramethyl is believed to be reliable. Spectrograms of lead tetramethyl! are 
‘shown in Figs. 2 and 3. By analyzing both ordinary and radiogenic lead it was 
possible to check the various isotopic contributions which gave rise to the ion 
currents at the different mass numbers. 

Lead was extracted from uranium samples by standard chemical methods, 
and the preparation of lead tetramethyl was carried out with the assistance of 
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. Normal Lead _Tetramethyl_ Spectrogram 
Pb(CHs)5 SS ee ea ete ee ae 


Sample source: Ethyl Corporation, Detroit 


Pb(CH;),* 


Pb(CH3)2* 


ar on g ses 
uw RANA 8&8 


x 
N 
Fic. 2. Normal lead tetramethyl spectrogram. 


Pb(CH3)3* Radiogenic Lead _Tetramethyl Spectrogram 
] Congo _pitchblende 


+ 
Pb(CH3)) 





Pb (CH3)>° 











Fic. 3. Radiogenic lead tetramethyl spectrogram. 
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Prof. G. F Wright and Mr. H. L. Cohen of the Department of Chemistry of the 
University of Toronto. The method of preparation has been reported (5). 
Our first preparations required about 50 mgm. of lead chloride, but a revised 
method now being investigated using less than 10 mgm. of lead chloride has 
already yielded adequate samples of lead tetramethyl for mass spectrometer 
analysis. 


ISOTOPIC CONSTITUTION OF RADIOGENIC AND ORDINARY LEADS 


The relative abundances of, the isotopes of common leads from various 
localities have been determined by the analysis techniques described in the 
preceding section. A number of these leads have been assigned ages from dated 
uranium minerals which occurred in the same geographical and geological 
areas, and which are believed to have been deposited at approximately the same 
time as the common leads. 

All the ages were determined from the ratio of radiogenic Pb**? to Pb?°* for 
lead extracted from uranium minerals, according to the well-known relation 
used by Nier (13) 

N(Pb7) = N(U?%) (exp A(U2*)¢ — 1) 


N(Pb?*) — N(U238) (exp A(U248)t — 1)° 


N(Pb**’), etc., refer to the number of atoms of the indicated isotope. \(U?**) 
and (U8) are the decay constants of the uranium isotopes. The ratio of the 
abundances of the uranium isotopes at the present time is 1/139. Numerical 
values for the ratio of Pb**’ to Pb?°* were taken from the isotopic analysis of 
uranium lead after correcting for the presence of ordinary lead on the basis of 
the abundance of Pb?%, which is not derived from uranium or thorium. This 
correction was small as may be seen in Table I, so that any uncertainty in the 
assumed isotopic constitution of ordinary lead would introduce a negligibly 
small error. 

In Table I, both our determinations and Nier’s have been listed for samples 
of uranium bearing minerals from Wilberforce, Katanga, and Great Bear Lake. 
The uraninite from Wilberforce, Ont., was exceptionally fine unaltered material 
free from ordinary lead; the agreement in age values is within 1%. Pitchblendes 
and uranium minerals formed by alteration are less suitable for age determina- 
tions. For seven samples from Katanga, Nier reported values from 620 to 665 
million years, as compared to a spread of from 615 to 655 million years for five 
samples analyzed in this laboratory. Nier’s single value for Great Bear Lake 
pitchblende falls within our age spread for samples from the same deposits. 
These comparisons provide a rough check between our values and Nier’s; 
they also show the spread in values that can be expected for determinations made 
on uraninites and on the less reliable pitchblende ores. 

All the ages listed in Table I have been used in an endeavor to date the 
ordinary lead samples whose isotopic analyses appear in Table II. Where Nier’s 
analyses of ordinary leads have been quoted, his sample numbers have been 
given. As the main source of uncertainty in calculations based on the isotopic 
constitution of lead ores is in the assignment of the age or time of deposition of 
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the ore, all Nier’s specimens from one locality believed to be the same age have 
been combined, as was done by Bullard and Stanley (4), to avoid giving excessive 
weight to lead from any one locality. The ages assigned by Bullard and Stanley 
to Nier’s leads of post-Cambrian age have been listed without change, and with 
them has been included sample No. 5 from Bolivia, believed like most deposits 
there to be Tertiary in age. Hence it has been given the same approximate 
mid-Tertiary age of 25 million years as sample No. 1 from the Casapalca Mine, 


Peru. 
TABLE II 
RELATIVE ISOTOPIC ABUNDANCES OF LEADS FROM LEAD ORES OF KNOWN AGE 


u : Isotopic abundances _ 
Geological (Pb? = 1,00) 
Nier’s |Our Source Age0N, eee eee et 
Nos. | No. 10° years 206 207 208 


Casapalca Mine, Peru 18.76 15.55 38.39 
Arizona, U.S.A. : 18.81 15. 38.63 
Freiberg, Saxony é 18.07 15.40 38.0 
Durango, Mexico : $97). 35. 38.5 
Bolivia ! 18.94 15. 39.70 
Washington, Idaho, U.S.A. j 17.67 42 37.38 
Saxony, etc. 220 18.01 5. 37.89 


16, 17 
11, 23 
24 
21 


9,10 
5, 6, 12, 
18, 19. 


Nook whe 


Yancey Co., N.C., U.S.A. 220 18.43 f 38.2 
Franklin, N.J., U.S.A. 400 17.15 ; 36.53 
Katanga, Belgian Congo 630 i: 38.34 
Tetraulte Mine, Que., Canada 1030 16. : 35.6 
Frontenac Co., Ont., Canada 1030 16.75 : 36.49 
Broken Hill, N.S. Wales 1500 16. .385 35.38 
Great Bear Lake, N.W.T., Canada 1400 § ; 35.3 
Negus Mine, Great Slave Lake, 1850 6 , 34. 
N.W.T., Canada 
Rex, Stark Lake, N.W.T., Canada 1850 : 5.65 36.34 
Sioux Lookout, Ont., Canada 2480 : a 33. 








Two methods have been used to date the Pre-Cambrian specimens. The com- 
mon lead samples from Katanga and Great Bear Lake came from veins or mines 
containing uranium ores, and the ages determined for the latter have been used. 
In the other cases no dates are known for the deposits themselves and it has been 
necessary to assume that the deposits are of approximately the same age as 
the geological provinces in which they occur. It has been pointed out for several 
Pre-Cambrian shield areas that patterns of foliation, major fault zones, changes 
in rock types, and kinds of mineralization enable different provinces to be 
distinguished. It has also been found for those provinces within which there are 
many localities dated by radioactive minerals that all the pegmatite deposits 
have ages lying within a range of a very few hundreds of millions of years. 
Although some mineralization has undoubtedly occurred at later dates, this 
method probably roughly fixes the age of the majority of ores of a province. 

Sample No. 11, Table II, from the Tetraulte Mine, Quebec, and sample No. 
12, from Frontenac Co., Ontario, are considered to be in the Grenville geological 
province whose best dated locality is Wilberforce, Ont., at 1030 million years. 
Sample No. 13 was dated from the Radium Hill age of 1500 million years (No. 
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22, Table 1). The value of 1850 million years, representing the maximum age 
of the uranium samples from Goldfields, Lake Athabaska, and from the Rex 
property on Great Slave Lake, has been taken as dating ordinary lead samples 
Nos. 15 and 16. The age value of 2480 million years for No. 17 has been estimated 
by averaging the two values for the Huron Claim reported by Nier with a value 
determined by Ahrens (1) using the rubidium-strontium method. 

The relative abundances of the isotopes of common leads from ores of un- 
known age are given in Table III. Several of these leads have isotopic constitu- 
tions quite different from any previously reported. They will be discussed in the 
last section of this paper. 


TABLE III 
RELATIVE ISOTOPIC ABUNDANCES OF COMMON LEADS OF UNKNOWN AGE 


Isotopic abundances 
Source (Pb? = 1.00) 

Source ni a a TE ie 
206 207 208 


j 
| 
| 





Ivigtut, Greenland 14.65 14.65 
Horseshoe Island, Yellowknife, N.W.T., 

Canada 14.16 .08 
Upper Canada Mine, Ont., Canada | 19.60 .22 
Rosetta Mine, South Africa 12.65 ae 
Box Mine, Goldfields, Sask., Canada 20.88 74 
Matachewan, Ont., Canada | 15.83 41 
Delhi Twp., Ont., Canada | 16.07 47 
Garson Mine, Sudbury, Ont., Canada | 23.35 
Worthington Mine, Sudbury, Ont., | 26.00 

Canada | 





A LIMIT ON THE TIME OF FORMATION OF THE ELEMENTS 

At the time Nier (12, 15) reported the isotopic constitution of 25 lead ores, 
he pointed out that the abundances of Pb?°*, Pb?°’, and Pb? relative to Pb?®4 
tended to be smaller in older ores, and suggested that this variation might be 
attributed to contributions of radiogenic lead from uranium and _ thorium. 
Later Alpher and Herman (2), by using certain simplifying assumptions, were 
able to deduce an approximate mathematical representation of this effect, and 
to evaluate the parameters involved. 

To permit such a calculation to be made, it must be assumed that since the 
time of formation of the elements lead produced from the decay of uranium 
and thorium has been continuously added to the original lead. It is further 
assumed that on the average the relative amounts of lead, uranium, and thorium 
in the earth have remained the same everywhere except for the change due to the 
steady decay of uranium and thorium to lead. At the time of a lead mineraliza- 
tion a sample of lead was concentrated from a body of rock. The resulting lead 

‘ores contain negligible amounts of uranium and thorium and so constitute a 
permanent record of the isotopic constitution of lead in the earth’s crust at the 
time of formation of the mineral. 

If these assumptions are true, the isotopic composition of lead at /, years 
before the present can be represented by the equations 
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Xm =a — aV(e™ — 1), 
[1] Ym =b— Vie ™—1), 


and Zn =O — we ™ — 1). 

Xm; Vm, and Z, are the abundances of Pb?%*, Pb2°7, and Pb? relative to Pb? 
at time f,,, and a, 6, and ¢ are the corresponding abundances of the respective 
lead isotopes in the earth’s crust at the present time. V and W are the present 
relative abundances of U?** and Th?%? with respect to Pb®°4. a is the present 
ratio of U?8* to U?*, and X, A’, and X” are the decay constants of U2#8, U2%5, and 
Th?* respectively. According to our assumptions, V and W are not functions 
of time or of the particular sample concerned. 

A least squares analysis has been used to calculate the values of a, b, c, V, 
and W which give the closest agreement between equations [1] and the isotopic 
constitution of the lead ores of known age listed in Table II. This was done in 
the following way. 

Define 

&6= DY {x, —a+aV(e™ —1)}? 


m 


+ Dlm—-b+ Ve *™-1)}" 


m 


+ > {en —c+ we ™ — 1)}*. 


For the best fit, 6 is a minimum, and therefore 


% _ 8 _ 0 _ % _ % _o 
da db dc OV OW 
The application of these conditions results in five independent linear equations 
in five unknowns: 


n 


>, Xm — na +aVy, (e™ —1) 


m=1 m=1 


4s CE ~2 


m m 


Hs —a+ WHC *—1) =0, 


m m 


ad xm(™ —1) + SS ymle®™ — 1) — aad (e™ — 1) 


m m m 


~ bo 6 * 14+ VE We = 17 4 = a ee 


m m 


Y mle — 1) - ee - 1) +WOe'- 1) =0. 


Using A = 0.154 X (10° years)-', 
’ = 0.980 X (10° years)~!, 
dN” = 0.0499 X (10° vears)—!, 

= 139, 
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and the isotopic abundances listed in Table II, it was found that 
= 18.45, 
15.61, 
38.40, 
0.0714, 
= 38.2(8). 


These values are for the present time. On substituting the above numerical 
values in equations [1], it was then possible to plot the curves shown in Fig. 4. 
These curves show the mean isotopic composition of lead distributed throughout 
the earth’s crust as a function of the time before the present. 

The mathematical calculations give no unique starting points for the curves 
shown in Fig. 4, and therefore no unique value for the time at which uranium 
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Fic. 4. Lead isotopic abundances relative to Pb?” versus time. These curves were obtained 
from a least squares fit of the data in Table II. 
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first began contributing to the lead isotopes. However, if we assume that 
equations [1] can be used to extrapolate the curves backward in time (as shown 
in the dotted portions of the curves of Fig. 4) a limit on the time of formation 
of the elements can be taken as the time at which the abundance of one of the 
three lead isotopes becomes negative, for this surely represents an impossible 
physical situation. It can be shown from either Fig. 4 or equations [1] that the 
abundance of Pb**? becomes zero 5.5 X 10° years ago. This figure may be con- 
sidered a limit to the time of formation of the elements. 

The figure of 5.5 billion years given as a limit to the time of formation of the 
elements may be in error owing to the use of incorrect physical constants, 
owing to the relatively few points used for determining the curves, or owing 
to a failure of the basic assumptions. The present uncertainty in the relative 
abundances of the uranium isotopes, U?38/U2%, is about 1%. The uncertainty 
in the decay constant of U?** is of the order of 1.5%, while the uncertainty in 
the U?*8 decay constant is.only about 0.25%. It is fortunate that the position 
and intercepts of the Pb?’ curve (Fig. 4) are less sensitive to changes in para- 
meters than those of the other two. This effect can be seen by comparison of the 
curves reported by Alpher and Herman (2). The data used to calculate the curves 
were appreciably different, as was the decay constant of U*** for which Alpher 
and Herman used 0.920 X (10° years)~!, and although the Pb?°* and Pb?° 
curves shown in their paper are quite different from those of Fig. 4, the Pb?°7 
curve is almost the same. Gross errors arising from the failure of basic assump- 
tions will be discussed in the last section of this paper. 

THE AGE OF THE EARTH’S CRUST 

In 1946 Holmes (9) suggested an ingenious method for estimating the age of 
the earth’s crust from the isotopic constitution of lead from lead ores of known 
age. More recently Bullard and Stanley (4) have modified Holmes’ method, 
treating the experimental data by a least squares analysis rather than by 
Holmes’ more laborious approach. Both authors used Nier’s isotopic analyses 
of common leads (12, 15), and reported values close to 3.3 billion years (Holmes 
—3.35 billion; Bullard and Stanley —3.29 billion years) for the age of the earth’s 
crust. This calculation has been repeated by the present authors using the least 
squares method of Bullard and Stanley (4), and basing the calculation on the 
experimental data given in Table II. These calculations all depend on the 
behavior of lead, uranium, and thorium in the outer part of the earth as set 
down by Holmes (9). 

Holmes’ method depends on several basic assumptions. It is assumed that the 
materials now found in the outer parts of the earth were initially so thoroughly 
mixed that the four isotopes of lead were present in the same relative proportions 
everywhere. When the crust solidified, physical and chemical processes produced 
small differences in the relative amounts of lead, uranium, and thorium from 
place to place. The lead which crystallized in the earth’s crust at that time has 
been called primeval lead. As time passed, part of the uranium and thorium 
decayed to lead, which added to primeval lead, and changed its isotopic con- 
stitution by an amount depending on the relative proportions of lead, uranium, 
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and thorium present in that particular part of the crust. At a later time a sample 
of lead was concentrated from a body of rock to form a lead mineral containing 
a negligible amount of uranium and thorium. This lead provided a permanent 
record of the isotopic constitution of lead in the crust at that time. 

In the ideal case, the isotopic constitution of lead in the crust at time ft, 
(and hence the isotopic composition of a lead mineral formed at time ¢,) can 
be represented by the mathematical expressions 


, Alo Atm 
Xm = Ao taVy(e®? — &”) 


; r 7d te Mtn 
and Vu bo + Vnn(e ==¢ ). 


where x,, and vy» respectively are the abundances of Pb?°* and Pb?® relative to 

Pb?°4 at ft, years before the present, a) and do are the corresponding abundances 

of primeval lead, fo is the time of solidification of the crust, and V,, is the abur- 

dance of U?*® with respect to Pb***. a is the present ratio of U*** to U**®, and A 

and }’ are the decay constants of U*** and U*** respectively. It is assumed that 

V,, is a function of the particular lead mineral, but is not a function of time. 
Equations [2] can be combined to obtain 


dt 
Xm — a ale ° 


[3] Ratt 


Ym — bo (e 
The unknown parameters do, }9, and f9 must be chosen so that equation [3] 
most nearly fits the experimental data in Table II. This is done in the following 
way. 


n ’ , 
Define = > Wm {a (Vm pe bo) (es" =e gt) a (Xm Ss ao) (e bei e 3 


m=1 


° 
so bis Wm €m 


m 


where the summation extends over the 2 samples. w,, is a weight that is calculated 
for each sample. As Bullard and Stanley (4) pointed out, the use of properly 
calculated w, provides unique solutions for do, bo, and ¢) that are independent of 
the manner in which 6 is defined, and that are the most probable solutions. 
For this calculation 


/ 2 r 2 7 2 
Wm = 1 (Xx m + Vn + Lx Ra 
> yz, K” fas 
where eo — (e * —e ) ixm, 


, lo Alm 
Yen a(e eg ) 8Ym, 


and Tm = {X'(%m — ao) os ad (Ym — bo) €'"} dtm, 


and where 6x, 5m, and dtm are the uncertainties in Xm, Vm, and tm. 

For the present calculation, 6x, and dym were taken to be + 0.1 in all cases, 
and dt, was taken as + 0.2 X 10° years in all but four cases. Samples Nos. 11, 
12, and 14 (Table II) were believed to be better dated than the others. For these 
samples dtm was taken to be + 0.15, + 0.10, and + 0.10, respectively, all in 
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10° years. The age of sample No. 15 was assigned an uncertainty of + 0.30 10° 
years. The assignment of these uncertainties is necessary in the least squares 
treatment of the data to avoid giving undue weight to younger samples.: This 
procedure can be justified physically on the basis of uncertainties in the isotopic 
abundance measurements, and in the assignment of ages to lead ores from 
dated uranium minerals supposedly deposited at the same time. 
€m is not linear in the unknowns. However, if approximate values do, b, and 
to are chosen, we can say that 
Oe Oz 
e2i+— Ant+=— 
+ Odo ‘ bo 
which is linear in Aap, Abo, and Ato. Where Ady = ay — G, Abp = by — 5, 
Aty = to — iy, and the derivatives are evaluated for the approximate solutions. 
Then 
f Oe Oe Oe \ 
6 = Wmye + =: Aa +—: AbD+ —: Ato( - 
Le mye t aq * 4a + 95° Ab + ay, Atos 


m 


0€ 
Abo + al, Ato 


For the best fit between the observed values and the theoretical equation 6 is a 
minimum and 


‘ ee Le 
Mi] (Ado) (Abo) (Ato) 
Using the data in Table II and taking 
a = 139, 
A = 0.154 X 10-° years“, 
and ’ = 0.980 X 10-° years“, 
equations [4] were solved to give 
ao = 11:38, 
by = 13.55, 
to 3.5 X 10° years. 
DISCUSSION 
The possibility of gross errors in the results which have been presented must 
be considered. Our value of 5.5 billion years as a limit on the time of formation 
of the uranium isotopes was obtained by extrapolation from the isotopic 
constitution of common leads, the oldest of which was dated at 2.5 billion years. 
Such an extrapolation is valid only if the relative amounts of lead, uranium, and 
thorium in the earth’s crust have remained the same, except for radioactive 
decay, from the time of formation of the elements to the present. This difficulty 
has already been pointed out by Alpher and Herman (2) who wrote as follows: 
“It is generally believed that there was a marked increase in the amounts of 
uranium and thorium relative to lead in the earth’s crust at the time of its 
formation. This implies the possibility of error in an extrapolation of the curves 
shown in Fig. 1 [our corresponding curves are shown in Fig. 4] to times earlier 
than that of crustal formation.” 
Our value of 3.5 billion years for the time of formation of the earth's crust is 
slightly higher than the 3.29 billion figure of Bullard and Stanley (4) and the 
3.35 billion value of Holmes (9), but the change is not large in view of the 
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additional data used in our calculation. From this comparison of estimates 
by different investigators using different methods of calculation and in: this 
paper additional experimental data, it seems probable that an answer of the 
order of 3.5 billion years is likely to be obtained as long as the calculation is 
based on the fundamental assumptions of Holmes’ method. Alpher and Herman 
(2) expressed doubts about the validity of Holmes’ method, but agreed that his 
figure of 3.35 billion years for the time of formation of the earth’s crust was 
reasonable in view of Behr’s (3) figure of 3.8 billion years from nebular recession 
data for the start of the expansion. Bullard and Stanley (4) pointed out that a 
. likely cause of gross error in Holmes’ method was a failure of the basic assump- 
tion that the rocks were isolated when the crust solidified, and then developed 
undisturbed until the time of deposition of the lead ore. The anomalous isotopic 
constitutions of the Sudbury leads (Nos. 25 and 26, Table III) support this 
criticism, although it is not easy to estimate the effect of such samples on the 
age of the earth’s crust calculation. 

Holmes assumed that all lead everywhere had the same isotopic constitution 
at the time of solidification of the earth’s crust. He called this lead primeval 
lead. According to this argument, no lead found at any subsequent time could 
have lower isotopic abundances relative to the 204 isotope than primeval lead. 
We have reported one common lead from the Rosetta Mine, South Africa, 
listed in Table IV, whose isotopic abundances fall within the limits set on the 
calculated isotopic constitution of primeval lead as shown in Table IV. Primeval 
lead must then have isotopic abundances relative to Pb®** not greater than those 
listed for the Rosetta sample. 

TABLE IV 
ISOTOPIC CONSTITUTION OF PRIMEVAL LEAD 





Isotopic abundances relative 
>p204 = 
Sample to Pb? = 1.00 


206 207 208 





Rosetta Mine, South Africa 12.65 14.27 32.78 
Primeval lead, calculated by Holmes (9) 10.945 13.51 
+8to10%+8to 10% 
Primeval lead, calculated by Bullard and 11.86 13.86 
Stanley (4) + 1.47 + 0.56 
Primeval lead, calculated by the present authors 11.33 13.55 
Isotopic constitution of lead from Fig. 4 at 11.35 13.46 31.1 
t = 3.5 X 10° years 





In Table IV are listed the abundances of the lead 206 and 207 isotopes 
calculated by the present authors for primeval lead at the time of formation of 
the earth’s crust, 3.5 billion years before the present. A comparison of the isotopic 
abundances of the same isotopes has been made by taking values from the curves 
of Fig. 4 at a time 3.5 billion years before the present. The good agreement 
shows that the primeval lead abundances derived from the age of the earth’s 
crust calculation fall on the extrapolated average curves drawn from the inde- 
pendent limit on the time of formation of the elements calculation. This check 
indicates that the two independent calculations are consistent with each other, 
but does not prove that either or both are valid. 
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The calculations carried out in this paper have been based on a set of leads 
of known isotopic constitution taken from lead ores of known age. While the 
two calculations are different in principle, they both depend on the total effect 
of all the determinations rather than on any individual one. Single samples such 
as Rosetta Mine, No. 22, Table III, can not be safely dated from the average 
curves of Fig. 4, or by subtracting the primeval lead abundances (which have 
rather large errors) from the isotopic abundances of the sample and then 
considering the Pb?°’?/Pb?"6 ratio of the remainder. This difficulty arises because 
of the varying amounts of lead, uranium, and thorium present in any one 
locality at different times in the earth’s history, a fact that is well illustrated 
by the isotopic constitution of several of the leads in Table III which do not 
fit the curves of Fig. 4. A recent attempt was made by McCrady (10) to set up 
a time scale from which the age of ordinary leads could be obtained from relative 
isotopic abundance versus time curves, on the assumption that the relative 
amounts of uranium, thorium, and lead have remained constant (except for 
radioactive decay) everywhere throughout the earth’s history. This assumption 
does not fit the experimental data given in Table III, as has already been pointed 
out in the preceding discussion. In the authors’ opinion, no general method of 
dating the time of mineralization of individual lead ores from their isotopic 
constitution is likely to give more than very approximate ages. But investigations 
based on a number of isotopic analyses of leads from different localities in any 
particular geological province would almost certainly provide a sound basis for 
fixing the time of mineralization of the lead. Such studies are also likely to 


provide specific information about the origin of lead ores, and about the chemical 
and physical processes by which galena is concentrated as a mineral in nature. 
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FISSION YIELDS OF THE STABLE AND LONG-LIVED 
ISOTOPES OF CESIUM, RUBIDIUM, AND STRONTIUM AND 
NUCLEAR SHELL STRUCTURE! 


By D. R. Wies,? B. W. Smitu,? R. Hors_ey,’ aNp H. G. THODE 


ABSTRACT 
The relative yields of the isotopes of cesium, rubidium, and strontium from 
thermal neutron fission of U** have been determined mass spectrometrically. 
The cesium isotope yields are combined with those obtained previously for the 
xenon isotopes to give high precision yields for mass chains from 131 to 137. 
In this work neutron capture reactions have been considered and corrections 
made where these take place to an appreciable extent. The results give further 
evidence of abnormal yields in the 82 neutron shell region. The half-life of Cs!7 

was determined and found to be 33 + 2 years. 


INTRODUCTION 

Previous mass spectrometer investigations (14, 8) of the stable and long- 
lived isotopes of xenon and krypton, formed from the fission of U***, revealed 
for the first time the existence of ‘‘fine structure’ in the mass yield curve. 
Kr**, Xe!88, and Xe!*4 were found to have abnormally high yields—Xe!*4 
having a yield about 30% above the smooth curve drawn through adjacent 
points. These high yields are undoubtedly due to the inass chains in question 
having members with 50 and 82 neutrons. It should be pointed out that “‘fine 
structure” in the Kr*‘ mass range is not certain, in view of the steepness of the 
yield curve at this point and the relatively low yields. Glendenin (3) has suggest- 
ed a chain branching mechanism to explain these yields. He has postulated 
that nuclides with 83 neutrons, one more than a closed shell, will have a certain 
probability of emitting neutrons when formed as a primary fission product. 
This mechanism does ‘explain qualitatively some of the yield values obtained 
between masses 131 and 136. The fission yields of the cesium, rubidium, and 
strontium isotopes, which fall in mass chains overlapping and adjacent to the 
mass chains of xenon and krypton studied previously, are of considerable 
interest. If the abnormal yields of certain mass chains are due completely or in 
part to chain branching, then adjacent chains will also have abnormal yields 
(vields differing from the interpolated values of the smooth fission yield curve). 
The fission chains in the two mass regions being studied are given in Tables | 
and I]. It is seen that mass chains 133, 135, and 137 end in stable or long-lived 
isotopes of cesium, and mass chains 85, 87, 88, and 90 end in stable or long-lived 
isotopes of rubidium and strontium respectively. 

Inghram, Hess, and Reynolds have already reported (7) on the abundances 
of the cesium isotopes from fission measured mass spectrometrically. Their results 
indicated high relative yields for Cs'** and Cs'*’. If their relative abundance data 
for Cs'*3 and Cs!*7 are normalized to 7.43% for the 133 mass chain yield determined 
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TABLE I 
FISSION PRODUCT CHAINS AND CHAIN BRANCHING IN THE HIGHER MASS RANGE 
NUCLEAR CHARGE 


50 51 52 53 54 55 56 
Sn™!___» Sh3!_.5, Tel81_5, ]18!__y, Xela (Stable) 





Sn Sb!3____», Te!3_» ]82__y Xe1® (Stable) 
(a) - 


| 
Sn'#___»Sb!85__» Te!#3____» ] '583_» Xe!83____y (5188 (Stable) 
(b) An 
Sb'31___» Te!44___-» ]181____y, Ye184 (Stable) 
(c) An 





alma citi < Muara rea yo Oe (Stable) 
( r ’ yr. 


Sb'95___ Te!#®___y 1186___y Xe1% (Stable) 
(e) AAn (g)** Delayed 
Xel37 
A 
/ 
Te!37___» [137 Xe!7—Primary yield 
i 
Xe!87___» Cs137___ Ba'37 (Stable) 
_ A(33 yr.) 
(f) | n 
Te!38__» ]38_ Xel38____(C5188_y Ba188 (Stable) 





* Vertical arrows lettered (a) to (f) indicate the branching of a chain by instantaneous emission 
of a neutron from a primary fission product having 83 neutrons. 
** Delayed neutron activity (g) is 0.2% of fission. 
TABLE II 
FISSION PRODUCT CHAINS AND CHAIN BRANCHING IN THE LOWER MASS RANGE 
NUCLEAR CHARGE 
34 35 36 37 38 39 40 
Se**_—_» Br**——_»Kr* (Stable) 
An* 
| 
| 
Se®__.»Br®_»Kr®_____» Rb® (Stable) 
An + 
[A Kr88 — 
| 
Se*__-+»Br*—_»Kr*® (Stable) Rb**——-»Sr® (Stable) 
Delayed An 


Br®?___» Kr*" Rb’?___» Sr87 (Stable) 
el (6.0 X 10!yr.) 


*» Kr87__ An 


Br®8. -»>Kr®&§_____» Rb®&____» Sr88 (Stable) 
An 
Kr89_____»Rb*___»Sr#*___» Y*®9 (Stable) 
An 


ee (Stable) 
9.9 yr. 


*Vertical arrows indicate the branching of a chain by emission of a neutron from a primary 
fission product having 51 neutrons. 5 
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in this laboratory for Xe!**, then the yields of these chains become 7.43, 9.18, and 
6.94% respectively (8, 1). (See footnote on page 423.) 

This high yield indicated for Cs'** cannot be satisfactorily explained by the 
Glendenin mechanism and it does not seem to be consistent with other fission 
yield data. Because of this, and because Inghram, Hess, and Reynolds were 
unable in their work to correct for the presence of normal Cs!** contamination, 
which may have been introduced in their extraction process, we have repeated 
their work on a number of different samples of uranium irradiated under different 
conditions. In addition we have measured the Rb**/Rb*’ and Sr**/Sr® ratios 
in fission product rubidium and strontium to determine whether there is any 
evidence of abnormal yields in this mass range. We have also had occasion to 
determine the half-life of Cs!” using the mass spectrometer abundance data from 
four samples of different ages. 

EXPERIMENTAL 
Cesium and Rubidium 

The procedure for the simultaneous isolation of carrier-free cesium and 
rubidium from irradiated uranium was as follows. The uranium was dissolved 
in nitric acid (6N) and the solution was evaporated to crystals. The crystals 
were dissolved in diethyl ether and the fission products separated from this 
solution by several extractions with small portions of dilute nitric acid. Ferric 
ion was added to the aqueous portion and precipitated as the hydroxide with 
ammonia gas, the supernatant liquid containing the rubidium, cesium, barium, 
and strontium and void of all the hydroxide forming elements. To the supernatant 
liquid was added ammonium ion and the subsequent precipitation of the 
ammonium chloroplatinate carried down only the cesium and rubidium. The 
supernatant solution at this point was used for subsequent isolation and purifica- 
tion of strontium. The chloroplatinate crystals were dissolved in hot water and 
the platinum was removed by reduction with formic acid. The carrier-free 
cesium and rubidium were then converted to the nitrate by repeated evaporation 
with nitric acid. The resultant nitrates were in a satisfactory form for mass 
spectrometer analysis. 

In order to determine the amount of contamination of natural cesium and 
rubidium a blank was prepared in the following manner. Two disks of equal 
weight were cut from an uranium rod and after irradiating one of them, both 
were dissolved in equal volumes of nitric acid. Half of the volume of one of 
these samples was combined with half of the volume of the other to make a 
synthetic sample having twice as much stable rubidium and cesium contamin- 
ation as the original active solution. 

The simultaneous separation of rubidium and cesium was then carried 
through with both samples as in the procedure outlined above. The comparison 
of the abundances of the isotopes in the two samples makes possible the deter- 
mination of the fission product isotope abundances. 


Strontium 
The isolation of carrier-free strontium was completed by the following process. 
The strontium solution, after the chloroplatinate precipitation just described, 
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was evaporated to dryness and fumed. Contaminants at this point were silica, 
platinic chloride, some organic matter, and small amounts of iron and uranium. 
The residue was dissolved in dilute hydrochloric acid, and 2 mgm. of ferric ion 
were added. This was precipitated as ferric hydroxide at pH 9 with ammonium 
carbonate. This caused the quantitative sorption of strontium carbonate and/or 
strontium hydroxide. In this way the desired component was separated from 
all substances soluble in basic solutions. 


Mass Spectrometer Analysis 

Two instruments were used in this work — (a) a double focusing mass spectro- 
meter (2), with a ‘Shaw type’”’ (12) ion source, and (6) a conventional sector 
type mass spectrometer with a heated filament ion source (6). The double 
focusing mass spectrograph, recently converted to a mass spectrometer, employs 
an ion source of the Shaw type in which a tungsten or columbium crucible is 
filled with the sample to be analyzed. The crucible is heated by electron bom- 
bardment from a tungsten filament which surrounds the crucible. Positive ions 
formed inside the crucible at elevated temperatures are accelerated through a 
pin-point hole in the cap covering the crucible to the hole in the shield by the 
same potential which attracts the electrons to the crucible. A further accelerat- 
ing voltage of 3000 to 6000 v. is applied between the shield and the entrance 
slits. 

Analyses were made by measuring singly-charged ions emitted from the oxide 
of the element under study. The oxides were used because they were found to 
give steadier emission than any other compounds investigated. The crucible 
temperature of cesium and rubidium analyses was estimated at between 600 
and 700°C. Strontium analyses required much higher temperatures in the range 
between 1700 and 1800°C. 

Identical analyses were conducted with the 90 degree mass spectrometer 
using a heated filament ion source. Ions of cesium and rubidium were produced 
from the oxides at filament temperatures estimated to be between 600 and 650°C. 
using an optical pyrometer. 


Sample Size 

It was found that satisfactory ion currents could be obtained with very small 
samples of material. For cesium and rubidium these samples amounted to the 
order of 0.1 and 0.02 ugm. respectively, whereas for strontium samples of the 
order of 0.5 to 1 ugm. were required. This high efficiency of ion production 
required that the crucibles and filaments used for heating the samples be free 
from small traces of impurities, which might give rise to ion currents in the mass 
range being studied. In practice each crucible and filament was tested under 
operating conditions for residual ion currents before the addition of a sample. 


Isotope Fractionation 

Experiments were carried out to determine the amount of isotope fractiona- 
tion, if any, which occurred during the thermal emission from the samples. 
Microgram samples of cesium and rubidium were analyzed continuously for six 
to eight hours using both the heated crucible source in the double focusing 
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mass spectrometer and the filament source in the conventional spectrometer. 
In one case using surface ionization from a heated filament, isotopic ratios 
Cs188/Cs!87 were obtained over a period of three days, the total running time 
being about sixteen hours. In this case there was an apparent isotope fractiona- 
tion. The average result obtained over the first two hours differed from that 
of the last two hours by about 2%. In all other cases with both the heated fila- 
ment and Shaw type crucible sources isotope fractionation was less than the 
error of the measurement. In general, the isotope ratios obtained over a six-hour 
period were averaged to give a final result. It seems reasonable to assume that 
under the conditions of our experiments any error in isotopic ratios due to 
isotope fractionation at the ion source is less than + 1%. 


RESULTS 
Cesium 
Table III gives the irradiation data available for the five samples investigated. 
Eight separate determinations were made, three on sample 1, two on sample 2, 


TABLE III 


URANIUM IRRADIATION DATA 


Thermal | | Gain in Xe!6 
neutron flux, |Xe 136/132 | due to neutron 
neutron /cm.?/sec.* ratio capture of Xe! 


Sample Reactor Time in reactor 


No. 


| 
| 
| 


| Oak Ridge | Nov. 17-Dec. 5, 5.9 & 10" 
| 1945, 18 days 
N.R.X. 


| Aug. 11, 1948 
Chalk River | 


N.R.X. | Aug. 9-18, 1950 2X 10" 


| Chalk River 

| N.R.X. 
Chalk River | 
N.R.X. 

| Chalk River 


201 hr. 
Dec. 6-27, 1949 
21 days 


April 3-6, 1951 


56.8 hr. 


1.430 


Ix 16° 1.565 


| 
} | 


*Flux was estimated from data contained in Chalk River Report No. C.R.E. 414—‘Some Data 
on Pile Irradiation” by I. L. Wilson. 


and one each on samples 3, 4, and 5—sample 4 being irradiated in a sheath of 
cadmium. The results are given in Table IV. Column 4 gives the cesium isotope 
ratios obtained with the two mass spectrometers described above. Columns 5 
and 6 give the same ratios corrected for 8 decay and neutron absorption, 
respectively. Finally, column 7 gives the absolute yields of the Cs'**, Cs!®, and 
Cs!*7 chains normalized at 7.43* % for the 133 chain. 
In correcting the yields for neutron absorption only the reaction 

[1] Xe !35 hb n— Xel36 


has been considered. The extent of this reaction can be determined from a 
study of the xenon isotope yields and will depend on the neutron flux density. 


* Values were normalized at 7.43 for 133 mass chain. The previous 6.3% yield for Xe'% (8) 
has been revised in accordance with the recently reported 3.28% I'*' yield (1,9) and a further small 
correction upward has been made to take care of some side reaction which occurs at high flux 
values. Unpublished work in this laboratory. 
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TABLE IV 
YIELDS OF CESIUM ISOTOPES IN THERMAL NEUTRON FISSION OF U2 





= 
Sample | Time after | Cesium | Mass spec.| Ratio corrected] Ratio corrected| Absolute 
run irradiation, | isotope isotope for Cs!87 | for neutron fission 
years ratio decay* capture | yield, t% 
— -|~ = — SS | - |— 
4.60 | le 215 1.110 .110 | 43 
.090 0.995 .045 10 


.000 1.000 .000 72 


215 1.110 .110 43 
. 100 1.005 .070 16 
.000 1.000 .000 3.72 


. 225 1.095 .095 43 
< eee 0.995 .060 16 
.000 1.000 3 5.78 


115 .070 ; 43 
045 .005 a ae 
.000 .000 








. 120 .075 
.400 
.000 


. 105 
.950 
.000 


.135 
. 100 
.000 


. 100 
975 


1.000 

















Absolute yields, % 


\ ee 


av verage and Rey nolds (7 


7.43 
9.18 
6.94 


135 
137 


7 
7 
6. 


This work a Inghram, Hess, 
| 
| 


| 
| 
133 | 
| 


*Correction based on Cs"? half-life of 35 years. 
t Values normaliz sed at 7.43 for 183 mass chain. See footnote on page 428. 
{Sample 2-2 was the came as sample 2-1 with added uranium and reagents to determine Cs'* 
contamination. 
§Cs isotopes extracted from uranium metal irradiated inside a sheath of cadmium. Neutron 
energies >1 ev. 
**Tt is now apparent that Inghram, Hess, and Reynolds (7) were unable to make satisfactory 
correction for neutron capture reactions. Private communication. 
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Insofar as the Xe!** yield is increased by reaction 1, the yield of Xe'*® and, 
therefore, of Cs'** will be decreased by an equivalent amount. Sample 4 irradiat- 
ed in a sheath of cadmium would not be expected to show much burning of 
Xe!* to Xe!46, The correction for neutron adsorption for this sample is, therefore, 
taken as zero. For the other samples the correction is made by comparing the 
ratio Xe!*6/Xe!8? obtained with that found for the cadmium shielded sample. 
In this way the percentage increase in Xe!** due to neutron absorption (reaction 
1) is determined. The Xe!**/Xe!* ratios obtained previously (15) are plotted in 
Fig. 1 as a function of neutron flux. The base level value for the cadmium 
shielded sample is also shown. In the case of samples 1 and 5 the extent of burn- 
ing of Xe!%> was determined directly from xenon isotope measurements. In the 
other cases the neutron flux used had to be estimated and the correction read 
from Fig. 1. 


CADMIUM SHIELDED 


NEUTRON FLUx x 10! 


Fic. 1. Xe®/Xe"? versus neutron flux. 


It might be argued that the above procedure for correcting the cesium yields 
for neutron absorption is not justified as cadmium shielded samples irradiated 
with epithermal and fast neutrons might have a totally different distribution of 
isotopes. However, a large number of samples have been investigated to date 
and only the relative yields Xe'** and Cs'** seem to change appreciably in the 
low flux range, depending on the extent of cadmium shielding.* The relative 
vields of the isotopes Kr**, Kr*4, Kr*, Xe!3!, Xe!8?, Cs!8%, and Cs!*7 remained 
constant within 1 or 2 %. 

Finally, if the extent of reaction 1 is determined from our measurements as 
indicated above then the thermal neutron capture cross section of Xe'*® can be 


*It has been observed recently that the yields of Kr® and Xe do vary to some extent with 
irradiation conditions (unpublished work of R. K. Wanless and H. G. Thode). 
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calculated if the flux is known. The fact that these calculations lead to reasonable 
values of the capture cross section is further evidence that the method used to 
correct the cesium isotope fission yields for neutron absorption is satisfactory. 


Cs!83 Contamination 

Because the Cs!** is used as the point of normalization for the 135 and 137 
chains, the possibility of Cs!** contamination in either or both of the reagents 
and the original uranium metal is extremely important. 

The results of the analyses of the blank samples (See Table IV, Samples 2-1, 
and 2-2) indicate no contamination in either the reagents or the metal used in 
this sample. In the other analyses the same sources of uranium metal and re- 
agents were used, but nevertheless the Cs'**/Cs!*7 ratios indicate a slightly 
higher Cs!** content. This is difficult to account for inasmuch as the entire set of 
results show a uniformity for the Cs!%*/Cs!*7 ratios (as corrected for Cs'*7 
decay, see column 5, Table IV), which would not be expected if the original 
uranium samples contained appreciable amounts of stable Cs!**. 

Further substantiation for the nonexistence of stable contamination is indicat- 
ed by the agreement between the Cs!*" half-life calculated from the mass 
spectrometer data of Table 1V as compared to that obtained by radiochemical 
means. Since the Cs!**/Cs!87 ratio is used in the calculations, any extraneous 
Cs!%8 will be reflected in the value obtained. For example, if the Cs!**/Cs!87 
ratio of sample 5-1 should be altered by a 3% increase in the Cs!** content 
(i.e. from 7.4 to 7.6) then the half-lives resulting would be of the order of 42 to 
46 years. These values differ markedly with the published half-life with which our 
other results agree. From our results we believe the upper limit for stable Cs'** 
contamination to be of the order of 1% of the absolute yield. 

Half-life of Cs'** 

Cs!87 is reported to have a half-life of from 33 to 37 years (9, 11). The Cs'%* 
Cs!57 ratio will, therefore, increase with time. The data of Table IV may be used 
to calculate the half-life of Cs'!*7. The results of these calculations are given 


in Table V. It should be pointed out that some variation in the 133 mass yield 


TABLE V 


Cs'87 HALF-LIFE FROM MASS SPECTROMETRIC ABUNDANCE DATA 


| | ear 
. : | Half-life 
Sample run | Decay time, | Cs 137/133 | Ratio N/No* | Elapsed time, T1/2, 
years ratio | | vears years 


64 . 8230 (a) /(c) .24 31 
5.65 .8163 (a) /(d) .39 31 
4 .9050 (b) /(c) 5.25 35 
.25 .9090 (b) /(d) 5.4 35 

| Average 33 +2 





See Table IV. 
: . C87 Cst37 Vv ' 
*From decay equation —— (a) / —— (c — = ¢~.6931/T, 7, where t = difference in decay 
C133 Csi83 No 


time between samples (a) and (c), etc. 
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may occur owing to neutron capture reactions, etc., which would introduce an 
error in the half-life determinations, particularly for samples irradiated at high 
flux values. (See footnote p. 425.) 


Rubidium 

The mass spectrometer abundance data for the rubidium isotopes: extracted 

from sample 2, Table IV, are given in Table VI. In the case of the cesium isotopes, 
TABLE VI 

RELATIVE YIELDS OF RUBIDIUM ISOTOPES FROM U*55 FISSION 


| Relative Relative | Relative Relative 
abundance abundance abundance abundance Fission 
Rb B¥. Pr. corrected corrected | yield, 
isotope Sample 2. |sample+blank| for normal for Kr® % 
contamination decay 


eehennepscseininsneesniisatinnatosscittais) eitemeencesia = } 


1.00 1.00 | 


| ‘ 0.506 0.374 


* Normalized at 1.3 for mass 85 from mass fission yield curve. 


blank experiments showed little or no contamination with normal Cs!** from the 
uranium or reagents. However, the fission product rubidium was found to con- 
tain normal rubidium as an impurity. The mass spectrometer data for sample 
2 and for sample 2 plus blank in columns 2 and 3 indicate the extent of this 
contamination. The corrected isotopic ratios given in column 4 were calculated 
as follows: 

Let the ratio of the fission yields of Rb** and Rb*’ be 1/a. 

If x is the amount of contamination of Rb*’, 

then 2.6 x* will be the amount of contamination of Rb*. 

Let R = ratio Rb**/Rb* from Fission Product Sample 2. 

Let R! = ratio Rb**/Rb* from Sample 2 + blank 

(blank = equal weight of unirradiated uranium + reagents). 

The measured ratio for F. P. Sample 2 then is given by 

ee EE een 
a+x 


= LEE. nen, 
a + 2x 
Hence 1/a = 0.374. 
The isotopic ratio given in column 5 has been corrected to allow for the decay 
of 9.4-year Kr*. It is assumed that 25% of the 85 mass chain decays through 
_this isomer and 75% through the independent 4.38-hr. Kr*® (16). 


Strontium 
Three fission product samples of strontium were investigated mass spectro- 
metrically, the results are tabulated in Table VII. It was not necessary to run 


*Normal isotopic abundance of Rb**/Rb*" taken as 2.6. 
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TABLE VII 





Age, Mass Ratios 
years Mass spectrometer corrected for 
ratio Sr® decay 





88 
90 


88 
90 


ie 
82 


-_ © 
“1 
_ 


S23 22 8 


88 
90 


a) 


Average 88 
90 


=o 











blanks to determine the extent of normal strontium in the fission product mate- 
rial since of the natural isotopes of strontium only one, Sr*’, is produced in 
fission. Thus an accurate measure of Sr** in the sample from normal contamina- 
tion can be calculated from the amount of Sr** present. The ratios of Sr** to Sr*° 
given in Table VII are corrected for normal contamination in this way. The 
relative yield of 19.9-year Sr®° (10) was obtained by extrapolating to zero time 
to correct for decay. 


DISCUSSION 

The relative fission yields of the rubidium and strontium isotopes are in fair 
agreement with interpolation values from the mass vs. yield curve. These 
results taken together with the yield data reported earlier for the krypton 
isotopes indicate little or no ‘‘fine structure’’ in the 50 neutron shell range for 
U**> slow neutron fission. 

It is just possible that some ‘“‘fine structure’ does exist in this mass range 
owing to closed shell effects, but that it does not show up because of the steepness 
of the curve and the relatively low yields at this point. In this case the ‘‘fine 
structure” should be more evident for U** slow neutron fission since the krypton, 
rubidium, and strontium isotopes then appear near the top of the yield curve. 


Nuclear Shell Structure in Fission 

The high precision mass spectrometer yields obtained in the xenon and 
cesium isotope range (masses 131 to 137) are given in Table VIII, column 3. 
It is clear from Fig. 2 that the mass spectrometer yield values do not fall on a 
smooth mass vs. yield curve. Masses 133 and 134 are for example highly 
favored in the fission process. If we assume that the distribution of both charge 
and mass fundamentally follow simple probability laws, then any ‘‘fine structure” 
is due to other effects superimposed on the elementary fission process. It has 
been suggested earlier that these abnormal yields are due to closed shell effects in 
atomic nuclei, since the mass chains in question have members with 82 neutrons, 
a most stable configuration. If this is true then either nuclides with 82 neutrons 
are favored in fission according to some rule depending on their stability, or 





MASS SPECTROMETER FISSION YIELDS IN NEIGHBORHOOD OF THE 82 NEUTRON SHELL (U?® + mn) 


Mass 
chain 
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TABLE VIII 


Departuref 
Isotope yields, from smooth 
re 


measured % curve 
(%absolute amount) 


429 


Departure 
predicted by 
Glendenin 


mechanism, 
4 








131 
132 


133 
134 
135 
136 
137 


*Values w 
Tt Smooth c 


Xe! 3.28* ais: 
Xe!l# 4.92 ee 
(Xess a i 
\ Csi33 43 1.3 high 
Xels 64 1.7 high 
Cs! .16 — 
Xel%s .10 0.2 low 
Cs!87 5.83 0.6 low 








vere normalized at 3.28 for the 131 mass chain (9, 1). 


urve arbitrarily drawn (solid line Fig. 2). 


(%) 


FISSION YIELDS 


YIELD CURVE PPR 
MASS SPECTROMETER YIELOS 


RADIOCHEMICAL YIELDS 





ee ————E—— EE 


128 130 132 134 136 138 140 142 
MASS 





Fic. 2. Mass spectrometer fission yields for masses 131 to 137 (U*% + n), 


nuclides with loosely held neutrons, one more than a closed shell, emit neutrons 
immediately after fission thereby changing the yields of various mass chains. 
This latter scheme was first suggested by Glendenin. According to his suggested 
mechanism, nuclides with 83 neutrons, or a certain fraction of these, when 
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formed as primary fission products will have sufficient energy to boil out the 
loosely held 83rd neutron. Table I shows the chain branching that would result in 
the mass region studied. 

In order to estimate the extent of these side reactions, ‘‘a”’ to “‘f’’, Table I, 
it is necessary, therefore, to know something about the primary, or independent, 
yields. Although little is known about the primary yields of the individual 
members of a chain, Coryell, Edwards, and Glendenin (4) have proposed a 
theory of equal charge displacement, which seems to fit the available primary 
yield data fairly well. Using this theory Glendenin has estimated the extent of 
reactions “‘a”’ to ‘‘f’’ (Table 1). As a result of this chain branching, the yields of 
the mass chains in question will be high or low according to column 5 (Table 
VIII). The departure of the experimental yields from an arbitrarily drawn solid 
line smooth curve (Fig. 2) is given in column 4 for comparison. Although the 
position of the smooth curve is arbitrary, it is quite evident that masses 133 and 
134 are high and masses 136 and 137 are low as predicted by the Glendenin 
mechanism. 

The Glendenin mechanism requires, however, that yields greater than the 
smooth curve values must be counterbalanced exactly in their totals, by yields 
in other chains, which are below the smooth curve. If the solid curve in Fig. 2 
were raised in the 135-137 mass range, then this would be very nearly true. 
However, the curve is already well above the smooth yield curve of the Pluto- 
nium Project Report and must fall sharply after mass 136. In this connection 
there has been a continual revision upward in the yield of I'*! (9, 1), which in 
turn results in the much higher yields for masses 131 to 137, whose relative 
values have been determined mass spectrometrically. The revised yields are 
given in Table VIII and plotted in Fig. 2. The yields in this mass range are now 
much higher than those given by the smooth curve in the Plutonium Project 
Report. The result is a decided hump on the side of the fission yield curve. 

To explain these facts it was suggested in this laboratory (17) that in addition 
to the Glendenin effect, nuclides with 82 neutrons are favored in the primary 
fission process, that is in the process which occurs before the fission neutrons are 
boiled off. Such an effect would only be appreciable for the masses 133, 134, 
135, and 136 since the increased yield of the primary fragment with 82 neutrons 
would certainly be a function of its displacement from the Bohr—Wheeler 
stability curve. 

One would expect that the Glendenin effect, which depends on the instability 
of nuclides with 83 neutrons, would decrease from mass chains 131 to 138, as the 
member of the chain with 83 neutrons approaches the Bohr-Wheeler stability 
curve (see Table I), whereas the primary fission process effect which depends 
on the extra stability of the 82 neutron shell would increase. 

The ‘‘fine structure” and high total yields found for the mass range in question 
would seem to result from a combination of the two effects discussed above. 
By establishing empirically the functional relationship between the two effects 
and nuclear stability, the results can be accounted for fairly well. 

Since this work (13, 17) was done, Glendenin, Steinberg, Inghram, and Hess 
(5) have reported abnormally high yields in the mass 98-100 range. They point 
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out that the high yield at Mo! suggests a preference for this mass in the fission 
act, perhaps as the complement of a preferred 82-neutron shell in the heavy 
fragment, since mass 100 is the complement of mass 133 and/or 134, both of 
which have high yields. This result then would seem to confirm our suggestion 
of a preference for nuclides with 82 neutrons in the initial fission process. If 
the high yields at masses 133 and 134 were wholly due to a neutron boil-off 
process after fission, the complement mass would be unaffected. 

Experiments are in progress to determine the extent of “fine structure” in 
the mass yield curve in the neighborhood of closed shells for the fission of Pu?*?, 
U3, U88, and Th?*? with thermal and fast neutrons. These experiments should 
give us a crucial test for the Glendenin mechanism, as any fine structure in the 
yield curve due to this mechanism will be shifted, depending upon the number of 
neutrons in the isotope undergoing fission, whereas fine structure due to closed 
shell effects in the primary process would not be shifted. On the other hand with 
very high energy neutron fission, one might expect the preferential formation of 
82-neutron fragments to decrease and the fission process to become less selective. 
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RESONANT SCATTERING OF SLOW NEUTRONS! 


By B. N. BROCKHOUSE 


ABSTRACT 


A method is described for measurement, as a function of neutron energy, of 
the scattering cross section for slow neutrons of materials having high absorption. 
The necessary corrections are discussed. The results are normalized to a vanadium 
scattering cross section of 5.07 barns inferred from the measured ‘total cross 
section, a7 = 5.07 + 5.29(0.025/E(ev.)). Measurements are presented for the 
low energy neutron resonances in cadmium, samarium, gadolinium, dysprosium, 
europium, rhodium, and indium. The single level Breit-Wigner formula is 
accurately obeyed for cadmium and samarium but not for gadolinium in which 
two resonances contribute. The scattering results indicate that the absorption 
cross sections of cadmium and samarium at resonance are higher than the values 
previously accepted. This has been confirmed by total cross section measure- 
ments on cadmium for which the Breit-Wigner parameters are Eo = 0.180 
(+ .003) ev.,f = 0.113 ev., cao = 7750 barns, cp = 4.7 barns, ag = 7 X 10-™" cm. 
The spin of the compound nucleus is unity. For Sm"™* the spin of the compound 
nucleus (0.096 ev. resonance) is i + 3. The scattering cross section of europium 
for neutrons of energy 0.03 ev. to 0.15 ev. is 8 barns, almost independent of the 
energy. The scattering cross section of dysprosium is 100 barns at 0.03 ev. and 
decreases slowly with increasing energy, thus confirming the existence of a level 
at negative neutron energy. The spin of the compound nucleus for the 1.26 ev. 
resonance in Rh! is found to be unity and that for the 1.45 ev. resonance in 
In"5 is found to be four. 


. INTRODUCTION 

Scattering and absorption of neutrons with energies near a resonance energy 
are described in terms of the single level Breit-Wigner dispersion formulas (10). 
These expressions have been compared with experiment individually on dif- 
ferent resonances (1, 27) but at the time this work was begun they had not 
been checked simultaneously and in detail on the same resonance. In this 
paper measurements of the total cross section and the ratio of scattering 
cross section to absorption cross section for the low energy cadmium resonance 
are presented which are in excellent agreement with the Breit-Wigner theory. 
Less detailed examination of several other resonances confirms this agreement. 

Preliminary accounts by Brockhouse, Hurst, and Bloom (7, 5) of part of 
this work have already appeared. Since these publications the apparatus and 
the method of treatment of the results have been improved and the measure- 
ments on cadmium, samarium, and gadolinium described in (7) and (5) have 
been repeated using the new arrangements. 

Tittman and Sheer (31) at Columbia University have measured the ratio 
of the scattering and total cross sections for the 4.87 ev. resonance in gold. 
Their measurements extended over a wider range of energy than was used in 
this work. They obtained good agreement with the Breit-Wigner formulas but, 
because of inadequate resolution, they were unable to measure the total cross 
sections in detail over the same resonance and had to rely on the self-indication 
cross section and resonance integrals, with which their measurements agreed 
well. 


1 Manuscript received December 8, 1952. 
Contribution from Physics Division, Chalk River Laboratories, Atomic Energy of Canada 


Ltd., Chalk River, Ontario, Canada. Issued as A.E.C.L. No. 36. This paper was presented at 
the Quebec meeting of the Royal Society of Canada, June, 1952. 
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The results of Brockhouse, Hurst, and Bloom and of Tittman and Sheer 
must be corrected for the change of neutron energy on scattering. The former 
authors gave no detailed resonance parameters but some of the parameters 
given by Tittman and Sheer are in error, probably by about 10%. The neces- 
sary correction is discussed below. 

The resonances investigated in this work are the low lying resonances which 
determine the properties of the nucleus for thermal neutrons. The detailed 
resonance parameters are therefore of some interest in themselves; in par- 
ticular, the spin of the compound nucleus is the spin of the capture state for 
thermal neutrons, of interest in such experiments as those of Kinsey and 
Bartholomew (16) on capture y-ray spectra. 


EXPERIMENT 


APPARATUS 


The apparatus shown in Fig. 1 was mounted on the arm of the spectrometer 
previously described (13). Monochromatic neutrons, reflected from the 200 
plane of a sodium chloride crystal, entered the scattering chamber through a 
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Fic. 1. Apparatus. The horizontal dimensions are compressed relative to the vertical in 
the ratio 2:1. 


collimating system. Neutrons scattered backwards from the specimen were 
detected in an annular bank of six BF; (enriched in B!°) counters which sur- 
rounded the scattering chamber. Neutrons which were not scattered were 
captured in a boron carbide sink. To reduce background the scattering chamber 
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and collimator were evacuated and the apparatus was shielded with thick 
layers:‘of boron carbide and paraffin. 

The counting rate with a specimen was compared with that from a vanadium 
standard scatterer. The background counting rate was taken with an empty 
sample holder in position and subtracted from the counting rates of the 
specimen and the standard. The background was principally due to neutrons 
leaking through the shield and to counter contamination and was independent 
of whether or not a sample was in position. A thin BF; counter monitored the 
neutron beam and all counting rates were normalized to a common beam 
intensity. 

Collimator and Sink 

The incoming neutron beam, which had an angular divergence of about 
12 min. of arc in the horizontal direction and 45 min. in the vertical direction, 
was stopped down to 3 in. diameter by an aperture in a boron carbide (B.C) 
plate (A). The beam passed through a second 3-in. aperture (B) at a distance 
of 19 in. from the first. Neutrons scattered by the boron carbide were pre- 
vented from reaching the counters by a third boron carbide stop (C) with an 
aperture 3? in. in diameter, through which the 4 in. beam passed without 
interference. 

At the specimen position the fraction of the beam lying outside a ? in. 
diameter circle was less than 1/30,000. Hence the aluminum sample holders 
which were ? in. inside diameter contributed little to the background. 

The neutron beam was caught in a sink of boron carbide contained in a 
thin (~ 0.010 in.) aluminum can. The sink was designed so that no singly 
scattered and few multiply scattered neutrons could reach the counter. 

As a result of these precautions the background from the beam was nearly 
unmeasurable, being about one per hundred thousand countable neutrons in 
the thermal region. 

The boron carbide parts were molded from boron carbide powder using a 
small amount of sodium silicate as binder. A mixture of boron carbide and 
just enough sodium silicate (diluted 50%) to wet it was kneaded into an alumi- 
num mold and then compressed by an arbor press or hammering. The mold 
was then heated slowly to 300° C. and held there for several hours. Parts made 
in this way have fair mechanical strength and good shielding properties and 
can be made in comparatively complex forms. 


THEORY OF THE EXPERIMENT 

In this experiment the scattering by a specimen of unknown scattering 
cross section (¢;) but known total cross section (o7), for which 0, & og “or, 
was compared with that by a standard specimen of vanadium which was almost 
‘“‘thin’’. It was impracticable to use a “‘thin’’ specimen of the unknown scatterer 
and hence it was necessary to take self-absorption into account. For the 
standard vanadium scatterer self-absorption and multiple scattering produced 
small effects which were calculated. 

Because of the cylindrical symmetry of the apparatus and scatterers it was 
possible to compute the scattering for a highly absorbing sample under the 
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condition that ¢,; oq. The approximation used is believed to take good 
account of self-absorption and to take partial account of multiple scattering. 
It is shown in the appendix that, for ¢, Koa, the counting rate from the 
absorbing sample is 
[1] = 6241-c.% \ acy} yl + Cy Sag 
where Cy, Ci, and C2 are constants which were computed from the geometry 
of the apparatus and the angular variation of the counter sensitivity, T is the 
transmission of the specimen, and H(T) is a function which is very slowly 
varying for T < 0.5 and which is also computed from the geometry. In taking 
account of self-absorption in the specimen consideration is given to the fact 
that the energy of the neutron after scattering is, in general, not the same as 
before scattering and hence that the absorption cross section (¢,’) for the 
neutron leaving the specimen after scattering is not equal to the absorption 
cross section (o,) for the neutrons entering the specimen. In the formula the 
quantity (o4’ — o¢a)/oa is averaged over the energy distribution of the neutrons 
after scattering, taking account of the energy dependence of the counter 
sensitivity. 

In an approximation which takes explicit account of multiple scattering, 
the counting rate from the standard scatterer is shown to be 


[2] Ny, = Co Gay G(E) 


where os, is the scattering cross section of vanadium, and C, is another con- 
stant computed from the geometry of the apparatus. G(£) is a slowly varying 
function, computed from the geometry of the apparatus, which takes account 
of self-absorption in the vanadium scatterer and has the value unity when the 
absorption is zero. Combining Equations [1] and [2], the ratio of the scattering 
and absorption cross sections of the absorbing sample is 


N j asl 


~ = o4\ 


(3] oS se Oss = v HT) GE) + G1 % y1 + Cs — 


Ta N Ga 


Equation [3] shows that, in addition to Heath the ratio V/N, and the 
transmission of the sample, it is necessary to know the absorption cross section 
as a function of energy and the distribution-in-energy of the neutrons after 


scattering, in order to obtain os/ea. 
CALIBRATION OF THE APPARATUS 


Cross Sections of Vanadium 

To determine the cross sections of vanadium the total cross section of the 
standard scatterer (which contained 1.784 X 10” atoms/cm.) was measured 
as a function of the neutron energy from 0.02 to 5 ev. and is shown in Fig. 2 
A least squares fit to the measurements gave 


= Ph ins. 
E 


The constant term was assumed to be the scattering cross section (o,»), and 
. 2 1 ° . . . 
the term proportional to E~? the absorption cross section (oay). Scattering 
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Fic. 2. Total cross section of vanadium plotted as a function of the inverse square root 
of the energy. The intercept on the ordinate axis gives the value of the free scattering cross 
section. 


cross sections of 5.0 (12) and 4.7 (8) barns have been previously reported 
without accuracy assignments. The uncertainty in the value 5.07 barns is 
about + 0.15 barns. 

The absorption cross section of 5.29 barns is somewhat higher than the 
value of 4.7 barns obtained by Pomerance (26) from pile oscillator measure- 
ments, the difference being outside the probable errors in the two measure- 
ments. The vanadium used in the standard scatterer was supplied by the 
Vanadium Corp. of America and had a nominal purity of 99.7%, the stated 
impurities being innocuous from the point of view of neutron cross sections. 
The material was examined spectrographically for boron and found to contain 
about 10 p.p.m., insufficient to affect the results. 

In the calculations for this experiment the cross sections os, = 5.07 and 
Say = 5.29 barns were used. Since the absorption cross section affects only 
the function G(£), small errors in oa, do not affect the results appreciably. 
The cross section ratio is, however, directly proportional to ¢,, and hence the 
uncertainty in os, is reflected in the ratio o5/oa. 


Counter Sensitivity Function 

By shielding the counters from the vanadium scatterer with a set of cad- 
mium rings which were removed one at a time, the counter sensitivity was 
determined experimentally as a function of angle. It was necessary to correct 
the measurements for the fact that normally some neutrons were scattered 
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by the walls of the counters and scattering chamber before being counted, 
and part of these were absorbed by the cadmium rings. The correction pro- 
duced about 7% change in the angular distribution. The constants are however 
relatively insensitive to the detailed form of the counter sensitivity function, 
and this correction produced only about 1% change in the final result. The 
counter sensitivity g(@) was found by dividing the corrected experimental 
angular distribution by the scattering 7,(@) expected from the vanadium 
scatterer. Neutrons of two different energies (0.05 and 0.11 ev.) were used 
and the counter sensitivity functions obtained for the two energies were in 
agreement within the experimental error. The function q(@) is shown in Fig. 3. 
It is approximately what would be predicted from the geometry of the appa- 
ratus. 
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Fic. 3. The relative counter sensitivity per unit solid angle g(@) as a function of the angle 
of scattering 0. 


Accuracy of Calibration 

The over-all accuracy of the calibration is estimated to be 5%, of which 
about half is due to the vanadium cross section. The remainder arises from 
errors in the angular sensitivity function and from the approximations used 
in obtaining the constants in Equation [3]. A further error, which varies from 
one specimen to another but is never large, is introduced by the correction for 
the change of neutron energy on scattering and by the effects of binding on 
the scattering by the standard and the specimen. 


TREATMENT OF RESULTS 


Correction for Background and Variations in Incident Flux 

_ Measured counting rates included a contribution from background. This 
background had two components, one proportional to the reactor flux and 
one independent of reactor operation. Counting was carried out for a con- 
venient time (tg) with the appropriate specimen, and the number of signal 
counts Ngo and monitor counts Nps registered. Similar data were taken with 
the standard scatterer (ts, Nso, Nsw) and with an empty cassette for back- 
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ground (tg, Ngo, New). The reactor independent background was measured 
with the reactor off. The signal and monitor counts were each corrected for 
reactor independent background and the signal counts normalized to unit 
monitor count. The ratio N/N,, corrected for background and for variations 







in incident flux, is then 


[4] (=) _SNao —xte__ Neo — sts \ /{Nso —xts _ Nao — xte \ 
N, 


o. New Tt Xtr Nem = Xa) pe re Xt = Neu ri Xty) 





where x is the reactor independent background counting rate in the signal 
counter and X in the monitor counter. 









Correction for Contamination in the Beam 

In addition to neutrons of the desired energy the beam diffracted by a 
crystal from the continuous reactor spectrum contains neutrons of 2nd and 
higher orders which satisfy the Bragg condition, i.e. neutrons of four times, 
nine times, etc., the desired energy. For primary energies of 0.05 to 0.125 ev. 
this order contamination was negligible. Below 0.05 ev. the second order 
content increased rapidly and below 0.02 ev. it became overwhelming. In 
this region the second order content (y) was determined by boron absorption 
as described by Sturm (30). From 0.150 to 0.40 ev. the order content (y) was 
determined by removing the first order neutrons with cadmium and extra- 
polating the residual transmission for various thicknesses of cadmium back to 
zero thickness. Above 0.40 ev. the order content was estimated from the 
residual transmission of samples of various elements at their resonance energies. 

The neutron beam was also contaminated with neutrons diffusely scattered 
from the crystal. These had predominantly energies in the thermal region 
with a mean at about 0.07 ev. The content (6) of diffusely scattered neutrons 
was determined by turning the monochromating crystal slightly out of the 
Bragg position. The ratio (V/N,)o was corrected for contamination by means 


















of the equation 


[5] N (N ‘N»)o = y(N N,)e2 sre 6(N IN ,) 
”» = = —_—-——-- ee 


N, Ly 6 






where (N/N,)2 is the value of the ratio at the second order energy and 
(N/N,)r is its value at the mean energy of the diffusely scattered neutrons. 
The quantities y and 6 in Equation [5] include effects of the energy dependence 
of the counter sensitivity and of the scattering power of the standard scatterer, 
and hence were determined in the scattering apparatus using the vanadium 
scatterer and the annular counters as the detector. In order to keep the cor- 
rection for contamination small, the specimens were made as thin as was con- 









sistent with a small correction for sample transmission. 






Correction for Change of Neutron Energy on Scattering and for Crystal Effects 
To apply formula [3] it is necessary to determine the distribution of neutron 
energies after scattering by the polycrystalline samples and standard. This 
question has been investigated experimentally and compared with the theory 
for Einstein crystals by Brockhouse and Hurst (6). The formula for the 
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Einstein crystal (Reference (6), Equation (3)) was found to be in agreement 
with experiment under conditions applicable in this work and was used to 
obtain the energy distributions over which the averages in Equation [3] were 
carried out. 

The procedure is illustrated by Fig. 4 in which the energy distributions 
computed for cadmium (Einstein temperature 129° K:) for incident neutron 
energies of 0.100 ev., 0.175 ev., and 0.275 ev. are shown. The absorption 
cross section of cadmium is also shown. This is the cross section seen by the 
neutron when entering the specimen. The cross section which is effective for 
the neutron leaving the specimen is the average of the absorption cross section 
over the energy distribution, shown as a dashed line. 
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Fic. 4. The energy distributions of neutrons of initial energy 0.10, 0.175, and 0.275 ev. 
after scattering by cadmium through an angle of 135°, computed according to an Einstein 
crystal model. The thick solid line shows the absorption cross section of cadmium as a function 
of energy. The dashed line shows the effective absorption for the scattered neutrons. 








Because both the mean loss of energy and the width of the distribution 
increase with neutron energy, the correction becomes difficult to handle at 
high neutron energies and this may set an upper limit on the neutron energy 
at which this method of measuring scattering cross sections can be used. 

Under the conditions holding in this experiment, the differential scattering 
cross section for the Einstein crystal is nearly equal to that of the free atom. 
Since the mean angle of scattering in this experiment was 130°-135° the 
measured differential cross section was somewhat lower (31) than 1/42 times 
the free cross section. This reduction, however, was almost exactly compen- 
sated by the increase of the counter sensitivity due to the decrease in mean 
neutron energy so that no corrections were made. The counters had very nearly 

a 1/v dependence on the neutron energy. 

' Vanadium was chosen as a standard scatterer because it gives virtually no 
coherent scattering and hence shows no crystal effects. Crystal effects were 
also of little importance in the specimen scatterers because of the compara- 
tively high temperature, large angle of scattering, and large solid angle. The 

first two result in reduction of the size of the Bragg peaks through the Debye 
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Waller temperature factor. The large solid angle includes many Debye- 
Scherrer rings, the resultant average scattering approaching that from a single 
bound atom. 
Corrections for Resolution and Doppler Effect 

The quantity measured in these experiments is the ratio of the scattering 
and absorption cross sections individually averaged over the Doppler motion, 
the ratio being itself averaged over the instrument resolution. Though o, and 
o, are rapidly changing functions of the energy in the neighborhood of a 
resonance, they behave in much the same way and their ratio is not nearly as 
energy dependent as the quantities themselves. Hence Doppler effect and 
resolution of the apparatus were taken into account only for resonances well 


above the thermal region. 
BREIT-WIGNER FORMULAS 
The absorption and scattering cross sections for neutrons having energies 
near a resonance energy are given by the single level Breit-Wigner dispersion 
formulas (10) 
; ee ve o gact 
[6] Og = gf (a/kko) Ve T,0 [(E —_ Eo) + §r /4]} : 


[7] Os; = Op = TE 7 [Teo /Re + AT, o(E — Eo) ar/kol [(E — E,)’ ts r’/4y-’. 


In these formulas FE is the neutron energy and k = V2ME/h is the neutron 
wave number, I, is the y-ray or absorption width, [, = I'no(Rk/ko) is the 
neutron width, Tf = T,, + Tq is the total width, f is the isotope abundance, 7 is 
the nuclear spin, dz is the potential scattering length of the resonant spin state 
of the resonant isotope, and gp is the so-called potential scattering cross section 
and includes the scattering due to the nonresonant isotopes, the nonresonant 
spin state of the resonant isotope, and the potential scattering by the resonant 
spin state of the resonant isotope. The subscript 0 on any quantity indicates 
its value at resonance. The spin weight factor g has two possible values 

; 1f ‘2 
eI e-store 
corresponding to 
[9] T=i+t+}3 


for the spin of the compound nucleus. 


The ratio o;/aq is 


[10] o s/o = (Rk/Ro) (Tno/Ta) + 4(E — Eo)k ar/Ta + op/oa. 


The first term of this expression is due to the purely resonant scattering, the 

last term to the potential scattering, and the middle term arises from inter- 

ference between the potential and resonant scattering and is zero at resonance. 
If o;/oq4 K 1, as was the case in the present work, then 

[1 1] Tro rs = ko Fal ‘(4g 7 ae 20 a0 Ro”) . 


Where the isotopic assignment of the resonance has been made and the total 
cross section has been determined, the only parameters available to fit the 
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scattering measurements are gp, dp, and g, which can therefore be determined 


unambiguously in most cases. 

The corrected scattering measurements and the theoretical formulas apply 
to free nuclei at rest while measured total cross sections apply to nuclei having 
the motion of atoms in the appropriate solid at room temperature. Parameters 
obtained from total cross section data must be corrected for Doppler effect 
before being used in the above equations. This was done using the methods of 
Lamb (18). In some cases a set of curves given by Wallace (33) was used. 


RESULTS AND DISCUSSION 
Measurements of ¢,/oq over a range of energies in the low energy region 
have been made on specimens of cadmium, samarium oxide, gadolinium oxide, 
europium oxide, and dysprosium oxide. In addition measurements have been 
made on isolated resonances in indium and rhodium. 


Cadmium 

A series of foils of different thicknesses were rolled from Johnson Matthey 
“specpure”’ cadmium metal. The foil having its transmission nearest 0.15 at 
the energy in question was used for the measurements. This thickness was 














0.2 
NEUTRON ENERGY — ev. 


Fic. 5. The ratio ¢3/oq for cadmium as a function of neutron energy. The solid line shows 
the best fitting Breit-Wigner curve, calculated with Ey = 0.180 ev., 'g = 0.113 ev., Tno/Ta 
= 0.00603, cp = 4.7 barns, ag = 7 X 10-" cm., J = 1. 
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about optimum from the point of view of the corrections for specimen trans- 
mission and beam contamination. 

The measured values of the ratio o;/o4 for cadmium, all corrections made, 
are shown in Fig. 5 as circles. An attempt was made to fit the measurements 
to a Breit-Wigner formula using the parameters of Rainwater et al. (27) 
obtained from total cross section measurements. It was not possible to obtain 
a fit within what was believed to be the accuracy of the measurements. The 
total cross section at resonance given by Rainwater et al. was 7200 barns. 
Early work at this laboratory (32) had indicated that the cross section was 
about 10% higher than this value and suggested the desirability of remeasuring 
the total cross section. 

The transmission of the cadmium foils was measured at an energy of 0.175 
ev., close to the resonance energy. The results are shown in Fig. 6, plotted on a 
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Fic. 6. Transmission of cadmium foils at 0.175 ev. plotted against cadmium thickness 
(lower scale). Presence of 2nd order component in beam is indicated. The same points are 
plotted on a more extended scale of cadmium thicknesses (upper scale) to show the effect of 
correction for the 2nd order contamination. " 


logarithmic scale as a function of the foil thickness in gm./cm.? The trans- 
missions fall along two well-defined straight lines corresponding to the primary 
neutron energy of 0.175 ev. and the higher order components in the beam. 
The transmissions of the thinner foils are also shown plotted on a more ex- 
tended scale of thickness together with the values obtained after correction 
for contamination in the beam. The corrected values correspond to a total 
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cross section of 7800 barns at 0.175 ev. Another, independent, set of measure- 
ments gave the same value. 

Similar though less extensive measurements were made at five energies 
below 0.175 ev. Above this energy the correction for contamination became 
large and difficult to handle but at 0.30 and 0.35 ev. the second order com- 
ponent was eliminated by rhodium and indium filters. Extensive sets of 
measurements were made at these two energies. After a small (~1% estimated) 
correction for instrument resolution was applied and the small scattering 
cross section was subtractéd, the measurements were fitted to a Breit-Wigner 
formula with the result 

pcan P 3 Tc a 
[12] 04 = 24.75>/ 0.180/E [(E — 0.180)" + 0.00319]"* x 10-*cm? 


corresponding to ky = 0.180 + .003 ev., 
ao = 7750 + 150 X 10°-*4 cm.?, 
r =0.113 + .002 ev. 


The abundance of the resonant isotope Cd" is 0. 1230 and its spin is 1/2 (25). 
Beeman (2) found that g = 0.75 corresponding to spin 1 for the compound 
nucleus. This is in agreement with the present work. After correcting oo for 
Doppler effect T'no/I'4 is found from Equation [11] to be 0.0060. Using these 
values for the parameters the measurements of Fig. 5 were fitted to Equation 
[10]. The best fitting expression of this form is shown (Fig. 5) and corresponds 


to 
Qe = 7.041.0 X 10°" cm., 
op = 4.7 + 0.5 barns. 


The value of az is close to the geometrical radius of cadmium (6.8 & 107!’ cm.) 
taken to be 1.42 X 107'%4! cm., where A is the atomic weight. The potential 
scattering cross section agrees satisfactorily with the value 5.2 + 0.7 barns 
found by Rainwater et a/. on the high energy wing of the resonance. 

The theoretical expressions for ¢, and a, (Equations [6], [7]) are compared 
with the values obtained from the measurements of ¢,/¢,4 and o7 in Fig. 7. 


Rare Earths 

Measurements were made on specimens of the sesquioxides of samarium, 
gadolinium, dysprosium, and europium (Johnson Matthey ‘‘specpure’’). The 
oxides were heated in air for two hours at 800° C. to remove water and carbon 
dioxide. They were then transferred in a dry box to small aluminum sample 
holders. The specimens and holders were then heated in vacuum at 350° C. 
for several hours. 

The oxides were contained between 0.002 in. aluminum fails. Since only 
the front foil contributed to the scattering, background measurements were 
made using a holder with one foil. 

In computing the energy distribution the Einstein temperature was assumed 
to be 0.025 ev., i.e. equal to room temperature. This is consistent with the 
very scanty specific heat data in the literature (15), and fortunately the energy 
distribution is not very sensitive to the Einstein temperature. The contri- 
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Fic. 7. Absorption and scattering cross sections of cadmium compared with theoretical 
curves. Values of oq shown for 0.225 and 0.25 ev. have little weight because of large corrections 
for contamination and low statistics. 


butions of the rare earth and oxygen ions to the energy distribution were 
computed separately and weighted according to their cross sections. 


Samarium 

The scattering ratio o;/o, was measured for neutron energies of 0.025 to 
0.16 ev. and is shown in Fig. 8, with all corrections made. The resonance at 
0.096 ev. belongs to Sm!'4° (19) which has an abundance of 0.1368 (the 
weighted mean of the two most recent determinations (14, 21)). The spin of 
Sm!*° has been variously reported as 5/2 (24) and 7/2 (3). Experimental (3) 
and theoretical (9) considerations support the value 7/2 which is also in 
better agreement with this experiment. 

The total cross section has been measured by Sturm (30) who found oo = 
15,500 barns and I = 0.074 ev. Using these values it has not been possible 
to obtain a fit within the precision of the measurements, though the spin of the 
compound nucleus is clearly i + }. To obtain a fit oa must be increased and 
this is in agreement with recent work of Sailor (28) who has observed a total 
cross section of 16,800 + 200 barns at 0.094 ev. Assuming this value for oa 
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6.020 


Fic. 8. Ratio of scattering to absorption cross sections for samarium oxide. The solid line 
is the best fitting Breit-Wigner curve, calculated with Ey = 0.096 ev., i = 7/2, J = 4, 
Ty = 0.074 ev., Pno/Ta = 0.00816, op(1/2(Sm:0;)) = 19 barns, ag = 5.0 X 107% cm. 


at room temperature and correcting for Doppler effect T'y»o/T. = 0.00816. If 
I, is assumed to be 0.074 ev. in agreement with Sturm and 7 is taken to be 
7/2, a good fit can be obtained with op = 19 barns and ag = 5 X 107" cm. 


Subtracting 5.6 barns for the oxygen contribution (3/2 ¢, for oxygen) o, for 
samarium is 13.4 barns, somewhat higher than the geometrical cross section 
of about 8 barns, while az is somewhat lower than the geometrical radius. A 
change in Eo or I’, could change these values appreciably however. 


Gadolinium 

Gadolinium has a very large slow neutron cross section which shows re- 
sonance behavior. According to Lapp, Van Horn, and Dempster (19) two 
isotopes contribute to the average cross section over a reactor neutron spec- 
trum, Gd! contributing about 78% and Gd" 22%. The spin of neither isotope 
is known. 

Measurements of ¢,/o,4 for gadolinium oxide for energies from 0.02 to 
0.16 ev., with all corrections made, are shown in Fig. 9. Combining these 
results with the total cross section data of Sturm (30) the scattering and 
absorption cross sections of Fig. 10 are obtained. Sturm, and Brill and Lichten- 
berger (4) were able to fit their total cross section measurements by the single 
level Breit-Wigner formula with a resonance energy of about 0.03 ev. The 
scattering cross section also shows resonance behavior at about this energy 
but does not obey a single level formula. The values expected for the scattering 
cross section at resonance, computed from Sturm’s data for various spins of 
target and compound nuclei on the assumption that the resonance is in Gd", 
are also shown. The disagreement with the measured scattering cross section 
is about as expected from the results of Lapp, Van Horn, and Dempster. 





CANADIAN JOURNAL OF PHYSICS. VOL. 31 

















10 
NEUTRON ENERGY -ev. 


The ratio os/aq for gadolinium oxide. 








CROSS SECTION 








Fic. 10. The scattering cross section of gadolinium computed from os/oq and the total 
cross section of Sturm (30). The bars at 0.03 ev. show the expected scattering cross section at 
this energy (assuming Sturm’s resonance parameters and assigning the resonance to Gd'*) for 
various spins of the target and compound nucleus. 
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Dysprosium 

Measurements of the total cross section of dysprosium by Sturm (30) sug- 
gested that the 1/v law is not obeyed and indicated a large resonance at a 
negative neutron energy of about | ev. Measurements of ¢,/04 for dysprosium 
oxide are shown in Fig. 11. The ratio is seen to be large ranging from 0.1 to 
0.2. Hence the theory of the experiment does not apply accurately and the 
results may be a few per cent in error. When these results are combined with 
Sturm’s total cross section and the oxygen contribution is subtracted, a 
scattering cross section of 100 barns at 0.03 ev. is obtained, which decreases 
almost linearly to 72 barns at 0.10 ev. It should be remarked, however, that 
a rough set of measurements made on the dysprosium sample indicate that 
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Fic. 11. The ratio ¢s/aq for dysprosium oxide. 
8 a d 


while the total cross section in the region of 0.03 ev. is the same as that given 
by Sturm the cross section at higher energies has very nearly a 1/v energy 
dependence. Use of this total cross section gives a scattering cross section of 
100 barns at 0.03 ev. decreasing to 90 barns at 0.10 ev. A large resonance at 
negative energy is again indicated but with the resonance energy considerably 
higher than | ev. (negative). 

Europium 

Results for europium are shown in Fig. 12. Combining these measurements 
with the total cross section data of Sturm (30) and subtracting 5.6 barns for 
the oxygen contribution, the scattering cross section of europium is found to 
be 8 + | barns, almost independent of the neutron energy in the range 0.03 to 
0.15 ev. 

This again is, qualitatively, as expected from the total cross section which 
shows the presence of a resonance at —0.011 ev. in Eu! and at 0.465 ev. in 
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Fic. 12. The ratio os/oq for europium oxide. 


Eu’, The first increases the scattering by Eu'! slightly above the potential 
scattering and the second lowers the scattering by Eu'®* slightly below the 
potential scattering. 
Rhodium and Indium 

Measurements were made on rhodium and indium foils in the vicinity of 
their resonance energies at 1.26 ev. and 1.45 ev. respectively. Because of low 
counting rates it was not possible to determine ¢,/a4 over the entire resonance 
and hence measurements were made only in the immediate vicinity of the 
resonance energies. In these cases the analysis was complicated by the necessity 
of taking into account the instrument resolution. In analyzing the measure- 
ments resonance parameters due to Sailor (28, 29) were employed. For indium 
the cross section at resonance is about 27,000 barns and the total width is 
0.11 ev. For rhodium the cross section at resonance is 4630 barns and the 
width is 0.176 ev. The cross section of indium is in agreement with the work 
of McDaniel (22). The cross section of rhodium is much higher than that given 
by Meijer (23) but is in better agreement with measurements by Lowry and 
Goldhaber (20). The parameters apply to the metal at room temperature and 
had to be corrected for Doppler effect following Lamb (18). 


Rhodium 

The 1.26 ev. level in Rh! is comparatively weak and hence potential 
scattering is important even near resonance. The scattering cross section for 
neutrons of 0.15 ev. energy was measured and found to be 4.1 barns. The 
potential cross section op was then computed from Equation [7] assuming 
de = 7 X10" cm., the geometrical radius of rhodium, and was found to be 
5.3 barns. 

The result of the scattering measurements was ['no/I'g = 0.004 + 0.003 
where the limits represent the estimated over-all uncertainty. This result is 
considerably lower than the value of 0.043 found by Harris et al. (11) from 
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integral scattering and absorption measurements. The spin of Rh! is 1/2 (17) 
and comparison with the total cross section of Sailor by means of Equation [11] 
shows that g is 3/4 and therefore the spin of the compound nucleus in unity. 
This result is consistent with the measurements of Lowry and Goldhaber (20) 
who found gl no/T'a = 0.0029. 

It should be remarked that the counting rate in this experiment was ex- 
tremely low, about 0.4 counts per minute, compared with a background of 
8 counts per minute, which fluctuated by about the amount of the signal. 
Four separate runs were taken and the results were statistically definite but 
the possibility of error cannot be excluded. 


Indium 

The resonance at 1.45 ev. in In™ is so strong that potential scattering has 
no effect. Using the parameters of Sailor (28) and taking 7 = 9/2 (25), the 
scattering measurements showed that I',o/I', = 0.043 + 0.006 and that the 
spin of the compound nucleus is 4. 
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APPENDIX A 
SCATTERING FROM AN ABSORBING SPECIMEN 

A beam of neutrons is incident normally on an infinite plane scatterer of 
thickness d. Per unit area and incident neutron the number of neutrons 
scattered by an element dx at depth x is no, exp( —no7x)dx. It is assumed that 
the scattering is isotropic and that only the absorption cross section is effective 
in reducing the scattered intensity on the way out of the material. If ¢(@) is 
the angular sensitivity of the counters (cylindrically symmetric), the number 
of neutrons scattered from an element dx at an angle @ into an element of solid 

angle 3 sin@ d@ with energy E’, and counted, is 

1 no, q(0) sin @exp{ — nx(o7 + o4(E’) sec 0)} (E’)* F(E’) dE’ do dx 


jf - “F(E") (E')"} dE! 


where F(E’) is the energy distribution and the factor (£’)~? is introduced to 
take account of the energy sensitivity of the counters, which were ‘‘thin’’ and 
hence obeyed the 1/v law. 

Integrating over the sample thickness, the number of neutrons scattered 
at an angle 6 with an energy E’ and counted is 


4[a,/(or + a,’ sec @)|[1 — exp{ — nd(o7r + o,’ sec 0)}] 


X g(0) sin 6 d@ F(E’) (E')* dE’ / f F(E’) (E')? dE’. 
0 
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If the energy distribution is assumed to be that at the mean scattering angle, 
the total number of neutrons counted is 


J 1 F(E’) (E')“aE’ f o;(o7 + o,’sec0) [1 — exp{ — nd(ar + o,'secd)} 


0 


x (0) sind ao/ J F(E’) (E’)} dE’. 
0 


Foro; Koa, ndor > 1, and aq’ — oa K oa, this can be written as Equation [1] 


where 
Co = 1 | q(@) sind (1 + secd)~* dé, 


<3 
ee 
= 


Cy " | g(@) sind (1 + seco)” dé, 
7/7 


ce row '(E’) “Ce wire \—} ~p 
=5 | AS —%« F(R!) (E')7! dE 


0 a 


x | q(@) sind (1 + secd)~” ao/ f F(E’) (E') dE’, 
70 0 


m C 1 oT ‘2 ‘ , Fie —1 
Wer) — E ~~ = - | (T)**8°"9 (8) sind (1 + secd)~ do} , 
= 0 
and T = e-": is the transmission of the specimen. 
APPENDIX B 
SCATTERING FROM THE STANDARD SCATTERER 
Assuming isotropic scattering, the fraction of incident neutrons scattered 
backwards into an angle @ from an element dx and emerging from the sample 
without further interference is 
4 ny Osy dx exp{ —NyXx o7y(1 + sec#)} sind dé 
where the change of neutron energy is neglected. Integrating this over the 
sample thickness yields 
lo,, [1 — exp{ — 2,d,o07,(1 + sec) : 
a [ pl — saeore(l +.0008)}}] sind dé, w/2<0< =. 
2er, (1 + sec@) 


If there were no scattering or absorption on the way out, the number emerging 


[B1] 7, d0 = 


would be 


[B2] I, do = 2 ox } 1 — exp( — en dé. 
2 oT» \ 


The difference between [B1] and [B2] represents the number interfered with 
on the way out, of which a fraction ¢,/o7 is scattered. 

For forward scattering (0 < 6 < 7/2) 
1 os, exp( — Mdror, secd) — exp( — Md ore) 
2 or, 1 — secd 


[B3] 1, dé = sin dé. 
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The total number rescattered is 
o Tr 
U = cus f (Lo = I;) dé. 
OT, V0 


Of these some get out, some are rescattered, and some absorbed. If the 
scatterer is thin, i.e. m»dyory K 1, the neutron density will not vary greatly 
with depth. Assume that the angular distribution and proportion of the multi- 
ply scattered neutrons getting out is the same as for the total scattered 
neutrons. The total number scattered into an angle @ is then 


IT, d@ = (1, + 3(/,/1o) U sin@) dé. 


Hence 


te ep 
[B4] 1 AS) * Fai at 


The number of neutrons per unit area and unit incident neutron which are 


counted is 


[BS] New J 1,(8) g(@) dé. 
0 


For the standard scatterer 3(U/J9) sin@ ~ 0.1 and n,d,o7r, ~ 0.1-0.2, and 
hence /,(@) is nearly proportional to os». Equation [B5] can therefore be 
written in the form of [2] where 


G =, f T .0(0) q(@) dé. 


= 
o 

Tyo is determined at E = © where ogy = O and ogy = ory, and 

on / 

G(E) = T,(8) q(@) a8 / c.. 

Tsy V0 

This form is used in order to emphasize that the results are normalized to a 
particular value of vanadium scattering cross section and will be changed in 


proportion with it. 
The constants Co, Ci, Co, and C, and the functions H(T) and G(£) were 


evaluated by numerical integration. 

The constant C, was also computed in an approximation of the type used 
for the absorbing sample (Appendix A) and found to differ by less than 2% 
from that found in this approximation. The effect of the finite beam diameter 
and specimen size was investigated and found to be negligible. 
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NOTES 


THE ABSORPTION SPECTRA OF FREE RADICALS FORMED IN 
AN ELECTRIC DISCHARGE* 


By P. J. Dynr** 


The study of the absorption spectra of free radicals formed in an electric 
discharge has, up to the present, not been systematic. This paper reports a 
search which has been made, using the technique for the spectra of some © 
diatomic and polyatomic radicals. Absorption spectra of the radicals OH, 
SH, CS, CN, and NH: have been observed in discharges through H.O, H.S, 
CS2, (CN)2, and NH; respectively. Discharges through a number of organic 
compounds were studied in the hope of finding absorption spectra .of such 
radicals as CH2 and CH;. However, no new spectra were observed. 

Oldenberg (9) developed a technique for the observation of OH radicals 
in a discharge through HO. Lewis and White (8) used a similar technique 
to make the first reported observation of the spectrum of the SH radical in 
a discharge through H:S. White (12) studied CN radicals formed in a dis- 
charge through (CN)... Laird, Andrews, and Barrow (7) observed free radicals 
of long life, e.g., CF2 and CS in discharges through a fluorocarbon and CS». 

In the present study a 5 kw. radiofrequency oscillator was used which had 
been used previously in studies of emission spectra (3). A glass discharge tube, 
1.2 meters long and 18 cm. in diameter, was placed along the axis of the 
oscillator coil. Brass end plates were sealed on to the glass tube. These carried 
two concentric extension tubes 5 cm. in diameter in which a set of multiple 
reflecting mirrors of the type described by White (13) could be placed (Fig. 1). 
The separation of the mirrors was 2 meters. The reflecting surfaces were 
normally about 50 cm. away from the luminous region of the discharge so 
that tarnishing was reduced. 





Mo ——VAPOUR IN 
7 - +» TO SPECTROGRAPH 
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TO PUMP ~— 





20 cm. 


Fic. 1. Diagram of the discharge tube showing the position of the set of multiple-reflecting 
mirrors M, and M2. The high frequency coil passed round the middle third of the discharge tube. 


In most experiments the absorption background was obtained with a 
“flash tube’”’ of the type used by Ramsay (10). A well exposed background 
could be obtained with one to five flasheson a Hilger El spectrograph in the 
region 45000-2800. A rotating shutter was used to reduce the intensity of 
emission from the continuously running discharge to a minimum. In some 
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experiments the oscillator was pulsed; the flash tube being synchronized to 
fire at the close of each pulse. 

The main search was carried out using a Hilger El Littrow spectrograph 
with both quartz and glass optics. The spectral region studied with a single 
traversal was \A8000-2500. With a multiple path of eight traversals the 
spectral range \7000~-2800 was studied. A more restricted search was made 
using the second and third orders of a 21 ft. grating between \A7000 and 3100, 
and one traversal of the discharge tube. 

The power level of the oscillator was fixed. A much lower power input 
could be obtained by placing the discharge tube in an auxiliary coil which was 
loosely coupled to the main coil. It was found that the optimum conditions 
for the observation of the OH radical were high pressure and high power 
input. The pressure used was the highest at which a stable discharge could 
be run. These conditions are similar to those found by Oldenberg for OH to 
give the strongest absorption (9); higher pressures could be used in the present 
experiments because of the different characteristics of the high frequency 
discharge as compared with the d-c. discharge used by Oldenberg. The maxi- 
mum pressure varied from compound to compound, but was of the order of 
10 mm. It was found that these conditions also gave the strongest absorption 
with other band systems, and they were therefore used in the search experi- 
ments. 

The observations on diatomic molecules essentially duplicate the results 
of previous workers. The SH band at A3237 was observed in a discharge 
through H.S. Bands of the violet CN system were observed in a discharge 
through cyanogen; the 0-0, 1-1, 2-2, 1-0, and 2-1 bands occur indicating a 
high vibrational temperature of the CN radicals. The bands of CS at 2575 
were observed in a discharge through CS»; this has been reported briefly else- 
where (5). 

From the results reported by Laird, Andrews, and Barrow (7), it was ex- 
pected that the spectrum of CF: should be observed relatively easily in a dis- 
charge through CF. A series of experiments in which the power input, pressure, 
and flow rate of the gas were varied, together with various pulsing techniques, 
all proved negative. The failure of these experiments must be ascribed to the 
use of CFy. Laird, Andrews, and Barrow used a perfluorocarbon (1) Average 
M. Wt. equivalent to C;Fy. which contains a number of CF: groups in each 
molecule. Dr. Barrow also noted that while CF, could be observed in a 10 liter 
bulb, it could not be observed in a discharge tube 4 cm. in diameter, and 
concluded that wall effects played a dominant part in the decay of the species. 
In the apparatus used here, the diameter of the discharge tube being 18 cm., 
wall effects would not be expected to predominate. 

The most interesting result is the observation in absorption of the so-called 
‘“‘e’’ bands of ammonia which have recently been shown to be due to the 
radical NH, (6). Using a single traversal, two lines of the strong NH: band at 
5700 were observed in the second order of the 21 ft. grating in a discharge 
through flowing ammonia at a pressure of ~10 mm. which was pale mauve 
in color. The total length of the luminous portion was about 70 cm. The in- 
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tensity of absorption could be greatly increased by increasing the absorbing 
path using the multiple traversal mirror system. For the spectral region of 
the ‘‘a’’ band (A8000-A4000) a tungsten light bulb gave a very convenient 
continuous background. The aluminum mirrors were not affected by the dis- 
charge under the above mentioned conditions. Satisfactory plates were ob- 
tained with a 5 to 15 min. exposure with a path length of 24 traversals of the 
discharge tube. 

The intensity of absorption in the discharge due to NHz is estimated to 
be greater by a factor of about 3 than that obtained by Herzberg and 
Ramsay (6) in the flash photolysis experiments. Three new bands lying at 
47904, 8160, and 8240 were observed which may, with the more highly de- 
veloped bands in the visible region, provide useful data for the analysis of 
this system. 

In choosing compounds for the search the possibility of observing band 
systems due to CH3, CHe, C3; (Douglas (2)),,CH, HCO (Dyne and Style (4)), 
and HCO, (Style and Ward (11)) was considered. Discharges through CH,, 
CoHa, C2H2, HoCO, CH3CHO, CH;COCHs3, and HCOOH were examined with 
the El spectrograph with one traversal of the discharge tube. Discharges 
through CHy, C2H4, C2H2, and HCOOH were examined with the El spectro- 
graph using 8 or 12 traversals of the discharge tube, and in the second and 
third orders of a 21 ft. grating using one traversal of the discharge tube. A 
discharge through CH, was also examined in the first order of the grating in 
the region \9000-A7000. 

The failure of these last experiments shows that it is not correct to suppose 
that a discharge in which a particular radical appears strongly in emission 
necessarily contains the radical in sufficient concentrations to be observed 
in absorption. For instance CH, NH, and OH bands appear strongly in 
emission in discharges through CH,4, NH3, and H.O. Only OH is observable 
in absorption under my experimental conditions. The a band of ammonia, which 
can be observed in absorption, is, on the other hand, not strikingly strong in 
emission in this discharge. Only in special cases, then, does a discharge contain 
a sufficient concentration of radicals to be observed in absorption. The in- 
tensity of a spectrum in absorption depends on the stationary concentration 
of radicals in the ground state; the intensity in emission depends on the rate 
at which radicals are formed in an excited state and their lifetime in that state. 
These factors are not necessarily correlated. 
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EFFECT OF A LARGE DIELECTRIC CONSTANT ON 
GROUND-WAVE PROPAGATION 


By JAMES R. WAIT AND L. L. CAMPBELL 


The rigorous theory of the propagation of electromagnetic waves around a 
spherical earth of infinite conductivity was first developed by Watson (6). His 
work has been extended to an earth of finite conductivity by Burrows (2), 
Norton (3), Bremmer (1), and others. Theoretically therefore, solutions are 
available for any value of the earth’s constants (dielectric constant and conduc- 
tivity) for either vertically polarized or horizontally polarized waves. 

Usually in the practical computations of these solutions the relative dielectric 
constant for the ground is assumed to be in the range from 2 to 15. This is quite 
justified in most cases for frequencies above 1 Mc. per sec. as indicated by 
typical measured values, given by Smith Rose (4). At lower frequencies how- 
ever the dielectric constant can be much greater. For example at 200 kc. per 
sec. Smith Rose has found the dielectric constant to be greater than 100 for 
damp soil. When the ground conductivity is poor this can be expected to have a 
pronounced effect on the propagation. Calculations show that this is the case. 

The electric field strength in millivolts per meter for a point on the earth at a 
distance d miles from the transmitter is shown plotted in Fig. 1. The antenna 
radiates vertically polarized waves at 200 kc. per sec. and the field intensity 
at 1 mile for a perfectly conducting earth would be 100 millivolts per meter. 
For this example the ground conductivity is taken as 0.1 X 10~" in electro- 
magnetic units. Various values of the dielectric constant ¢ of the ground are 
shown. For values of ¢ from 2 to 200 Norton’s formulas can be used directly for 
distances up to about 85 miles which is the limit of the flat earth approximation 
at this frequency. For distances up to 2000 miles Norton gives formulas that 
can be applied to the case when ¢ is not greater than 15. The complete curve 
of « = 200 is calculated from Bremmer’s expansion formula given on page 45 
in his book. The remaining portions of the curves for e equal to 30 and 80 are 
interpolated. 

The curves indicate that an effective dielectric constant of the order of 100 
can at great distances from the transmitter give rise to field strengths up to 
100 times the field strength for the case when the dielectric constant is of the 
order of 3. Charts (5) similar to this one have been prepared for other typical 
values of the ground conductivity in the frequency range from 15 to 500 kc. 
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Fic. 1. The field strength of a vertical radiation on the surface of the earth showing the 


effect of changes of the dielectric constant of the ground. 


From this work it appears that the dielectric constant as well as the conduc- 
tivity has a pronounced effect on ground-wave propagation at low frequencies, 
particularly when the ground conductivity is less than 1.0 XK 107" e.m.u. in 
value. 
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